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Abstract 

We compute the 0{as) virtual QCD corrections to the 7 * —> qqg matrix element arising from the 
interference of the two-loop with the tree-level amplitude and from the self-interference of the one-loop 
amplitude. The calculation is performed by reducing all loop integrals appearing in the two-loop ampli¬ 
tude to a small set of known master integrals. Infrared and ultraviolet divergences are both regularized 
using conventional dimensional regularization, and the ultraviolet renormalization is performed in the MS 
scheme. The infrared pole structure of the matrix elements agrees with the prediction made by Catani 
using an infrared factorization formula. The analytic result for the finite terms of both matrix elements 
is expressed in terms of one- and two-dimensional harmonic polylogarithms. 



1 Introduction 


Among jet observables, the three-jet production rate in electron-positron annihilation plays an outstanding 
role. The initial experimental observation of three-jet events at PETRA |Q, in agreement with the theoretical 
prediction provided first evidence for the gluon, and thus strong support for the theory of Quantum 
Chromodynamics (QCD). Subsequently the three-jet rate and related event shape observables were used for 
the precise determination of the QCD coupling constant as (see [|| for a review). Especially at LEP, three-jet 
observables were measured to a very high precision and the error on the extraction of as from these data is 
dominated by the uncertainty inherent in the theoretical next-to-leading order (NLO) calculation of the 
jet observables. The planned TESLA linear e+e“ collider will allow precision QCD studies at even higher 
energies than at LEP. Given the projected luminosity of TESLA, one again expects the experimental errors 
to be well below the uncertainty of the NLO calculation. 

Related to e+e“ —> 3 jets by crossing symmetry are (2-|-l)-jet production in deep inelastic ep scattering and 
vector-boson-plus-jet production at hadron colliders. The experimental data from HERA on ep —> (2 + 1) jets 
and related event shape observables have already reached a level of precision demanding predictions beyond 
the present NLO accuracy; a further improvement on these data is expected soon from the HERA high 
luminosity programme. Similarly, vector-boson production at large transverse momentum is a classic test of 
QCD in hadron-hadron collisions and demands the theoretical prediction to be as precise as possible. In this 
case, it is also an important background in searches for new physics at the Tevatron and the LHC. 

Besides its phenomenological importance, the three-jet rate has also served as a theoretical testing ground 
for the development of new techniques for higher order calculations in QCD: both the subtraction Q and the 
phase-space slicing [|| methods for the extraction of infrared singularities from NLO real radiation processes 
were developed in the context of the first three-jet calculations. The systematic formulation of phase-space 
slicing as well as the dipole subtraction [|| method were also first demonstrated for three-jet observables, 
before being applied to other processes. It is very likely that similar techniques at higher orders will first 
be developed in the context of jet production in e“'"e“ annihilation, which in contrast to hadron-hadron 
collisions or electron-proton scattering does not pose the additional difficulty of the regularization of initial 
state singularities. 

The calculation of next-to-next-to-leading order (NNLO), i.e. 0{al), corrections to the three-jet rate in 
e“'"e“ annihilation has been considered as a highly important project for a long time ||l^ . In terms of matrix 
elements, it requires the computation of three contributions: the tree level 7 * ^ 5 partons amplitudes 
the one-loop corrections to the 7 * ^ 4 partons amplitudes El- 0 ^ and the two-loop (as well as the one-loop 
times one-loop) corrections to the 7 * ^ 3 partons matrix elements. While the former two contributions have 
been known for some time already, the two-loop amplitudes have presented an obstacle that prevented further 
progress on this calculation up to now. 

This calculation has now become tractable owing to various technical developments over the last two 
years. In particular, the systematic application of integration-by-parts |[l^-|^ and Lorentz invariance [ pl[| 
identities allowed the large number of Feynman integrals appearing in two-loop four-point matrix elements 
to be reduced to a small number of so-called master integrals. The use of these techniques already allowed 
the calculation of two-loop QED and QCD corrections to many 2 —> 2 scattering processes with massless on- 


shell external particles 
with all legs on-shell 



-|27|] , which require master integrals corresponding to massless four-point functions 
The master integrals relevant in the context of the present work are massless 
four-point functions with three legs on-shell and one leg off-shell. The complete set of these integrals was 
computed in Q , earlier partial results had been presented in |^, ^ . 

The master integrals in are expressed in terms of two-dimensional harmonic polylogarithms (2dHPLs). 
The 2dHPLs are an extension of the harmonic polylogarithms (HPLs) of All HPLs and 2dHPLs that 
appear in the divergent parts of the planar master integrals have weight < 3 and can be related to the 
more commonly known Nielsen generalized polylogarithms |^,|^| of suitable arguments. The functions of 
weight 4 appearing in the finite parts of the master integrals can all be represented, by the very definition, as 
one-dimensional integrals over 2dHPLs of weight 3, hence of Nielsen’s generalized polylogarithms of suitable 
arguments according to the above remark. A table with all relations is included in the appendix of 
Numerical routines providing an evaluation of the HPLs and 2dHPLs are available. 

In this paper, we present the 0{ag) corrections to the 7 * ^ qqg matrix element. At this order, two combi- 
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nations of amplitudes contribute: the interference of two-loop and tree amplitudes and the self-interference of 
the one-loop amplitude. We work in conventional dimensional regularization with d = 4 —2e space- 

time dimensions, where all external particles are d-dimensional. Ultraviolet renormalization is performed in 
the MS scheme. The infrared pole structure of the two-loop corrections to the 7 * —> qqg matrix element was 
predicted by Catani [Q, using an infrared factorization formula. We confirm Catani’s prediction with our 
explicit calculation, and we use the formalism introduced in ||4^ to present the infrared poles and the finite 
parts of the 7 * —> qqg matrix elements in a compact form. 

The paper is structured as follows. In Section H, we define the notation and kinematics used in the paper. 
Section ^ explains the method we used to reduce the two-loop diagrams to master integrals. The result for 
the two-loop QCD contribution to the 7 * —*■ qqg matrix element, decomposed into infrared-divergent and 
infrared-finite parts according to the prescription derived in |Q|, is given in Section Finally, Section ^ 
contains the conclusions and an outlook on future steps needed for the completion of a full NNLO calculation 
of three-jet production in e“'"e“ annihilation. 


2 Notation 


We consider the decay of a virtual photon into a quark-antiquark-gluon system: 

7*('3') —> q{pi) + q{p2) + g{p3) ■ 

The kinematics of this process is fully described by the invariants 

S 12 = (pi +P 2 )^ , Si3 = (pi +P3)^ , S23 = {P2+P3f, 


which fulfil 


q^ — S12 + S13 -I- S23 = 5123 • 

It is convenient to define the dimensionless invariants 

X = S12/S123 , V = S13/S123 , z = S23/S123 , 


( 2 . 1 ) 


( 2 . 2 ) 

(2.3) 

(2.4) 


with X + y + z = 1 . 

Our calculation is performed in conventional dimensional regularization with d = 4 — 2e, and all 

external particle states are taken to be d-dimensional. Renormalization of ultraviolet divergences is performed 
in the MS scheme. The renormalized amplitude can be written as 


\M) = VI TTaeq\/47Ta's 


+ (^) |A(<«) + lA)®) + 0(0?) 


(2.5) 


where a denotes the electromagnetic coupling constant, e, the quark charge, as the QCD coupling constant 
at the renormalization scale p, and the \M.^V are the i-loop contributions to the renormalized amplitude. 
They are vectors in colour space. 

The squared amplitude, summed over spins, colours and quark flavours, is denoted by 

{M\M) = ^ |AI(7* ^ qqg)? ='^{x,y,z). ( 2 . 6 ) 

The perturbative expansion of T {x, y, z) at renormalization scale V = q^ = S 123 reads: 


T{x,y,z) = leVa^elasiq"^) 


T^‘^\x,y,z) -f '^^^Hx,y,z) 


^) ri«>(x.„,.) + 0(o?(,2)) 


(2.7) 
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where 


r(2) (x, y, z) = = 4F(1 - e) 


(l-e)(^ + - 


2{1 - y - z) - 2eyz 
yz 


( 2 . 8 ) 


r('‘)(x,2/,z) = , (2.9) 

T^^\x,y,z) = , ( 2 . 10 ) 

where F = N'^ — 1, with N the number of colours. T^^\x^y, z) was first derived in |0,H; we quote an 


( 2 . 11 ) 


explicit expression for it in Section 4.1. In the following, we present the contribution to (a:, y, z) from the 
interference of two-loop and tree diagrams 

^(6.[2x0])(2;,y,z) = (A4(°)|A4(2)) -f (M(2)|_y\^(0)^ ^ 
as well as the one-loop self-interference 

7-(6.[ixi])(^^y^^) = (M(i)|7W(i)) . (2.12) 

At the same order in as, one finds also a contribution to three-jet final states from the self-interference of the 
7 * —> ggg amplitude. The matrix element for this process does not contain infrared or ultraviolet divergences; 
it was computed long ago and can be found in 145, ^. 

For the remainder of this paper we will set the renormalization scale =9^- The full scale dependence 
of the perturbative expansion is given by 


T{x,y,z) = l6Tr^a^elasin‘^)lT^^\x,y, z) 

(1 ^ 

{x, y, z) + {x, y, z) In 


as(y^) 


27r 


27T 




(x, y, z) + (^2boT^'^'> {x, y, z) + (x, y, z)^ In 


+62T(2)(x,y,z)ln 


2 / M 


+ C’(as) 


(2.13) 


3 Method 


The Feynman diagrams contributing to the i-loop amplitude {i = 0,1,2) were all generated using 

QGRAF 1^^. There are two diagrams at tree-level, 13 diagrams at one loop and 229 diagrams at two loops. 
We then project by and by and perform the summation over colours and spins 

using the computer algebra programs MAPLE , FORM2 and FORM3 [Q . When summing over the 
polarizations of the external gluon and off-shell photon, we use the Feynman gauge: 

^ efer = (3.1) 

spins 

This is valid because the gluon always couples to a conserved fermionic current, which selects only the physical 
degrees of polarization. The use of an axial polarization sum to project out the transverse polarizations (as 
applied in ||^,|^]) is therefore not needed. 

The one-loop self-interference contribution is computed by reducing all tensorial loop integrals 

according to the standard Passarino-Veltman procedure to scalar one-loop two-point, three-point and 
four-point integrals. It has been known for a long time that those three-point integrals can be further re¬ 
duced to linear combinations of two-point integrals using integration-by-parts identities. After this reduction, 
7'(6,[ixi]) jg expressed as a bilinear combination of only two integrals: the one-loop box and the one-loop 
bubble, which are listed in Appendix 

The computation of is by far less straightforward. It is explained in detail in the following two 

subsections. 
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3.1 Master Integrals 


The integrals appearing in the individual two-loop diagrams contain up to seven propagators in the denom¬ 
inator, and up to four irreducible scalar products in the numerator (i.e. scalar products which can not be 
expressed as linear combinations of the occurring propagators). Each integral is classified by the number of 
different denominators t, the total number of denominators r and the number of irreducible scalar products 
s. The set of t different denominators defines the topology of the diagram. 


Using the reduction procedure described in Section 3.2 below, all of the two-loop Feynman diagrams can 


be reduced to a basis set of master integrals. Owing to the presence of the extra scale, there are considerably 
more master topologies than in the on-shell case. Altogether there are 14 planar topologies and 5 non-planar 
topologies requiring a total of 24 master integrals, as five topologies possess two master integrals (see below). 
For each topology, we take one of the master integrals (or the master integral, when it is unique) to be equal 
to the full scalar amplitude with the first power on all propagators. The simpler master integrals are the 
single scale integrals [p2|-^, which can be written in terms of T functions, 


Sunrise(si2) = 


Glass(si 2 ) = 


Darti(si2) = 


Pl2 




as well as the more complicated 


Xtrii(si2) — 



The two-scale master integrals can be written in terms of T functions, 

P 3 

^ Pl2/^~'''^k^^~^Pl23 


Tglass(si 2 ,si 23 ) 




or as generalized polylogarithms or one-dimensional harmonic polylogarithms (see e.g. 

—*—Pl2 


Dart2(si2,si23) = 


Dart2(si23,S23) = 


Plane(si2, S123) = 


Xtri 2 (si 23 ,Sl 2 ) = 


P123, 



P123. 



P123 



Pl 23 y/\ 





-P 3 

-Pi 

-P23 

P12 

P 3 

P12 

P 3 
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The three-scale master integrals can be written in terms of two-dimensional harmonic polylogarithms 
There are the planar graphs, 


Aboxi(s23) Sia, S123) — 


Ab 0 X 2 (s 23 ,Sl 3 ,Sl 23 ) = 


Cb0Xi(s23, S13, S123) = 


Cb0X2(s23j Sis, S123) — 


Tb0Xi(s23, Si3, S123) = 


Pb 0 Xi(s 23 ,Si 3 ,Si 23 ) = 


Bbox(s 23 ,Si 3 ,Si 23 ) = 


and the non-planar graphs, 


Eboxi(si 2 , S13, S123) = 


Xbmoi(si 2 , Si3, S123) — 


Xbmii(si 2 , Si3, S123) — 



Pl 23 -*~ 

Pi~*- 


P123-*- 

Pi~*- 


P123 -o 

Pl~*— 


-P2 

-P 3 



■*-P2 

-*~P 3 


-P2 

-P 3 



It turns out that for the Cbox 2 , Pbox, Ebox, Xbmo and Xbmi topologies, a second master integral is required 
for expressing all the occurring integrals. One is free to choose the precise form of these additional master 
integrals in terms of additional powers of a propagator or as a tensor integral, provided the second master 
integral is not directly related by LI and IBP identities to the first one. As in |^] we choose as second master 
integral 


Cb0X2A(s23j Si3, S123) 
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Pb0X2(s23,Sl3,Sl23) 


Eb0X2(si2,Sl3,Sl23) 


Xbm02(si2,Si3,Si23) 


Xbmi 2 (si 2 ,Sl 3 ,Sl 23 ) 


Pl23- 

Pl- 

Pl23- 


( 2 ) 


-P2 

-P3 

-P2 


P3-* - V \ ». Pi 

P2 



Pi 

P2 

-Pi 


3.2 Reduction to master integrals 

The reduction of the two-loop diagrams to master integrals using integration-by-parts (IBP) ||l^-|^ and 
Lorentz invariance (LI) [0 identities was performed with two independent programs, which allowed a com¬ 
parison diagram by diagram. 

One program was based on the reduction algorithm already described in . In this program, reduction 
equations are derived by applying the differentiation operators of the IBP and LI identities to the following 
integrands: t = 7,r = 7, s < 4; 3 < t < 6,t < r < 7,s < 3 (note that it is not necessary to include 
s = 4 integrands with t < 7). As a result, one obtains a linear system of equations, which contains, besides 
integrals of the class t,r,s and simpler integrals, also more complicated integrals of the classes <,r -|- l,s 
and t, r -I- 1 , s -I- 1. The number of equations does, however, exceed the number of integrals, thus yielding 
an apparently overconstrained system of equations. By solving this system of equations using the computer 
algebra programs MAPLE [Q and FORM2 one expresses integrals of a given class t,r,s as linear 
combinations of simpler integrals. The irreducible integrals, which are left over after this procedure, are the 
master integrals. Compared to the algorithm used in one major modification was made. For integrals 
corresponding to those topologies which do not possess an independent master integral (i.e. can be fully 
reduced to integrals with one less propagator by the IBP and LI identities), it is straightforward to obtain a 
symbolic solution of the IBP and LI identities. This solution takes the form of a single equation decreasing 
the power on one or more of the propagators (see | |^ for an example). These symbolic reduction equations 
for all non-master topologies were derived automatically using a FORM2 program. The use of these 
symbolic identities did lead to a considerable speed-up of the reduction of the amplitudes of such diagrams 
to master integrals. 

The second approach derived the reduction equations using the IBP and LI identities for three auxiliary 
diagrams — one planar and two non-planar. These auxiliary diagrams are obtained by mapping the nine 
possible dot-products involving the three independent external momenta (pi, p 2 and ps) and the two loop 
momenta as propagators. Scalar products in the numerator are replaced by propagators raised to negative 
powers. For illustration, the planar auxiliary diagram is shown in Fig. 1, with each propagator i raised to the 
arbitrary power Vi. An advantage of this method is that the 10 IBP and 3 LI identities are written in terms 
of arbitrary and are valid for all subtopologies of the auxiliary diagram {vi 0 ) and all tensor integrals 
(vi < 0). Typically we encounter subtopologies with up to 7 propagators {t = 7) and tensor integrals up 
to fourth rank (s = 4). For a given value of t, and assignment of propagators, we generate equations using 
tensor integrals with s < 4 as seeds for the generic IBP and LI equations. These seed integrals are those that 
appear in the actual calculation of Feynman diagrams. For example, for t = 7, there are 15 seed integrals 
yielding 195 equations. The equations are stored and solved one by one according to a ranking similar to that 
proposed by Laporta Q using MAPLE [Q and FORMS [Q. At each stage, more complicated integrals are 
decomposed in terms of simpler integrals of the same topology and pinched integrals. As in the first method, 
the remaining irreducible integrals are master integrals. Of course, the system of equations contains many 
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P3 





Figure 1: The planar and two non-planar auxiliary diagrams. All nine possible dot-products involving the 
loop momenta and the external particle momenta are mapped onto the nine propagators. 


more complicated integrals (with, for example, r > 7 or s > 4), which cannot be reduced to master integrals 
with this approach. However, they are not needed for the evaluation of the Feynman diagrams and can be 
eliminated from the system of equations. 

In both methods, one observes that on the face of it, the planar and non-planar graphs appear to be of 
equal difficulty — requiring the same depth of tensor reduction with the same numbers of seed integrals and 
the same number of IBP and LI equations. However the complexity of the non-planar graphs is first revealed 
in the number of terms in each equation and then in the work done to eliminate the unwanted integrals from 
the system of equations. 

Both methods have their advantages. The major advantage of the first method is that the reduction 
equations for each topology containing master integrals have to be evaluated only once, while the second 
method requires their repeated evaluation each time a topology appears as subtopology of the auxiliary 
diagram. As a consequence, the first method needed less computer time for the derivation of all required 
reduction equations. An advantage of the second method is that it needs less human intervention, since 
mapping of subtopologies to master topologies is not necessary, and the generalization to more complicated 
graphs with additional energy scales requires only the modification of the 13 IBP and LI equations. The 
solution of the system and identification of master integrals is only limited by computer resources. 

3.3 Expansion of master integrals 

The two-loop master integrals relevant to the 7 * ^ qqg matrix element are two-loop four-point functions with 
one leg off-shell. These functions were all computed in in the framework of dimensional regularization with 
d = 4—2e space-time dimensions. The results of take the form of a Laurent series in e, starting at e~'^, with 
coefficients containing one- and two-dimensional harmonic polylogarithms |^ , which are a generalization of 
Nielsen’s polylogarithms All master integrals in ||3^ were given for one particular configuration of the 
external momenta. They were expressed in a form where the argument of the 2dHPLs was always y, while ^ 
appeared in the index vector of the 2dHPLs; y and z appeared also as argument of the HPLs. 

Each master integral can occur in six kinematic configurations (corresponding to the permutations of 
(pi,P 2 ,P 3 ))- To avoid hidden zeros (arising from cancellations occurring in the combinations of HPLs and 
2dHPLs with different arguments and different variables in the index vector), we express the master integrals 
for all kinematic configurations in a unique form, which is the same as in [p4| : the argument of the 2dHPLs 
is always y, the variable in their index vector is z, which appears also as argument of the HPLs. To obtain 
this unique form, we apply the following relations to the 2dHPLs, which are all derived following the same 
lines as the formulae for the interchange of arguments discussed in the appendix of p^ : 

1. Interchange of arguments: G(m(j/); z) ^ G(m(z); ?/) -I- H(n; z) . 

2. Replacement of argument variable: G(m(z); 1 — y — z) ^ G(m(z); j/) - 1 - H(n; z). 

3. Replacement of index variable: G(to(1 — y — z);y) ^ G(m(z); y) -|- H( 77 ; z). 
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The relations required for the two remaining combinations G(m(?/); 1 — y — z) and G(m(l — y — z); z) are 
obtained by applying the interchange of argument relations followed by the appropriate replacement relations. 

The master integrals in were derived in the kinematical situation of a (space-like) 1 —> 3 decay, which 
corresponds to the 7 * ^ qqg, such that the only analytic continuation of them required here is the expansion 
of the overall factor in the time-like region 

Re(-l)-2" = 1 - 2Tr^e^ + 2/37r^e'‘ -h 0{e^) . (3.2) 

The analytic continuation of the master integrals to other kinematical regions is discussed in the appendix 
of ||. 


3.4 Ultraviolet renormalization 


The renormalization of the matrix element is carried out by replacing the bare coupling ao with the renor¬ 
malized coupling as = as(/i^), evaluated at the renormalization scale 




e V27r/ 



(3.3) 


where 


Ss = ( 47 r)'^e with Euler constant 7 = 0.5772 ... 


and /Tg is the mass parameter introduced in dimensional regularization M,^,43| to maintain a dimensionless 


coupling in the bare QGD Lagrangian density; /3o and f3i are the first two coefficients of the QCD /3-function: 


/3o — 


IICa - ^TrNf 


6 


Pi = 


17C\ - IQCaTrNf - 6CfTrNf 
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with the QGD colour factors 


Ca =N, Cf = 


- 1 
2N ’ 


Tr=1- 


(3.4) 


(3.5) 


We denote the i-loop contribution to the unrenormalized amplitudes by using the same normal¬ 

ization as for the decomposition of the renormalized amplitude (2.5). The renormalized amplitudes are then 
obtained as 


4 The matrix element 

We further decompose the renormalized one- and two-loop contributions to as a sum of two terms 

7 '(6.[*xi])(a;^y^2;) = Voles^^^^\x,y, z) + Tinite^^^^\x,y, z). (4.1) 

Voles contains infrared singularities that will be analytically cancelled by those occurring in radiative processes 
of the same order (ultraviolet divergences are removed by renormalization). Tinite is the renormalized 
remainder, which is finite as e —> 0. In this section we first give explicit expressions for the infrared pole 
structure using the procedure advocated by Catani Q] and then give the analytic results for the finite 
remainders. For simpli city we set the renormalization scale ^ = S 123 and restore the renormalization scale 
dependence using Eq. dUi). 












4.1 Infrared factorization 

Catan i has shown how to organize the infrared pole structure of the two-loop contributions renormalized 
in the MS scheme in terms of the tree and renormalized one-loop amplitudes, and respectively, 

as 


= 23? 




r(l-e) ve 


and 


=3? 


2(7W(^)|7(^)(e)|At(°)) - (M(°)|7^^)1'(e)7(^^(e)|At(°)) 


where the constant K is 


K = {— -—\ Ca- —TrNf. 
' 18 6 y 9 


(4.2) 


(4.3) 


(4.4) 


It should be noted that, in this prescription, part of the finite terms in are accounted for by the 

0{e^) expansion of Voles^^^^\ 

For this particular process, there is only one colour structure present at tree level which, in terms of the 
gluon colour a and the quark and antiquark colours i and j, is simply Tjt. Adding higher loops does not 
introduce additional colour structures, and the amplitudes are therefore vectors in a one-dimensional space. 
Similarly, the infrared singularity operator 7*-^^ (e) is a 1 x 1 matrix in the colour space and is given by 


7«(e) = - 


an 


2 r(l - e) 


N 


3 l3o \ , ^ 1/1 


2 e 


where (since we have set = S123) 


S., = I - 


S123 


(4.5) 


(4.6) 


Note that on expanding S^-, imaginary parts are generated, the sign of which is fixed by the small imaginary 
part +i0 of Sij. Other combinations such as (Al^°^|7^^^1(e) are obtained by using the hermitian conjugate 
operator 7^^^1(e), where the only p rac tical change is that the sign of the imaginary part of S is reversed. The 
origin of the various terms in Eq. (O) is straightforward. Each parton pair ij in the event forms a radiating 


antenna of scale Sy. Terms proportional to Sij are cancelled by real radiation emitted from leg i and absorbed 
by leg j. The soft singularities C>(l/e^) are independent of the identity of the participating partons and are 
universal. However, the collinear singularities depend on the identities of the participating partons. For each 
quark we find a contribution of 3/(4e) and for each gluon we find a contribution of /3o/(2e) coming from the 
integral over the collinear splitti ng fu nction. 

Finally, the last term of Eq. (|4.2[) that involves produces only a single pole in e and is given by 


(M(0)|H'(2)(g)|_y^(0)) ^ 


4er(l-e) 




(4.7) 


where the constant 77^^^ is renormalization-scheme-dependent. As with the single pole parts of the 

process-dependent can be constructed by counting the number of radiating partons present in the event. 
In our case, there is a quark-antiquark pair and a gluon present in the final state, so that 


77^2) = 277^) -f 77f) 


(4.8) 
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where in the MS scheme 
'1 


7j(2) 

9 


^( 2 ) 

9 




7, 409 

-Cs H- 

4^^ 864 


5 

IItt^ ^ 

) -H — iV| -b ( 

^ TT^ 

89 \ 

12 

144 / 

1 27 \ 

V 72 

108 / 




NNf - 


Np 

W’ 


(4.9) 




1 41 TT^ 

~4^^ ~ 108 ~ % 


, 3 . 3 

' 2*’^ 32 


1 
8 J 


7r‘ 25 \ {N^ - 1)Nf 


,48 216 


N 


(4.10) 


so that 


F( 2 ) = (4C3 + 


589 llTT^ 


432 


72 




41 TT^ \ 


2^^ 54 48 / 


+ 




1 


19 TT^ 


18 36 


1 


-777 + 7777 NNp + "TT - 777 Kf + 73/V^- 


54 24 / TV 


Np 


27 


(4.11) 


The factors Hq^'^ and are directly related to those found in gluon-gluon scattering psl], quark-quark 


scattering |23| and quark-gluon scattering 
factor 1^. 


(which each involve four partons) as well as in the quark form 


59| . We also note that (on purely dimensional grounds) one might expect terms of the type Sf^ 
to be present in Of course such terms are 1 + C>(e) and therefore leave the pole part unchanged and 

only modify the finite remainder. At present it is not known how to systematically include these effects. 

The renormalized interference of tree and one-loop amplitudes also appears in Eq. ( |4.2| ). This can be 
written to all orders in e using the relation 


where 


26 

(M(o).-|M(i).-) = y (^Nf,iy, z) + ^My, z) + {y^ z)^ . 


(4.12) 


The functions fi{y, z) and f 2 {y, z) can be written in terms of the one-loop bubble integral and the one-loop 
box integral in d = 6 — 2e dimensions, Box®: 


hiy,z) = 

1 

yz 

+ y- 

z 


(-3 -I- e -I- 2e^)Bub(si23) + (-h 12 - 8e )Bub(?/si23) 


-h 8 — lOe -|- 3e e j Bub(zsi 23 ) T (—3 -I- 4e -I- e — 2e )Bub(si 23 ) 


-b 


y 


+ 


(l-z)2 

y 

'(l-z) 

1 

(T^' 


+ ( - - + 8 - lOe -b 4e ) Bub(ysi23) 


-12 -b 9e - e jBub(zsi23) -b (6 - 2e - 4e )Bub(si23) + 

(1 - e) ^Bub(zsi23) - Bub(si23) 

(3 - 5e -b 2e®)Bub(7Si23) + (-3 -b 4e -b - 2e®)Bub(si23) 


- 12 -b 8e )Bub(j/Si23) 


(4 - 3e - 3e^ - 2e^) ( Bub(si23) - Bub(7Si23) 
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+ (4 - 9e + 6 e^ - e^)Bub(zsi23) 

+si 23 Box®(j/si 23 , 2Si23, si 23)(1 — 2e) 8 (—1 + e) + —(—2 + 4e — 2e^) + —(6 — 8 e + 2e^) 

\z z z 

1 


+z(—2 + 2e — 86 ) + (4 — 3e + 3e ) H-2(1 — c) ) i 


f2{y,z) = 

1 

yz 

y 


yz 


(3 - 6 - 2e^)Bub(si23) + ( - - 6 + 4e )Bub((l - y- z)si 23 ) 


(4.13) 


+ 7 ( - 6 (1 - e)Bub(zsi23) + (3 - 4e - 6 + 2e'^)Bub(si23) 

+ - 8 + lOe - 4e^^ Bub((l - y - 2 ) 5123 )^ 

+ ^ ( e(l - e)Bub(zsi23) + (-6 + 2e + 4e^)Bub(si23) + ^ - ^ + 12 - 8e^Bub((l - y - z)si23) 

|■2^Bub((l - y- 2 ) 3123 ) - Bub(si23)^ + ^ 2e^Bub(si23) - 2Bub((l - y- 2 ) 3123 ) 


{y + z)' 


( 1 - 2)2 


(1 - e) Bub( 3 i 23 ) - Bub( 23 i 23 ) 


+ ^ ^(3 - 4e - 6^ + 2e^)Bub(3i23) + (-3 + 5e - 2e^)Bub(23i23)^ 

+ ■^^3^(2 + e - 5e^ - 2e^) ^Bub(2Si23) - Bub(3i23)^ 

+ (2 - 7e + 2e^ + 3e^)Bub(23i23) + (-4 + lOe - 4e^)Bub((l - y- 2 ) 3123 ) 

+3i23Box®((1 -y- 2 ) 3123 ,23123, 3i23)(l “ 2e) ^(8 - 4e) - ■^^—^4(1 - e) 

+ (?/ + 2)(-4 + 4e - 66^ - 26^) + ^(-2 + 4e - 26^)^ . (4.14) 

Explicit formulae for the bubble and box integrals are given in the appendix. The square of the Born amplitude 
is given in Eq. (^). 


4.2 The finite part 


The finite remainders of the one- and two-loop contributions to can be decomposed according to their 
colour structure and to their dependence on the number of quark flavours Np. In the two-loop contribution, 
one finds moreover a term proportional to the charge-weighted sum of the quark flavours Np^j] this equals, 
in the case of purely electromagnetic interactions: 


Np^j = 


V e 2 

Z-^q 


(4.15) 


This term originates from diagrams containing a closed quark loop coupling to the virtual photon and which 
first appear at the two-loop level. 

The tree-level combination of invariants 


T = 


2 

yz 


2 2 

y z 


(4.16) 
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frequently occurs in the finite part. We therefore extracted this combination by expressing l/(?/z) by T 
according to the above equation. 


4.2.1 Two-loop contribution to 

The finite remainder of the interference of the two-loop amplitude with the tree-level amplitude is decomposed 
as 




(^ 20 ( 2 /; z) -f ^ 20 ( 2 , y)) + {B 2 o{y, z) + B 2 o{z, y)) 


+ ^ (<^ 20 ( 2 /, z) + C2o{z, y)) + NNp iD2o{y, z) + £> 20 (^, 2 /)) 
+ ^ (£20(2/, z) + E2o{z, y)) + Np {F2o{y, z) + F2oiz, y)) 
+Nf,j f ~ (^20(2/) z) + G2o{z, y)) 


(4.17) 


with 


A2oiy,z) = 


12y 


[2tt^ + 67 r^H( 0 ; z) - 127r2G(l; y) - 12C,3 + 8H(0; z) - 36H(0; 2 )G( 1 , 0; y) - 36H(0,1, 0; z) 


-b39H(l, 0; z) + 39G(1,0; y) + 72G(1,1,0; y)] 

-247r2H(0; z) + 487 r 2 G(l; y) + 288 C 3 + 457 - 84H(0; z) - 36H(0; z)G(0; y) + 144H(0; z)G(l, 0; y) 




-bl44H(0,1,0; z) - 306H(1, 0; z) - 192G(0; y) - 234G(1,0; y) - 288G(1,1, 0; y)] 


36(1-y)= 




-b 67 r 2 H( 0 ; z) + z) - 67 r 2 G(l -z;y) + 187 r 2 G( 0 ; y) - 127r2G(l; y) + 36 C 3 - 36H(0; z)G(l - z,0; y) 

-b60H(0; z)G(0; y) + 72H(0; z)G(0, 0; y) + 36H(0,1,0; z) - 36H(1, 0; z)G(l -z;y) + 36H(1, 0; z)G(0; y) 
-b36H(l, 1,0; z) + 36G(1 - z, 1,0; y) - 355G(0; y) + 270G(0,0; y) - 108G(0,1,0; y) + 6G(1,0; y) 

-72G(1,0, 0; y) + 72G(1,1,0; y)] + + 187 r 2 H( 0 ; z) + 187 r"H(l; z) - 187 r"G(l - z; y) 

-b 547 r 2 G( 0 ; y) - 367 r 2 G(l; y) -f IO 8 C 3 - 277 -f 60H(0; z) - 108H(0; z)G(l - z, 0; y) -b 216H(0; z)G(0; y) 
-b216H(0; z)G(0,0; y) -b 108H(0,1,0; z) -b 36H(1,0; z) - 108H(1, 0; z)G(l - z; y) -b 108H(1, 0; z)G(0; y) 
-bl08H(l, 1, 0; z) -b 108G(1 - z, 1, 0; y) - 615G(0; y) -b 594G(0, 0; y) - 324G(0,1,0; y) -b 198G(1,0; y) 

z IItt^ 

-216G(1, 0,0; y) + 216G(1,1,0; y)] + 2/) - 33H(0; z)G(l - z, 0; y) 

-33H(0,1,0; z) - 33H(1, 0; z) - 33H(1,0; z)G(l - z; y) -b 33H(1,0; z)G(0; y) -b 33H(1,1, 0; z) 

z 227r^ 

+33G(1 - z, 1,0; y) + 33G(0,1, 0; y) - 33G(1,0; y)] + " ^^(1; z) 


-b^—G(1 - z; y) -b 33H(0; z) -b 44H(0; z)G(l - z, 0; y) - 66H(0; z)G(0; y) -b 44H(0,1, 0; z) 

o 

-22H(1,0; z) -b 44H(1, 0; z)G(l - z; y) - 44H(1, 0; z)G(0; y) - 44H(1,1,0; z) - 44G(1 - z, 1, 0; y) 

Z 11TT^ 

-33G(0; y) - 44G(0,1, 0; y) + 110G(1, 0; y)] + ^ - 11H(0; z) + 11H(0; z)G(0; y) 

2 (y-bz)'- 6 
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+11H(1,0; 0 ) + 11G(0; y) - 11G(1,0; y)] + [ - ^H(l; 2 ) + ^G(l - z; y) 

+33H(0; z)G(l - z, 0; y) + 33H(0,1,0; z) + 33H(1, 0; z)G(l -z;y)- 33H(1, 0; z)G(0; y) - 33H(1,1, 0; z) 

z^ Htt^ IItt^ 11 

-33G(1 - z, 1,0; y) - 33G(0,1, 0; y)] + ^^3 [—^ + —^H(l; z)-2^) 

-22H(0; z)G(l - z, 0; y) + 33H(0; z)G(0; y) - 22H(0,1, 0; z) + 33H(1, 0; z) - 22H(1,0; z)G(l - z; y) 
+22H(1,0; z)G(0; y) + 22H(1,1,0; z) + 22G(1 - z, 1, 0; y) + 22G(0,1, 0; y) - 33G(1, 0; y)] 

z 2 1 

+ 2 (y + z)^ [ - - 11H(0; z)G(0; y) - 11H(1, 0; z) + 11G(1, 0; y)] + [ + 2371^ - 127r2H(0; z) 

-127r2H(l; z) + 127r2G(l - z;y) - 367 r 2 G( 0 ; y) + 247r2G(l; y) - 72(3 + 72H(0; z)G(l - z, 0; y) 

-120H(0; z)G(0; y) - 144H(0; z)G(0, 0; y) - 72H(0,1, 0; z) - 18H(1, 0; z) + 72H(1, 0; z)G(l - z; y) 
-72H(1, 0; z)G(0; y) - 72H(1,1,0; z) - 72G(1 - z, 1, 0; y) + 515G(0; y) - 432G(0, 0; y) + 216G(0,1, 0; y) 

-138G(1,0; y) + 144G(1,0, 0; y) - 144G(1,1, 0; y)] + " ^^(1; z) + ^G(l - z; y) 

+14H(0; z)G(l - z, 0; y) + 14H(0,1,0; z) + 14H(1, 0; z)G(l -z;y)- 14H(1, 0; z)G(0; y) - 14H(1,1, 0; z) 

-14G(1 - z, 1,0; y) - 14G(0,1, 0; y)] + + ^H(l; z) - ^G(l -z;y)- 22 

+11H(0; z) - 28H(0; z)G(l - z, 0; y) + 56H(0; z)G(0; y) - 28H(0,1, 0; z) + 56H(1, 0; z) 

-28H(1,0; z)G(l -z;y) + 28H(1,0; z)G(0; y) + 28H(1,1, 0; z) + 28G(1 - z, 1, 0; y) + 11G(0; y) 

Ttt^ 

+28G(0,1,0; y) - 56G(1,0; y)] + \ - 1045 + 147H(0; z) + 36H(0; z)G(l -z;y) + 108H(0; z)G(0; y) 

216 

-36H(0; z)G(l; y) + 72H(0,1; z) + 54H(1; z) + 72H(1; z)G(l -z;y)- 72H(1; z)G(-z; y) 

-36H(1; z)G(l; y) + 72H(1,0; z) + 36H(1,1; z) - 186G(1 - z; y) + 36G(1 - z, 0; y) - 72G(1 - z, 1; y) 
+72G(-z, 1 - z; y) - 72G(0, l-z;y) + 147G(0; y) - 72G(0,1; y) + 36G(1,1 - z; y) + 132G(1; y) 

T 90931 1327r‘* 

-108G(1,0; y) + 72G(1,1; y)] + — [ - + 4776C3 - 216C3H(0; z) + 1080C3H(1; z) 

-864C3G(1 - z; y) - 216C3G(0; y) - 216C3G(1; y) + 304H(0; z) - 1116H(0; z)G(l - z, 0; y) 

-216H(0; z)G(l - z, 1,0; y) + 432H(0; z)G(-z, 1 - z, 0; y) - 432H(0; z)G(0,1 - z, 0; y) 

-144H(0; z)G(0; y) + 1512H(0; z)G(0, 0; y) - 216H(0; z)G(0,1, 0; y) + 216H(0; z)G(l, 1 - z, 0; y) 

-36H(0; z)G(l, 0; y) - 432H(0; z)G(l, 0,0; y) + 1920H(0, 0; z) + 1512H(0,0; z)G(0; y) 

+432H(0,0; z)G(0,0; y) + 864H(0,0,1,0; z) + 1008H(0,1,0; z) - 216H(0,1, 0; z)G(l - z; y) 

+432H(0,1, 0; z)G(-z; y) + 216H(0,1,0; z)G(0; y) - 216H(0,1,0; z)G(l; y) + 432H(0,1,1, 0; z) 
-1095H(1,0; z) - 1116H(1, 0; z)G(l - z; y) + 216H(1,0; z)G(l - z, 0; y) + 432H(1,0; z)G(-z, 1 - z; y) 
-432H(1,0; z)G(-z, 0; y) - 432H(1, 0; z)G(0,1 - z; y) + 1152H(1, 0; z)G(0; y) 

+216H(1,0; z)G(l, 1 - z; y) - 216H(1, 0; z)G(l, 0; y) + 1512H(1,0, 0; z) + 432H(1, 0,0; z)G(0; y) 
+864H(1,0,1,0; z) + 324H(1,1, 0; z) + 432H(1,1,0; z)G(l - z; y) - 432H(1,1,0; z)G(-z; y) 

-216H(1,1, 0; z)G(l; y) + 432H(1,1, 0,0; z) + 216H(1,1,1, 0; z) + 216G(1 - z, 0,1, 0; y) 
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+1116G(1 - z, 1,0; y) + 432G(1 - z, 1,1, 0; y) - 432G(-z, 1 - z, 1,0; y) - 432G(-z, 0,1,0; y) 

+432G(0,1 - z, 1, 0; y) + 304G(0; y) + 1920G(0, 0; y) - 432G(0,0,1, 0; y) - 1008G(0,1, 0; y) 

+432G(0,1,1, 0; y) - 216G(1,1 - z, 1,0; y) + 1095G(1,0; y) - 1512G(1,0, 0; y) + 648G(1,0,1,0; y) 

-792G(1,1,0; y) + 432G(1,1,0, 0; y) - 432G(1,1,1,0; j/)] + ^ [ - ^ + 1 - 14H(0; z)G(0; y) 

-14H(1,0;z) + 14G(1,0;2/)] , (4.18) 

B2oiy,z) = 

2 

^ [ - 3H(0; z)G(l - z]y)- 3H(1; z)G(-z; y) + 3G(-z, 1- z;y)] + ^ [H(0; z)G(1 - z; y) 

+H(1; z)G(-z; y) - G(-z, 1 - z; y)] + ^ [2H(0; z)G(l - z;y) + 2H(1; z)G(-z; y) - 2G(-z, 1 - z; y)] 

2 

ZTT^ 

+ — [3H(0; z) - 24H(0; z)G(l -z;y) + 21H(1; z) - 12H(1; z)G(l -z;y) + 12G(1 - z, 1 - z;y) 

18y 

+G(1 - z; y) + 12G(1 - z, 1; y) + 6G(0,1 - z; y) + 36G(1; y)] + A [ _ 27(3 - 90C3G(1 - z; y) 

9y 

-36H(0; z) - 84H(0; z)G(l - z, 1 - z; y) + 18H(0; z)G(l - z, 1 - z, 0; y) 

+18H(0; z)G(l - z, -z, 1 - z; y) + 152H(0; z)G(l - z; y) - 57H(0; z)G(l - z, 0; y) 

+36H(0; z)G(-z, 1 - z, 1 - z; y) + 78H(0; z)G(-z, 1 - z; y) - 18H(0; z)G(-z, 0,1 - z; y) 

-18H(0; z)G(0,1 - z, 1 - z; y) - 3H(0; z)G(0,1 - z; y) + 36H(0; z)G(0,1 - z, 0; y) 

+18H(0; z)G(0, -z, 1 - z; y) + 54H(0; z)G(l, 0; y) - 18H(0, 0; z)G(l - z, 1 - z; y) 

-36H(0,0; z)G(l - z; y) - 18H(0, 0; z)G(l - z, 0; y) - 18H(0, 0; z)G(0,1 - z; y) 

-9H(0, 0,1; z) + 54H(0,0,1; z)G(l - z; y) - 72H(0,0,1; z)G(-z; y) - 9H(0,1; z) 

+54H(0,1; z)G(l - z, -z; y) + 9H(0,1; z)G(l - z; y) - 18H(0,1; z)G(l - z, 0; y) 

-108H(0,1; z)G(-z, -z; y) + 6H(0,1; z)G(-z; y) + 18H(0,1; z)G(-z, 0; y) + 18H(0,1; z)G(0,1 - z; y) 
-36H(0,1; z)G(0, -z; y) + 9H(0,1,0; z) - 36H(0,1, 0; z)G(l - z; y) + 18H(0,1, 0; z)G(-z; y) 

+72H(1; z)G(l - z, -z, -z; y) - 75H(1; z)G(l - z, -z; y) - 18H(1; z)G(l - z, -z, 0; y) 

— 18H(1; z)G(l — z, 0, —z; y) + 36H(1; z)G(—z, 1 — z, —z; y) — 84H(1; z)G(—z, 1 — z; y) 

+36H(1; z)G(—z, —z, 1 — z; y) — 108H(1; z)G(—z, —z, —z; y) + 84H(1; z)G(—z, —z; y) 

+36H(1; z)G(-z, -z, 0; y) + 143H(1; z)G(-z; y) + 24H(1; z)G(-z, 0; y) - 18H(1; z)G(0, -z, 1 - z; y) 
-18H(1; z)G(0, -z, -z; y) - 3H(1; z)G(0, -z; y) + 36H(1; z)G(0, -z, 0; y) + 9H(1; z)G(0; y) 

+27H(1; z)G(l, 0; y) + 36H(1, 0; z)G(l - z, 1 - z; y) - 18H(1, 0; z)G(l - z, -z; y) 

+57H(1, 0; z)G(l - z; y) - 18H(1, 0; z)G(l - z, 0; y) - 18H(1, 0; z)G(-z, 1 - z; y) 

-78H(1, 0; z)G(-z; y) + 36H(1,0; z)G(0,1 - z; y) - 18H(1, 0; z)G(0, -z; y) - 18H(1,0, 0; z)G(l - z; y) 
-36H(1, 0,1; z) - 18H(1, 0,1; z)G(l - z; y) - 36H(1,1; z)G(-z, -z; y) + 84H(1,1; z)G(-z; y) 

+18H(1,1; z)G(0, -z; y) - 18H(1,1,0; z) - 54H(1,1,0; z)G(l - z; y) + 18H(1,1, 0; z)G(-z; y) 

— 18G(1 — z, 1 — z, 1,0; y) + 75G(1 — z, —z, 1 — z; y) + 18G(1 — z, —z, 1 — z, 0; y) 

—72G(1 — z, —z, —z, 1 — z; y) + 18G(1 — z, —z, 0,1 — z; y) + 18G(1 — z, 0, —z, 1 — z; y) — 9G(1 — z, 0; y) 


14 



+18G(1 - z, 0,1,0; y) + 54G(1 - z, 1,0; y) - 36G(1 - z, 1,1,0; y) + 84G(-z, 1 - z, 1 - z; y) 

—36G(—z, 1 — z, —z, 1 — z; y) — 143G(—z, 1 — z; y) — 24G(—z, 1 — z, 0; y) + 18G(—z, 1 — z, 1, 0; y) 

—36G(—z, —z, 1 — z, 1 — z; y) — 84G(—z, —z, 1 — z; y) — 36G(—z, —z, 1 — z, 0; y) 

+108G(—z, —z, —z, 1 — z; y) — 36G(—z, —z, 0,1 — z; y) — 24G(—z, 0,1 — z; y) + 18G(—z, 0,1, 0; y) 
-9G(0,1 - z; y) + 18G(0, -z, 1 - z, 1 - z; y) + 3G(0, -z, 1 - z; y) - 36G(0, -z, 1 - z, 0; y) 

+18G(0, -z, -z, 1 - z; y) - 36G(0, -z, 0,1 - z; y) + 27G(0,1, 0; y) - 27G(1,1 - z, 0; y) 

-27G(1,0,1 - z; y) + 9G(1,0; y) - 108G(1,1, 0; y)] + j [2H(0; z)G(l - z; y) + 2H(1; z)G(-z; y) 

-2G(-z, 1 - z; y)] + , \ , [2H(0; z)G(l - z; y) + 6H(0; z)G(-z, 1 - z; y) - 6H(0; z)G(0,1 - z; y) 

y(y + z)'- 

-2H(0,1; z) - 6H(0,1; z)G(-z; y) + 6H(1; z)G(-z, 0; y) - 6H(1; z)G(0, -z; y) + 2H(1; z)G(0; y) 

+6H(1, 0; z)G(l - z; y) - 6H(1, 0; z)G(-z; y) - 6H(1,1,0; z) - 2G(1 - z, 0; y) + 6G(1 - z, 1, 0; y) 
-6G(-z, 1 - z, 0; y) - 6G(-z, 0,1 - z; y) - 2G(0,1 - z; y) + 6G(0, -z, 1 - z; y) + 6G(0,1,0; y) 

7^2 

+2G(1,0; y)l + — r - 3 + 3H(0; z) + 24H(0; z)G(l - z; y) - 15H(1; z) + 12H(1; z)G(l - z; y) 

9y 

1 45 

-12G(1 - z, 1 - z; y) - 7G(1 - z; y) - 12G(1 - z, 1; y) - 6G(0,1 - z; y) - 24G(1; y)] + — [yH(l; z) 
45 

-y G(1 - z; y) - 54C3 + 180C3G(1 - z; y) + 139 + 57H(0; z) + 132H(0; z)G(l - z, 1 - z; y) 

-36H(0; z)G(l - z, 1 - z, 0; y) - 36H(0; z)G(l - z, -z, 1 - z; y) - 250H(0; z)G(l - z; y) 

+96H(0; z)G(l - z, 0; y) - 72H(0; z)G(-z, 1 - z, 1 - z; y) - 210H(0; z)G(-z, 1 - z; y) 

+36H(0; z)G(-z, 0,1 - z; y) + 36H(0; z)G(0,1 - z, 1 - z; y) + 96H(0; z)G(0,1 - z; y) 

-72H(0; z)G(0,1 - z, 0; y) - 36H(0; z)G(0, -z, 1 - z; y) - 9H(0; z)G(0; y) - 72H(0; z)G(l, 0; y) 
+36H(0,0; z)G(l - z, 1 - z; y) + 144H(0, 0; z)G(l - z; y) + 36H(0,0; z)G(l - z, 0; y) 

+36H(0, 0; z)G(0,1 - z; y) - 18H(0, 0,1; z) - 108H(0,0,1; z)G(l - z; y) + 144H(0,0,1; z)G(-z; y) 
+36H(0,1; z) - 108H(0,1; z)G(l - z, -z; y) + 36H(0,1; z)G(l - z, 0; y) + 216H(0,1; z)G(-z, -z; y) 
+42H(0,1; z)G(-z; y) - 36H(0,1; z)G(-z, 0; y) - 36H(0,1; z)G(0,1 - z; y) + 72H(0,1; z)G(0, -z; y) 
+18H(0,1, 0; z) + 72H(0,1, 0; z)G(l - z; y) - 36H(0,1,0; z)G(-z; y) - 144H(1; z)G(l - z, -z, -z; y) 
+132H(1; z)G(l - z, -z; y) + 36H(1; z)G(l - z, -z, 0; y) + 36H(1; z)G(l - z, 0, -z; y) 

-72H(1; z)G(-z, 1 - z, -z; y) + 132H(1; z)G(-z, 1 - z; y) - 72H(1; z)G(-z, -z, 1 - z; y) 

+216H(1; z)G(—z, —z, —z; y) — 168H(1; z)G(—z, —z; y) — 72H(1; z)G(—z, —z, 0; y) 

-214H(1; z)G(-z; y) - 30H(1; z)G(-z, 0; y) + 36H(1; z)G(0, -z, 1 - z; y) + 36H(1; z)G(0, -z, -z; y) 
+96H(1; z)G(0, -z; y) - 72H(1; z)G(0, -z, 0; y) - 27H(1; z)G(0; y) - 18H(1; z)G(l, 0; y) - 18H(1,0; z) 
-72H(1,0; z)G(l - z, 1 - z; y) + 36H(1, 0; z)G(l - z, -z; y) - 168H(1,0; z)G(l - z; y) 

+36H(1,0; z)G(l - z, 0; y) + 36H(1, 0; z)G(-z, 1 - z; y) + 210H(1, 0; z)G(-z; y) 

-72H(1,0; z)G(0,1 - z; y) + 36H(1,0; z)G(0, -z; y) + 36H(1,0, 0; z)G(l - z; y) + 36H(1,0,1; z) 
+36H(1,0,1; z)G(l - z; y) + 72H(1,1; z)G(-z, -z; y) - 132H(1,1; z)G(-z; y) - 36H(1,1; z)G(0, -z; y) 
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+72H(1,1, 0; z) + 108H(1,1, 0; z)G(l -z;y)- 36H(1,1, 0; z)G{-z-, y) + 36G(1 -z,l-z,l, 0; y) 

— 132G(1 — z, —z, 1 — z;y) — 36G(1 — z, —z, 1 — 2 , 0;y) + 144G(1 — z, —z, —z, I — z;y) 

-36G(1 - z, -z, 0,1 - z; y) - 36G(1 - z, 0, -z, l-z;y) + 27G(1 - z, 0; y) - 36G(1 - z, 0,1,0; y) 
-108G(1 - z, 1, 0; y) + 72G(1 - z, 1,1,0; y) - 132G(-z, 1 - z, 1 - z; y) + 72G(-z, 1 - z, -z, 1 - z; y) 
+214G(-z, 1 - z; y) + 30G(-z, 1 - z, 0; y) - 36G(-z, 1 - z, 1, 0; y) + 72G(-z, -z, 1 - z, 1 - z; y) 
+168G(—z, —z, 1 — z; y) + 72G(—z, —z, 1 — z, 0; y) — 216G(—z, —z, —z, 1 — z; y) 

+72G(-z, -z, 0,1 - z; y) + 30G(-z, 0,1 - z; y) - 36G(-z, 0,1, 0; y) + 27G(0,1 - z; y) 

—36G(0, —z, 1 — z, 1 — z; y) — 96G(0, —z, 1 — z; y) + 72G(0, —z, 1 — z, 0; y) — 36G(0, —z, —z, 1 — z; y) 
+72G(0, -z, 0,1 - z; y) - 48G(0; y) - 54G(0,1, 0; y) + 18G(1,1 - z, 0; y) + 18G(1,0,1 - z; y) 

z 277r^ 277r^ 277r^ 

-18G(l,0;y) + 144G(1,1,0;y)] + _ ^^2 [" -^H(0;z)-^H(l;z) + G(1 - z;y) 

01 01 01 

- 67 r 2 G( 0 ; y) + 67 r 2 G(l; y) + 225(3 + yH(l; z)G(0; y) - y G(1 - z, 0; y) - y G(0,1 - z; y) 

+9H(0; z)G(l - z, 0; y) - 72H(0; z)G(0,1 - z; y) + 12H(0; z)G(0; y) - 9H(0; z)G(0,0; y) + 81H(0,0,1; z) 
+36H(0,1; z) - 81H(0,1; z)G(l - z; y) - 9H(0,1; z)G(0; y) - 81H(0,1, 0; z) - 81H(1; z)G(l - z, -z; y) 
+36H(1; z)G(-z; y) - 81H(1; z)G(0, -z; y) + 27H(1; z)G(0, 0; y) - 36H(1, 0; z) + 81H(1,0; z)G(l - z; y) 
-9H(1,0; z)G(0; y) + 81H(1,0,1; z) - 81H(1,1,0; z) + 81G(1 - z, -z, 1 - z; y) - 27G(1 - z, 0,0; y) 
-36G(-z, 1 - z; y) - 27G(0,1 - z, 0; y) + 81G(0, -z, 1 - z; y) + 37G(0; y) - 27G(0, 0,1 - z; y) 

-45G(0, 0; y) + 45G(0,1,0; y) + 30G(1,0; y) + 36G(1,0,0; y) - 36G(1,1,0; y)] 

H--- r — 97 r^H( 0 ; z) — 97 r^H(l; z) + 97 r^G(l — z;y) — 367 r^G( 0 ;y) + 367 r^G(l; y) + 54(3 + 62 

18(1-y)^ 

+60H(0; z) - 144H(0; z)G(l - z; y) + 54H(0; z)G(l - z, 0; y) - 108H(0; z)G(0; y) - 54H(0; z)G(0,0; y) 
+54H(0, 0,1; z) - 18H(0,1; z) - 54H(0,1; z)G(l - z; y) - 54H(0,1; z)G(0; y) - 54H(0,1, 0; z) 

+117H(1; z) - 54H(1; z)G(l - z, -z; y) - 162H(1; z)G(-z; y) - 54H(1; z)G(0, -z; y) + 99H(1; z)G(0; y) 
+162H(1; z)G(0,0; y) - 18H(1,0; z) + 54H(1, 0; z)G(l - z; y) - 54H(1, 0; z)G(0; y) + 54H(1, 0,1; z) 
-54H(1,1,0; z) + 54G(1 - z, -z, 1 - z; y) - 117G(1 - z; y) - 99G(1 - z, 0; y) - 162G(1 - z, 0,0; y) 
+162G(-z, 1 - z; y) - 99G(0,1 - z; y) - 162G(0,1 - z, 0; y) + 54G(0, -z, 1 - z; y) + 132G(0; y) 
-162G(0,0,1 - z; y) - 198G(0,0; y) + 270G(0,1, 0; y) + 18G(1, 0; y) + 216G(1,0, 0; y) - 216G(1,1, 0; y)] 

+ 7 —r - 27r2H(l; z) + 27r2G(l - z; y) + 12H(0; z)G(l - z, 0; y) + 12H(0,1, 0; z) + 12H(1, 0; z) 

[y + z^ 

+12H(1, 0; z)G(l - z; y) - 12H(1, 0; z)G(0; y) - 12H(1,1,0; z) - 12G(1 - z, 1, 0; y) - 12G(0,1,0; y) 
z 4:7r^ 47 r^ 

+12G(1,0;y)] + + —H(l; z)-^G(l - z;y)- 6H(0; z) - 8H(0; z)G(l - z, 0; y) 

+12H(0; z)G(0; y) - 8H(0,1,0; z) + 4H(1,0; z) - 8H(1, 0; z)G(l - z; y) + 8H(1, 0; z)G(0; y) 

+8H(1,1, 0; z) + 8G(1 - z, 1,0; y) + 6G(0; y) + 8G(0,1,0; y) - 20G(1, 0; y)] + ^ ^ 

y + z 3 
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r2 O-TT^ 

+2H(0; z) - 2H(0; z)G(0; y) - 2H(1, 0; 2 ) - 2G(0; y) + 2G(1,0; y)] + 2 ) 

27r^ 27r^ 

+—H(l; z) - —G(1 -z;y)- UCs + 4H(0; z)G(0, l-z;y)- 4H(0, 0,1; z) + 4H(0,1; 2)G(1 - z; y) 
+4H(0,1, 0; z) + 4H(1; z)G(l - z, -z; y) + 4H(1; z)G(0, -z; y) - 4H(1,0; z)G(l -z;y)- 4H(1,0,1; z) 

+4H(1,1, 0; z) - 4G(1 - z, -z, 1- z;y)- 4G(0, -z, 1 - z; y)] + ^ [ + 4H(0; z)G(l - z; y) 

+4H(1; z)G(-z; y) - 4G(-z, 1 - z; y)] + ^ [ + 2H(0; z)G(l - z; y) + 2H(1; z)G(-z; y) 

1 -2/ 

-2G(-z, 1 - z; y)] + [27r2H(l; z) - 2n'^G{l - z-y) - 12H(0; z)G(l - z, 0; y) - 12H(0,1,0; z) 

( 2 / + zp 

-12H(1, 0; z)G(l - z; y) + 12H(1, 0; z)G(0; y) + 12H(1,1,0; z) + 12G(1 - z, 1, 0; y) + 12G(0,1,0; y)] 

47r^ 47r^ 

+ . . X 3 [ - 27r2-^H(l; z) + —G(1 - z; y) + 8H(0; z)G(l - z, 0; y) - 12H(0; z)G(0; y) 

(y + z)'^'- 3 3 

+8H(0,1, 0; z) - 12H(1, 0; z) + 8H(1,0; z)G(l - z; y) - 8H(1, 0; z)G(0; y) - 8H(1,1, 0; z) 

2 2 

-8G(1 - z, 1,0; y) - 8G(0,1,0; y) + 12G(1, 0; y)] + [y + 2H(0; z)G(0; y) + 2H(1,0; z) 

-2G(1, 0; y)] + [y + H(0; z)G(0; y) + H(l, 0; z) - G(l, 0; y)] + + 127r2H(0; z) 

+97r2H(l; z) - 97r2G(l - z; y) + 127r2G(0; y) - 2l7r2G(l; y) - 234C3 + 54H(0; z)G(0,1 - z; y) 

+21H(0; z)G(0; y) + 18H(0; z)G(0, 0; y) - 18H(0; z)G(l, 0; y) - 90H(0, 0,1; z) - 63H(0,1; z) 

+90H(0,1; z)G(l - z; y) + 36H(0,1; z)G(0; y) + 54H(0,1,0; z) + 90H(1; z)G(l - z, -z; y) 

-63H(1; z)G(-z; y) + 90H(1; z)G(0, -z; y) - 18H(1; z)G(0; y) - 72H(1; z)G(0,0; y) + 63H(1, 0; z) 
-54H(1, 0; z)G(l - z; y) - 90H(1, 0,1; z) + 54H(1,1,0; z) - 90G(1 - z, -z, 1 - z; y) + 18G(1 - z, 0; y) 
+72G(1 - z, 0,0; y) - 36G(1 - z, 1, 0; y) + 63G(-z, 1 - z; y) + 18G(0,1 - z; y) + 72G(0,1 - z, 0; y) 
-90G(0, -z, 1 - z; y) + 55G(0; y) + 72G(0,0,1 - z; y) + 108G(0, 0; y) - 108G(0,1,0; y) + 3G(1,0; y) 

-90G(1,0, 0; y) + 126G(1,1,0; y)] + [ + 18H(0; z)G(l - z, 1 - z; y) + 27H(0; z)G(l - z; y) 

-18H(0; z)G(-z, 1 - z; y) - 72H(0,0,1; z) + 39H(0,1; z) + 36H(0,1; z)G(l - z; y) 

-90H(0,1; z)G(-z; y) + 18H(0,1; z)G(0; y) + 18H(0,1, 0; z) - 18H(0,1,1; z) + 230H(1; z) 

+54H(1; z)G(l - z, -z; y) - 102H(1; z)G(l - z; y) - 18H(1; z)G(l - z, 0; y) + 36H(1; z)G(-z, 1 - z; y) 
-108H(1; z)G(-z, -z; y) + 66H(1; z)G(-z; y) + 18H(1; z)G(-z, 0; y) - 18H(1; z)G(0,1 - z; y) 

+18H(1; z)G(0, -z; y) - 27H(1; z)G(0; y) - 27H(1,0; z) - 18H(1,0; z)G(l - z; y) + 18H(1,0; z)G(-z; y) 
-36H(1, 0,1; z) + 102H(1,1; z) - 36H(1,1; z)G(-z; y) + 18H(1,1; z)G(0; y) + 18H(1,1,0; z) 

+102G(1 - z, 1 - z; y) + 18G(1 - z, 1 - z, 0; y) - 54G(1 - z, -z, 1 - z; y) - 230G(1 - z; y) 

+18G(1 — z, 0,1 — z; y) + 27G(1 — z, 0; y) — 36G(—z, 1 — z, 1 — z; y) — 66G(—z, 1 — z; y) 

— 18G(—z, 1 — z, 0; y) + 108G(—z, —z, 1 — z; y) — 18G(—z, 0,1 — z; y) + 18G(0,1 — z, 1 — z; y) 
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1 

+27G(0, 1- z;y)- 18G(0, -z, 1 - z; y)] + [-^^(1; z) + —G(1 - z;y) - 170 

-18H(0; z)G(l - z; y) + 9H(0; z)G(l - z, 0; y) - 72H(0,1; z) + 9H(0,1,0; z) - 123H(1; z) 

-90H(1; z)G(-z; y) + 18H(1; z)G(0; y) + 18H(1,0; z) + 9H(1,0; z)G(l -z;y)- 9H(1,0; z)G(0; y) 

-9H(1,1, 0; z) + 123G(1 - z-,y) - 18G(1 - z, 0; y) - 9G(1 - z, 1,0; y) + 90G(-z, 1 - z; y) 

Ttt^ 

-18G(0,1 - z; y) - 9G(0,1,0; y)] + — [ - 115 - 24H(0; z)G(l - z; y) - 12H(0; z)G(l; y) - 12H(0,1; z) 
+7H(1; z) - 48H(1; z)G(l - z; y) + 36H(1; z)G(0; y) - 12H(1; z)G(l; y) - 12H(1,0; z) + 12H(1,1; z) 
+48G(1 - z, 1 - z; y) + 19G(1 - z; y) - 24G(1 - z, 0; y) - 36G(0,1 - z; y) + 48G(0,1; y) 

T 37r^ 1 5251 357 

+12G(1,1 — z; y) — 26G(1; y) + 36G(1,0; y) — 48G(1,1; y)] + ^ [-^—I- ^2 - 

+108C3H(1; z) - 270C3G(1 - z; y) + 162C3G(1; y) - 360H(0; z) - 198H(0; z)G(l - z, 1 - z; y) 

+108H(0; z)G(l - z, 1 - z, 0; y) + 108H(0; z)G(l - z, -z, 1 - z; y) + 78H(0; z)G(l - z; y) 

+54H(0; z)G(l - z, 0,1 - z; y) - 180H(0; z)G(l - z, 0; y) - 108H(0; z)G(l - z, 0, 0; y) 

+54H(0; z)G(l - z, 1,0; y) + 108H(0; z)G(-z, 1 - z, 1 - z; y) + 297H(0; z)G(-z, 1 - z; y) 

-108H(0; z)G(-z, -z, 1 - z; y) - 108H(0; z)G(0,1 - z, 1 - z; y) - 180H(0; z)G(0,1 - z; y) 

+216H(0; z)G(0, -z, 1 - z; y) - 216H(0; z)G(0; y) + 54H(0; z)G(0,1, 0; y) - 54H(0; z)G(l, 1 - z, 0; y) 
-108H(0; z)G(l, 0,1 - z; y) + 9H(0; z)G(l, 0; y) + 108H(0; z)G(l, 0, 0; y) - 378H(0, 0; z)G(l - z; y) 
-108H(0, 0; z)G(l - z, 0; y) - 108H(0,0; z)G(0,1 - z; y) + 18H(0,0,1; z) + 216H(0, 0,1; z)G(l - z; y) 
-432H(0,0,1; z)G(-z; y) - 108H(0, 0,1; z)G(0; y) + 54H(0, 0,1; z)G(l; y) + 216H(0,0,1,0; z) 

+348H(0,1; z) + 324H(0,1; z)G(l - z, -z; y) - 279H(0,1; z)G(l - z; y) - 162H(0,1; z)G(l - z, 0; y) 
+108H(0,1; z)G(-z, 1 - z; y) - 540H(0,1; z)G(-z, -z; y) - 63H(0,1; z)G(-z; y) 

+ 108H(0,1; z)G(-z, 0; y) + 54H(0,1; z)G(0,1 - z; y) - 108H(0,1; z)G(0, -z; y) - 117H(0,1; z)G(0; y) 

+ 108H(0,1; z)G(0,0; y) - 54H(0,1; z)G(l, 1 - z; y) + 54H(0,1; z)G(l, 0; y) - 252H(0,1,0; z) 

-216H(0,1, 0; z)G(l - z; y) + 108H(0,1,0; z)G(-z; y) - 54H(0,1,0; z)G(0; y) - 54H(0,1, 0; z)G(l; y) 
+198H(0,1,1; z) - 108H(0,1,1; z)G(-z; y) + 108H(0,1,1, 0; z) + 17H(1; z) 

+432H(1; z)G(l - z, -z, -z; y) - 477H(1; z)G(l - z, -z; y) - 108H(1; z)G(l - z, -z, 0; y) 

+297H(1; z)G(l - z; y) - 108H(1; z)G(l - z, 0, -z; y) + 198H(1; z)G(l - z, 0; y) 

+162H(1; z)G(l - z, 1,0; y) + 216H(1; z)G(-z, 1 - z, -z; y) - 396H(1; z)G(-z, 1 - z; y) 

— 108H(1; z)G(—z, 1 — z, 0; y) + 216H(1; z)G(—z, —z, 1 — z; y) — 648H(1; z)G(—z, —z, —z; y) 

+234H(1; z)G(-z, -z; y) + 108H(1; z)G(-z, -z, 0; y) + 426H(1; z)G(-z; y) 

-108H(1; z)G(-z, 0,1 - z; y) + 108H(1; z)G(-z, 0, -z; y) - 297H(1; z)G(-z, 0; y) 

-54H(1; z)G(0,1 - z, -z; y) + 198H(1; z)G(0,1 - z; y) - 108H(1; z)G(0, -z, 1 - z; y) 

+ 108H(1; z)G(0, -z, -z; y) - 297H(1; z)G(0, -z; y) + 216H(1; z)G(0, -z, 0; y) - 78H(1; z)G(0; y) 
+108H(1; z)G(0, 0, -z; y) + 378H(1; z)G(0,0; y) - 54H(1; z)G(0,1, 0; y) - 54H(1; z)G(l, 1 - z, -z; y) 
-54H(1; z)G(l, 0, -z; y) - 81H(1; z)G(l, 0; y) - 108H(1; z)G(l, 0,0; y) - 81H(1,0; z) 
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+216H(1,0; z)G(l -z,l-z]y)- 108H(1,0; z)G(l - 2 , -z; y) + 117H(1, 0; z)G(l - z; y) 

-162H(1,0; z)G(l - z, 0; y) - 108H(1, 0; z)G(-z, 1 - z; j/) + 108H(1, 0; z)G(-z, -z; y) 

-297H(1,0; z)G(-z; y) - 216H(1, 0; z)G(0, -z; y) + 171H(1,0; z)G(0; y) 

+108H(1,0; z)G(0,0; y) + 54H(1,0; z)G(l, 1 - z; j/) + 54H(1,0; z)G(l, 0; y) - 108H(1, 0,0; z)G(l - z; y) 
-54H(1,0, 0,1; z) + 360H(1, 0,1; z) - 162H(1, 0,1; z)G(l - z]y) - 108H(1, 0,1; z)G(-z; y) 

+54H(1,0,1; z)G(0; y) + 54H(1,0,1; z)G(l; y) + 108H(1, 0,1, 0; z) - 297H(1,1; z) 

-216H(1,1; z)G(-z, -z; y) + 396H(1,1; z)G(-z; y) + 108H(1,1; z)G(-z, 0; y) + 108H(1,1; z)G(0, -z; y) 
-198H(1,1; z)G(0; y) + 81H(1,1,0; z) - 378H(1,1, 0; z)G(l - z; y) + 108H(1,1, 0; z)G(-z; y) 

+108H(1,1, 0; z)G(0; y) - 54H(1,1, 0; z)G(l; y) + 108H(1,1, 0,1; z) + 162H(1,1,1, 0; z) 

-297G(1 - z, 1 - z; y) - 198G(1 - z, 1 - z, 0; y) - 216G(1 - z, 1 - z, 1, 0; y) + 477G(1 - z, -z, 1 - z; y) 
+108G(1 — z, —z, 1 — z, 0; y) — 432G(1 — z, —z, —z, 1 — z; y) + 108G(1 — z, —z, 0,1 — z; y) 

-17G(1 -z;y)- 198G(1 - z, 0,1 - z; y) + 108G(1 - z, 0, -z, l-z;y) + 78G(1 - z, 0; y) 

-378G(1 - z, 0, 0; y) + 108G(1 - z, 0,1,0; y) - 162G(1 - z, 1,1 - z, 0; y) - 162G(1 - z, 1,0,1 - z; y) 
+81G(1 - z, 1, 0; y) + 108G(1 - z, 1, 0,0; y) + 396G(-z, 1 - z, 1 - z; y) + 108G(-z, 1 - z, 1 - z, 0; y) 
-216G(-z, 1 - z, -z, 1 - z; y) - 426G(-z, 1 - z; y) + 108G(-z, 1 - z, 0,1 - z; y) + 297G(-z, 1 - z, 0; y) 

—216G(—z, —z, 1 — z, 1 — z; y) — 234G(—z, —z, 1 — z; y) — 108G(—z, —z, 1 — z, 0; y) 

+648G(—z, —z, —z, 1 — z; y) — 108G(—z, —z, 0,1 — z; y) + 108G(—z, 0,1 — z, 1 — z; y) 

+297G(-z, 0,1 - z; y) - 108G(-z, 0, -z, 1 - z; y) - 198G(0,1 - z, 1 - z; y) + 54G(0,1 - z, -z, 1 - z; y) 

+78G(0,1 - z; y) - 378G(0,1 - z, 0; y) + 162G(0,1 - z, 1,0; y) + 108G(0, -z, 1 - z, 1 - z; y) 

+297G(0, -z, 1 - z; y) - 216G(0, -z, 1 - z, 0; y) - 108G(0, -z, -z, 1 - z; y) - 216G(0, -z, 0,1 - z; y) 
-360G(0; y) - 378G(0,0,1 - z; y) - 108G(0,0, -z, 1 - z; y) + 108G(0,0,1,0; y) + 54G(0,1,1 - z, 0; y) 
+54G(0,1,0,1 - z; y) + 333G(0,1, 0; y) - 216G(0,1,1, 0; y) + 54G(1,1 - z, -z, 1 - z; y) 

+81G(1,1 - z, 0; y) + 108G(1,1 - z, 0, 0; y) + 81G(1, 0,1 - z; y) + 108G(1,0,1 - z, 0; y) 

+54G(1,0, -z, 1 - z; y) + 3G(1, 0; y) + 108G(1,0, 0,1 - z; y) + 378G(1, 0,0; y) - 216G(1,0,1,0; y) 

+117G(1,1,0; y) - 216G(1,1,0, 0; y) + 216G(1,1,1,0; y)] + ^ [H + 9H(0; z) - 24H(0; z)G(l - z; y) 

-6H(0,1; z) + 8H(1; z) - 36H(1; z)G(l - z; y) + 24H(1; z)G(0; y) + 6H(1; z)G(l; y) + 24H(1,1; z) 

+24G(1 - z, 1 - z; y) + 2G(1 - z; y) - 24G(1 - z, 0; y) + 12G(1 - z, 1; y) - 6G(0,1 - z; y) + 9G(0; y) 

+12G(0,1; y) - 6G(1,1 - z; y) - 10G(1; y) + 24G(1, 0; y) - 6G(1,1; y)] + ^ [ 288 C 3 + 180C3H(1; z) 

-360C3G(1 - z; y) + 180C3G(1; y) - 188H(0; z) - 150H(0; z)G(l - z, 1 - z; y) 

+72H(0; z)G(l - z, 1 - z, 0; y) + 72H(0; z)G(l - z, -z, 1 - z; y) + 295H(0; z)G(l - z; y) 

-300H(0; z)G( 1 - z, 0; y) - 36H(0; z)G(l - z, 0, 0; y) + 72H(0; z)G(-z, 1 - z, 1 - z; y) 

+216H(0; z)G(-z, 1 - z; y) - 72H(0; z)G(-z, -z, 1 - z; y) - 36H(0; z)G(0,1 - z, 1 - z; y) 

-156H(0; z)G(0,1 - z; y) + 36H(0; z)G(0,1 - z, 0; y) + 144H(0; z)G(0, -z, 1 - z; y) + 78H(0; z)G(0; y) 
-72H(0; z)G(0,0,1 - z; y) + 18H(0; z)G(0,0; y) + 36H(0; z)G(0,1, 0; y) - 72H(0; z)G(l, 1 - z, 0; y) 
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-36H(0; z)G(l, 0,1 - z; y) + 132H(0; z)G(l, 0; y) + 36H(0; z)G(l, 0,0; y) + 36H(0; z)G(l, 1,0; y) 
+36H(0,0; z) - 36H(0, 0; z)G(l - z, 1 - z; y) - 108H(0, 0; z)G(l - z;y) - 36H(0, 0; z)G(l - z, 0; y) 
-36H(0, 0; z)G(0,1 - z; y) + 18H(0, 0; z)G(0; y) - 108H(0, 0,0,1; z) + 84H(0,0,1; z) 

+180H(0,0,1; z)G(l - z; y) - 288H(0, 0,1; z)G(-z; y) + 36H(0, 0,1; z)G(l; y) + 108H(0,0,1,0; z) 
-36H(0, 0,1,1; z) + 289H(0,1; z) + 216H(0,1; z)G(l - z, -z; y) - 222H(0,1; z)G(l - z; y) 

-72H(0,1; z)G(l - z, 0; y) + 72H(0,1; z)G(-z, 1 - z; y) - 360H(0,1; z)G(-z, -z; y) 

+192H(0,1; z)G(-z; y) + 72H(0,1; z)G(-z, 0; y) - 72H(0,1; z)G(0, -z; y) - 60H(0,1; z)G(0; y) 
+36H(0,1; z)G(0,0; y) - 36H(0,1; z)G(l, 1 - z; y) - 114H(0,1,0; z) - 144H(0,1,0; z)G(l - z; y) 
+72H(0,1, 0; z)G(-z; y) - 72H(0,1, 0; z)G(0; y) + 36H(0,1, 0; z)G(l; y) - 72H(0,1,0,1; z) 

+150H(0,1,1; z) - 72H(0,1,1; z)G(-z; y) + 36H(0,1,1; z)G(0; y) + 36H(0,1,1, 0; z) - 376H(1; z) 
+288H(1; z)G(l - z, -z, -z; y) - 372H(1; z)G(l - z, -z; y) - 72H(1; z)G(l - z, -z, 0; y) 

+204H(1; z)G(l - z; y) - 72H(1; z)G(l - z, 0, -z; y) + 150H(1; z)G(l - z, 0; y) 

+36H(1; z)G(l - z, 0, 0; y) + 36H(1; z)G(l - z, 1, 0; y) + 144H(1; z)G(-z, 1 - z, -z; y) 

-300H(1; z)G(-z, 1 - z; y) - 72H(1; z)G(-z, 1 - z, 0; y) + 144H(1; z)G(-z, -z, 1 - z; y) 

-432H(1; z)G(-z, -z, -z; y) + 408H(1; z)G(-z, -z; y) + 72H(1; z)G(-z, -z, 0; y) 

+584H(1; z)G(-z; y) - 72H(1; z)G(-z, 0,1 - z; y) + 72H(1; z)G(-z, 0, -z; y) - 216H(1; z)G(-z, 0; y) 
-36H(1; z)G(0,1 - z, -z; y) + 150H(1; z)G(0,1 - z; y) + 36H(1; z)G(0,1 - z, 0; y) 

-72H(1; z)G(0, -z, 1 - z; y) + 72H(1; z)G(0, -z, -z; y) - 216H(1; z)G(0, -z; y) 

+144H(1; z)G(0, -z, 0; y) - 295H(1; z)G(0; y) + 36H(1; z)G(0, 0,1 - z; y) - 36H(1; z)G(0,0, -z; y) 

+ 108H(1; z)G(0, 0; y) - 36H(1; z)G(l, 1 - z, -z; y) - 36H(1; z)G(l, 0, -z; y) - 72H(1; z)G(l, 0, 0; y) 
+15H(1,0; z) + 108H(1,0; z)G(l - z, 1 - z; y) - 72H(1,0; z)G(l - z, -z; y) + 6H(1, 0; z)G(l - z; y) 
-72H(1,0; z)G(l - z, 0; y) - 72H(1,0; z)G(-z, 1 - z; y) + 72H(1,0; z)G(-z, -z; y) 

-216H(1,0; z)G(-z; y) + 108H(1, 0; z)G(0,1 - z; y) - 144H(1, 0; z)G(0, -z; y) + 168H(1,0; z)G(0; y) 
+36H(1,0; z)G(0,0; y) - 36H(1,0; z)G(l, 1 - z; y) + 72H(1, 0; z)G(l, 0; y) + 18H(1, 0,0; z) 

-36H(1,0, 0; z)G(l - z; y) - 72H(1,0, 0,1; z) + 222H(1, 0,1; z) - 36H(1,0,1; z)G(l - z; y) 

-72H(1,0,1; z)G(-z; y) + 36H(1,0,1; z)G(0; y) + 36H(1,0,1; z)G(l; y) + 72H(1, 0,1, 0; z) 

-204H(1,1; z) - 144H(1,1; z)G(-z, -z; y) + 300H(1,1; z)G(-z; y) + 72H(1,1; z)G(-z, 0; y) 

+72H(1,1; z)G(0, -z; y) - 150H(1,1; z)G(0; y) - 36H(1,1; z)G(0, 0; y) + 54H(1,1, 0; z) 

-216H(1,1, 0; z)G(l - z; y) + 72H(1,1,0; z)G(-z; y) + 36H(1,1,0; z)G(0; y) + 36H(1,1, 0; z)G(l; y) 
+108H(1,1,1,0; z) - 204G(1 - z, 1 - z; y) - 150G(1 - z, 1 - z, 0; y) - 36G(1 - z, 1 - z, 0,0; y) 
-108G(1 - z, 1 - z, 1, 0; y) + 372G(1 - z, -z, 1 - z; y) + 72G(1 - z, -z, 1 - z, 0; y) 

-288G(1 - z, -z, -z, 1 - z; y) + 72G(1 - z, -z, 0,1 - z; y) + 376G(1 - z; y) - 150G(1 - z, 0,1 - z; y) 
-36G(1 - z, 0,1 - z, 0; y) + 72G(1 - z, 0, -z, 1 - z; y) + 295G(1 - z, 0; y) - 36G(1 - z, 0,0,1 - z; y) 
-108G(1 - z, 0, 0; y) + 108G(1 - z, 0,1,0; y) - 36G(1 - z, 1,1 - z, 0; y) - 36G(1 - z, 1, 0,1 - z; y) 
+144G(1 - z, 1, 0; y) + 72G(1 - z, 1, 0,0; y) - 72G(1 - z, 1,1,0; y) + 300G(-z, 1 - z, 1 - z; y) 
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+72G{—z, 1 — z,l — z,0;y) — 144G(— 2 , 1 — z, —z, 1 — z;y) — 584G(—z, 1 — z;y) 

+72G(—z, 1 — z, 0,1 — z;y) + 216G(—z, 1 — z, 0; y) — 144G(—z, —z, 1 — z, 1 — z; y) 

—408G(—z, —z, 1 — z; y) — 72G(—z, —z, 1 — z, 0; y) + 432G(—z, —z, —z, 1 — z; y) 

-72G(-z, -z, 0,1 - z; y) + 72G(-z, 0,1 - z, 1 - z; y) + 216G(-z, 0,1 - z; y) 

-72G(-z, 0, -z, 1 - z; y) - 150G(0,1 - z, 1 - z; y) - 36G(0,1 - z, 1 - z, 0; y) 

+36G(0,1 - z, -z, 1 - z; y) + 295G(0,1 - z; y) - 36G(0,1 - z, 0,1 - z; y) - 108G(0,1 - z, 0; y) 

+72G(0,1 - z, 1,0; y) + 72G(0, -z, 1 - z, 1 - z; y) + 216G(0, -z, 1 - z; y) - 144G(0, -z, 1 - z, 0; y) 
-72G(0, -z, -z, 1 - z; y) - 144G(0, -z, 0,1 - z; y) - 188G(0; y) - 36G(0, 0,1 - z, 1 - z; y) 

-108G(0,0,1 - z; y) + 36G(0,0, -z, 1 - z; y) + 36G(0,0; y) + 108G(0,0,1,0; y) + 6G(0,1,0; y) 

-72G(0,1,1, 0; y) + 36G(1,1 - z, -z, 1 - z; y) + 72G(1,1 - z, 0,0; y) + 72G(1,1 - z, 1,0; y) 

+72G(1,0,1 - z, 0; y) + 36G(1,0, -z, 1 - z; y) - 310G(1,0; y) + 72G(1,0,0,1 - 2 ; y) 

+90G(1,0,0; y) - 144G(1, 0,1,0; y) + 60G(1,1,0; y) - 108G(1,1, 0,0; y) + 36G(1,1,1, 0; y)] , (4.19) 

C2o{y,z) = 

2 

+ ^ [6H(0; z)G(l - z; y) + 6H(1; z)G(-z; y) - 6G(-z, 1 - z; y)] + ^ [ “ 2H(0; z)G(l - z; y) 

-2H(1; z)G(-z; y) + 2G(-z, 1 - z; y)] + ^ [ - 4H(0; z)G(l - z; y) - 4H(1; z)G(-z; y) 

V 

Z'k‘^ 

+4G(-z, 1 - z; y)] + — [ - 1 + 4H(0; z) + 10H(0; z)G(l - z; y) + H(l; z) + 6H(1; z)G(l - z; y) 

+10H(1; z)G(-z; y) - 16G(1 - z, 1 - z; y) + G(1 - z; y) + 2G(1 - z, 0; y) + 4G(1 - z, 1; y) 

-10G(-z, 1 - z; y)] + ^ [ - 36 C 3 - 16C3G(1 - z-,y) - 4H(0; z) + 6H(0; z)G(l - z, 1 - z; y) 

-40H(0; z)G(l - z, 1 - z, 0; y) + 9H(0; z)G(l - z;y) - 8H(0; z)G(l - z, 0,1 - z; y) 

-24H(0; z)G(l - z, 0; y) - 16H(0; z)G(-z, 1 - z, 1 - z; y) + 46H(0; z)G(-z, 1 - z; y) 

-40H(0; z)G(-z, 1 - z, 0; y) - 10H(0; z)G(0,1 - z; y) + 12H(0; z)G(l, 0; y) 

+16H(0,0; z)G(l - z, 1 - z; y) + 40H(0, 0; z)G(l - z, 0; y) - 4H(0, 0,1; z) + 17H(0,1; z) 

+24H(0,1; z)G(l - z, 1 - z; y) - 14H(0,1; z)G(l - z; y) - 8H(0,1; z)G(l - z, 0; y) 

-8H(0,1; z)G(-z, 1 - z; y) - 24H(0,1; z)G(-z, -z; y) + 14H(0,1; z)G(-z; y) - 8H(0,1, 0; z) 

+32H(0,1,0; z)G(l - z; y) - 40H(0,1,0; z)G(-z; y) - 6H(0,1,1; z) + 8H(0,1,1; z)G(-z; y) 

+24H(1; z)G(l — z, 1 — z, —z; y) — 8H(1; z)G(l — z, —z; y) — 16H(1; z)G(l — z, 0, —z; y) 

-6H(1; z)G(l - z, 0; y) - 24H(1; z)G(-z, 1 - z, -z; y) + 12H(1; z)G(-z, 1 - z; y) 

+8H(1; z)G(-z, 1 - z, 0; y) - 24H(1; z)G(-z, -z, 1 - z; y) - 24H(1; z)G(-z, -z, -z; y) 

+60H(1; z)G(-z, -z; y) + 24H(1; z)G(-z, -z, 0; y) + 26H(1; z)G(-z; y) + 8H(1; z)G(-z, 0,1 - z; y) 
+18H(1; z)G(-z, 0; y) - 6H(1; z)G(0,1 - z; y) - 10H(1; z)G(0, -z; y) - 17H(1; z)G(0; y) 

+24H(1; z)G(l, 0; y) - 13H(1, 0; z) - 40H(1, 0; z)G(l - z, 1 - z; y) + 4H(1, 0; z)G(l - z; y) 

-24H(1, 0; z)G(-z, 1 - z; y) - 46H(1,0; z)G(-z; y) + 40H(1,0; z)G(-z, 0; y) + 40H(1,0, 0; z)G(l - z; y) 
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-10H(1,0,1; z) - 16H(1,0,1; z)G(l -z-,y) + 8H(1, 0,1; z)G{-z; y) + 24H(1,1; z)G{-z, -z; y) 
-12H(1,1; z)G(-z; y) - 8H(1,1; z)G{-z, 0; y) + 6H(1,1; z)G(0; y) - 20H(1,1,0; z) 

+8H(1,1, 0; z)G(l -z]y) + 24H(1,1, 0; z)G(-z; y) - 24G(1 - z, 1 - z, -z, 1 - z; y) 

+6G(1 — z, 1 — z, 0; y) + 48G(1 — z, 1 — z, 1,0; y) + 8G(1 — z, —z, 1 — z; y) + 6G(1 — z, 0,1 — z; y) 
+16G(1 - z, 0, -z, 1 - z; y) + 17G(1 - z, 0; y) - 48G(1 - z, 0,1, 0; y) + 18G(1 - z, 1,0; y) 

— 16G(1 — z, 1,1,0; t/) — 12G(—z, 1 — z, 1 — z; t/) — 8G(—z, 1 — z, 1 — z, 0; y) 


+24G(—z, 1 — z, —z, 1 — z; y) — 26G(—z, 1 — z; y) — 8G(—z, 1 — z, 0,1 — z; y) — 18G(—z, 1 — z, 0; y) 
+64G(—z, 1 — z, 1, 0; y) + 24G(—z, —z, 1 — z, 1 — z; y) — 60G(—z, —z, 1 — z; y) 

—24G(—z, —z, 1 — z, 0; y) + 24G(—z, —z, —z, 1 — z; y) — 24G(—z, —z, 0,1 — z; y) 


-8G(-z, 0,1 - z, 1 - z; y) - 18G(-z, 0,1 - z; y) + 64G(-z, 0,1, 0; y) + 6G(0,1 - z, 1 - z; y) 
+17G(0,1 - z; y) + 10G(0, -z, 1 - z; y) + 18G(0,1, 0; y) - 24G(1,1 - z, 0; y) - 24G(1,0,1 - z; y) 

-30G(1, 0; y)] + ^ [ - 2H(0; z)G(l - z; y) - 2H(1; z)G(-z; y) + 2G(-z, 1 - z; y)] 


1 


■ ‘ - 6H(0; z)G(l - z, 1 - z; y) - 17H(0; z)G(l - z; y) + 18H(0; z)G(-z, 1 - z; y) 


2y(y + z) 

-18H(0; z)G(0,1 - z; y) + 17H(0,1; z) + 6H(0,1; z)G(l - z; y) - 18H(0,1; z)G(-z; y) - 6H(0,1,1; z) 
-6H(1; z)G(l - z, 0; y) + 18H(1; z)G(-z, 0; y) - 6H(1; z)G(0,1 - z; y) - 18H(1; z)G(0, -z; y) 
-17H(1; z)G(0; y) - 17H(1,0; z) + 24H(1, 0; z)G(l - z; y) - 18H(1, 0; z)G(-z; y) - 6H(1,0,1; z) 
+6H(1,1; z)G(0; y) - 24H(1,1,0; z) + 6G(1 - z, 1 - z, 0; y) + 6G(1 - z, 0,1 - z; y) + 17G(1 - z, 0; y) 
+18G(1 - z, 1,0; y) - 18G(-z, 1 - z, 0; y) - 18G(-z, 0,1 - z; y) + 6G(0,1 - z, 1 - z; y) 

+17G(0,1 - z; y) + 18G(0, -z, 1 - z; y) + 18G(0,1, 0; y) - 34G(1, 0; y)] + ^ [ - 4H(0; z)G(l - z; y) 

+H(1; z) - 6H(1; z)G(l - z; y) - 4H(1; z)G(-z; y) + 10G(1 - z, 1 - z; y) - 6G(1 - z; y) 


—2G(1 — z, 0; y) — 4G(1 — z, 1; y) + 4G(—z, 1 — z; y) + G(l; y)] + ^ [ — 8^3 + 80C3G(1 — z; y) + 19 

4y 

+8H(0; z) - 12H(0; z)G(l - z, 1 - z; y) + 32H(0; z)G(l - z, 1 - z, 0; y) + 26H(0; z)G(l - z; y) 

+16H(0; z)G(l - z, 0,1 - z; y) + 40H(0; z)G(l - z, 0; y) + 32H(0; z)G(-z, 1 - z, 1 - z; y) 

-108H(0; z)G(-z, 1 - z; y) + 32H(0; z)G(-z, 1 - z, 0; y) + 68H(0; z)G(0,1 - z; y) - 16H(0; z)G(l, 0; y) 
-32H(0,0; z)G(l - z, 1 - z; y) - 32H(0,0; z)G(l - z, 0; y) - 8H(0, 0,1; z) - 26H(0,1; z) 

-48H(0,1; z)G(l - z, 1 - z; y) + 36H(0,1; z)G(l - z; y) + 16H(0,1; z)G(l - z, 0; y) 

+16H(0,1; z)G(-z, 1 - z; y) + 48H(0,1; z)G(-z, -z; y) - 28H(0,1; z)G(-z; y) + 16H(0,1,0; z) 

-16H(0,1, 0; z)G(l - z; y) + 32H(0,1, 0; z)G(-z; y) + 12H(0,1,1; z) - 16H(0,1,1; z)G(-z; y) - 2H(1; z) 
-48H(1; z)G(l - z, 1 - z, -z; y) + 24H(1; z)G(l - z, -z; y) + 32H(1; z)G(l - z, 0, -z; y) 

+12H(1; z)G(l - z, 0; y) + 48H(1; z)G(-z, 1 - z, -z; y) - 24H(1; z)G(-z, 1 - z; y) 

— 16H(1; z)G(—z, 1 — z, 0; y) + 48H(1; z)G(—z, —z, 1 — z; y) + 48H(1; z)G(—z, —z, —z; y) 

-136H(1; z)G(-z, -z; y) - 48H(1; z)G(-z, -z, 0; y) - 16H(1; z)G(-z, 0,1 - z; y) - 20H(1; z)G(-z, 0; y) 
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+12H(1; z)G(0, l-z]y) + 68H(1; z)G(0, -z; y) + 30H(1; z)G(0; y) - 24H(1; z)G(l, 0; y) + 26H(1,0; z) 
+32H(1,0; z)G(l - z, 1 - z; y) - 56H(1, 0; z)G(l - z; y) + 108H(1, 0; z)G(-z; y) - 32H(1,0; z)G(-z, 0; y) 
-32H(1,0, 0; z)G(l - z;y) + 4H(1,0,1; z) + 32H(1,0,1; z)G(l - z]y) - 16H(1, 0,1; z)G(-z; y) 

-48H(1,1; z)G(-z, -z; y) + 24H(1,1; z)G(-z; y) + 16H(1,1; z)G(-z, 0; y) - 12H(1,1; z)G(0; y) 

+56H(1,1,0; z) - 16H(1,1, 0; z)G(l - z;y) + 48G(1 - z, 1 - z, -z, 1 - z; y) - 12G(1 - z, 1 - z, 0; y) 
-48G(1 - z, 1 - z, 1,0; y) - 24G(1 - z, -z, l-z;y) + 2G(1 - z;y) - 12G(1 - z, 0,1 - z; y) 

-32G(1 - z, 0, -z, l-z;y)- 30G(1 - z, 0; y) + 48G(1 - z, 0,1,0; y) - 36G(1 - z, 1,0; y) 

+32G(1 - z, 1,1,0; y) + 24G(-z, 1 - z, 1 - z; y) + 16G(-z, 1 - z, 1 - z, 0; y) 

—48G(—z, 1 — z, —z, 1 — z; y) + 16G(—z, 1 — z, 0,1 — z; y) + 20G(—z, 1 — z, 0; y) — 80G(—z, 1 — z, 1,0; y) 
—48G(—z, —z, 1 — z, 1 — z; y) + 136G(—z, —z, 1 — z; y) + 48G(—z, —z, 1 — z, 0; y) 

—48G(—z, —z, —z, 1 — z; y) + 48G(—z, —z, 0,1 — z; y) + 16G(—z, 0,1 — z, 1 — z; y) + 20G(—z, 0,1 — z; y) 
-80G(-z, 0,1,0; y) - 12G(0,1 - z, 1 - z; y) - 30G(0,1 - z; y) - 68G(0, -z, 1 - z; y) - 36G(0,1, 0; y) 
+24G(1,1 - z, 0; y) + 24G(1,0,1 - z; y) + 60G(1, 0; y) - 8G(1,1, 0; y)] + ^(^4^ [ + 

1 4-7r^ 147r^ 147r^ 27r^ 47r^ 

+ ^H(0; z) + ^H(l; z) - ^G(l - z; y) + —G(0; y) - —G(l; y) - 92C3 

+28H(0; z)G(0,1 - z; y) + 12H(0; z)G(0; y) + 4H(0; z)G(0, 0; y) - 32H(0, 0,1; z) - 16H(0,1; z) 

+32H(0,1; z)G(l - z; y) + 4H(0,1; z)G(0; y) + 28H(0,1,0; z) + 32H(1; z)G(l - z, -z; y) 

-16H(1; z)G(-z; y) + 32H(1; z)G(0, -z; y) + 6H(1; z)G(0; y) - 4H(1; z)G(0, 0; y) + 16H(1,0; z) 

-28H(1,0; z)G(l - z; y) - 32H(1, 0,1; z) + 28H(1,1, 0; z) - 32G(1 - z, -z, 1 - z; y) - 6G(1 - z, 0; y) 
+4G(1 - z, 0, 0; y) - 4G(1 - z, 1,0; y) + 16G(-z, 1 - z; y) - 6G(0,1 - z; y) + 4G(0,1 - z, 0; y) 

-32G(0, -z, 1 - z; y) + 23G(0; y) + 4G(0,0,1 - z; y) - 10G(0, 0; y) - 8G(0,1, 0; y) - 14G(1,0; y) 

z 1 Itt^ 

-8G(1, 0,0; y) + 8G(1,1, 0; y)] + [—^ + 67r2H(0; z) + 67r2H(l; z) - - z; y) 

+27r2G(0; y) - y) - 132C3 + 17 + 4H(0; z) + 28H(0; z)G(l - z; y) + 36H(0; z)G(0,1 - z; y) 

+12H(0; z)G(0,0; y) - 48H(0, 0,1; z) + 4H(0,1; z) + 48H(0,1; z)G(l - z; y) + 12H(0,1; z)G(0; y) 

+36H(0,1,0; z) - 2H(1; z) + 48H(1; z)G(l - z, -z; y) + 32H(1; z)G(-z; y) + 48H(1; z)G(0, -z; y) 
-22H(1; z)G(0; y) - 12H(1; z)G(0, 0; y) - 36H(1, 0; z)G(l - z; y) - 48H(1, 0,1; z) + 36H(1,1,0; z) 

-48G(1 - z, -z, 1 - z; y) + 2G(1 - z; y) + 22G(1 - z, 0; y) + 12G(1 - z, 0, 0; y) - 12G(1 - z, 1,0; y) 
-32G(-z, 1 - z; y) + 22G(0,1 - z; y) + 12G(0,1 - z, 0; y) - 48G(0, -z, 1 - z; y) + 39G(0; y) 

+12G(0,0,1 - z; y) - 22G(0,0; y) - 24G(0,1,0; y) - 26G(1,0; y) - 24G(1, 0,0; y) + 24G(1,1, 0; y)] 

+ -—[ - eTT^Hfl; z) + GTr^Gfl - z; y) + 36H(0; z)G(l - z, 0; y) + 36H(0,1, 0; z) + 36H(1, 0; z) 

[y + zY 

+36H(1, 0; z)G(l - z;y) - 36H(1, 0; z)G(0;y) - 36H(1,1,0; z) - 36G(1 - z, 1, 0; y) - 36G(0,1,0;y) 

+36G(1,0; y)l + , ^ z) - Att'^GG - z;y) - 18H(0; z) - 24H(0; z)G(l - z, 0; y) 

[y + zy^ 

+36H(0; z)G(0; y) - 24H(0,1, 0; z) + 12H(1, 0; z) - 24H(1,0; z)G(l - z; y) + 24H(1,0; z)G(0; y) 
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+24H(1,1, 0; z) + 24G(1 - 0 ,1,0; y) + 18G(0; y) + 24G(0,1,0; y) - 60G(1, 0; y)] + —— [ - 

y + z 

27r2 

+6H(0; z) - 6H(0; z)G(0; y) - 6H(1, 0; z) - 6G(0; y) + 6G(1,0; y)] + 

27r^ 27r^ 

H(l; z) + —G(1 -z;y) + I 2 C 3 - 4H(0; z)G(0,1 - z;y) + 4H(0, 0,1; z) - 4H(0,1; z)G(l - z; y) 
-4H(0,1, 0; z) - 4H(1; z)G(l - z, -z; y) - 4H(1; z)G(0, -z; y) + 4H(1,0; z)G(l -z;y) + 4H(1,0,1; z) 

-4H(1,1, 0; z) + 4G(1 - z, -z, l-z;y) + 4G(0, -z, 1 - z; y)] + [ - 4H(0; z)G(l - z; y) 

-4H(1; z)G(-z; y) + 4G(-z, 1 - z; y)] + — [ - 2H(0; z)G(l - z-,y) - 2H(1; z)G(-z; y) 


+2G(-z, 1 - z; 2 /)l + --^ r + GTr^Hfl; z) - GTr^Gfl - z;y) - 36H(0; z)G(l - z, 0; y) - 36H(0,1,0; z) 

[y + zp 

-36H(1,0; z)G(l - z; y) + 36H(1, 0; z)G(0; y) + 36H(1,1,0; z) + 36G(1 - z, 1, 0; y) + 36G(0,1,0; y)] 
z^ 

+ 7 - -7 r — — 47 r^H(l; z) + 47 r^G(l — z;y) + 24H(0; z)G(l — z,0; y) — 36H(0; z)G(0; y) 

(y + z)'^ 

+24H(0,1,0; z) - 36H(1,0; z) + 24H(1,0; z)G(l - z-,y) - 24H(1,0; z)G(0; y) - 24H(1,1, 0; z) 

-24G(1 - z, 1,0; y) - 24G(0,1, 0; y) + 36G(1,0; y)] + [tt^ + 6H(0; z)G(0; y) + 6H(1, 0; z) 

[y + zy 

-6G(1,0; y)] + - - - -[ - ^ - 2H(0; z)G(0; y) - 2H(1,0; z) + 2G(1,0; y)] + , [ - 271^ 

Stt^ 47r^ 47r^ 27r^ Stt^ 

z) - —H(l; z) + —G(1 -z;y) + —G(0; y) + ^G(l; y) + - 8H(0; z)G(l - z, 0; y) 

-16H(0; z)G(0,1 - z; y) + 6H(0; z)G(0; y) + 4H(0; z)G(0,0; y) + 8H(0; z)G(l, 0; y) + 20H(0,0,1; z) 
+12H(0,1; z) - 20H(0,1; z)G(l - z;y) - 4H(0,1; z)G(0; y) - 8H(0,1, 0; z) - 20H(1; z)G(l - z, -z; y) 
+12H(1; z)G(-z; y) - 20H(1; z)G(0, -z; y) - 2H(1; z)G(0; y) + 4H(1; z)G(0, 0; y) - 10H(1,0; z) 

+8H(1, 0; z)G(l - z; y) + 8H(1, 0; z)G(0; y) + 20H(1,0,1; z) - 8H(1,1,0; z) + 20G(1 - z, -z, 1 - z; y) 
+2G(1 - z, 0; y) - 4G(1 - z, 0, 0; y) + 12G(1 - z, 1,0; y) - 12G(-z, 1 - z; y) + 2G(0,1 - z; y) 

-4G(0,1 - z, 0; y) + 20G(0, -z, 1 - z; y) - 5G(0; y) - 4G(0,0,1 - z; y) - 4G(0, 0; y) - 16G(1,1, 0; y)] 

+ 2 (y^z )2 y) - 2H(0; z)G(1 - z, 1 - z; y) + 3H(0; z)G(l - z; y) 

-44H(0; z)G(l - z, 0; y) - 2H(0; z)G(-z, 1 - z; y) + 2H(0; z)G(0,1 - z; y) + 3H(0,1; z) 

-2H(0,1; z)G(l - z; y) - 2H(0,1; z)G(-z; y) - 44H(0,1,0; z) + 2H(0,1,1; z) + 27H(1; z) 

-4H(1; z)G(l - z, -z; y) + 6H(1; z)G(l - z; y) + 2H(1; z)G(l - z, 0; y) - 4H(1; z)G(-z, 1 - z; y) 

-4H(1; z)G(-z, -z; y) + 6H(1; z)G(-z; y) + 2H(1; z)G(-z, 0; y) + 2H(1; z)G(0,1 - z; y) 

+2H(1; z)G(0, -z; y) - 3H(1; z)G(0; y) - 3H(1, 0; z) - 44H(1,0; z)G(l - z; y) + 2H(1,0; z)G(-z; y) 
+44H(1, 0; z)G(0; y) + 2H(1, 0,1; z) - 6H(1,1; z) + 4H(1,1; z)G(-z; y) - 2H(1,1; z)G(0; y) 

+44H(1,1,0; z) - 6G(1 - z, 1 - z; y) - 2G(1 - z, 1 - z, 0; y) + 4G(1 - z, -z, 1 - z; y) - 27G(1 - z; y) 
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—2G(1 — z, 0,1 — z]y) + 3G(1 — z, 0; y) + 46G(1 — z, 1,0; y) + 4G(—2, 1 — z,l — Z]y) — 6G{—z, 1 — z;y) 
-2G(-z, 1- z,0;y) + 4G(-2, -z, 1- z;y)- 2G(-z, 0,1- z;y)- 2G(0, 1- z,l- z;y) 

+3G(0, 1- z-,y)- 2G(0, - 2 :, 1- z;y) + 46G(0,1,0; y)] + 2(y\ z) ^ ~ ~ 2tj:'^Y1{1-, z) 

+27r2G(l - z; y) - 17 - 5H(0; z) - 4H(0; z)G(l - z; y) + 12H(0; z)G(l - z, 0; y) - 44H(0; z)G(0; y) 
-4H(0,1; z) + 12H(0,1, 0; z) + 4H(1; z) - 4H(1; z)G(l - z; y) - 8H(1; z)G(-z; y) + 4H(1; z)G(0; y) 
-42H(1, 0; z) + 12H(1, 0; z)G(l - z; y) - 12H(1, 0; z)G(0; y) + 4H(1,1; z) - 12H(1,1, 0; z) 

+4G(1 - z, 1 - z; y) - 4G(1 - z; y) - 4G(1 - z, 0; y) - 12G(1 - z, 1,0; y) + 8G(-z, 1 - z; y) 

Ttt^ 

-4G(0,1 - z; y) - 5G(0; y) - 12G(0,1,0; y) + 46G(1,0; y)] + — [ + 29 + 6H(1; z) 

+8H(1; z)G(l — z;y) — 16G(1 — z, 1 — z;y) + 12G(1 — z;y) + 8G(1 — z, 1;y) + 8G(0,1 — z;y) 

T 227r^ 2^^ 

-8G(0,1; y) - 18G(1; y) + 8G(1,1; y)] + - [ - ^ + _ - 6OC3 - 16C3H(1; z) - 16C3G(1 - z; y) 

+32C3G(1; y) + 18H(0; z)G(l - z, 1 - z; y) + 15H(0; z)G(l - z; y) - 8H(0; z)G(l - z, 0,1 - z; y) 

— 16H(0; z)G(—z, 1 — z, 1 — z; y) + 42H(0; z)G(—z, 1 — z; y) — 16H(0; z)G(—z, —z, 1 — z; y) 

+16H(0; z)G(-z, 0,1 - z; y) - 42H(0; z)G(0,1 - z; y) + 16H(0; z)G(0, -z, 1 - z; y) 

-16H(0; z)G(0,0,1 - z; y) + 16H(0, 0; z)G(l - z, 1 - z; y) - 16H(0, 0,1; z)G(l - z; y) 

+8H(0, 0,1; z)G(l; y) + 16H(0,0,1,1; z) + 15H(0,1; z) + 32H(0,1; z)G(l - z, 1 - z; y) 

+6H(0,1; z)G(l - z; y) + 8H(0,1; z)G(l - z, 0; y) - 16H(0,1; z)G(-z, 1 - z; y) 

-16H(0,1; z)G(-z, -z; y) + 42H(0,1; z)G(-z; y) + 8H(0,1; z)G(0,1 - z; y) - 8H(0,1; z)G(l, 1 - z; y) 
-8H(0,1; z)G(l, 0; y) - 8H(0,1, 0; z)G(l - z; y) + 16H(0,1,0,1; z) - 18H(0,1,1; z) 

+16H(0,1,1; z)G(-z; y) - 16H(0,1,1; z)G(0; y) + 16H(0,1,1, 0; z) - 52H(1; z) 

+32H(1; z)G(l - z, 1 - z, -z; y) + 24H(1; z)G(l - z, -z; y) - 18H(1; z)G(l - z; y) 

-18H(1; z)G(l - z, 0; y) - 16H(1; z)G(l - z, 0,0; y) - 8H(1; z)G(l - z, 1,0; y) 

—32H(1; z)G(—z, 1 — z, —z; y) + 36H(1; z)G(—z, 1 — z; y) + 16H(1; z)G(—z, 1 — z, 0; y) 

—32H(1; z)G(—z, —z, 1 — z; y) — 32H(1; z)G(—z, —z, —z; y) + 84H(1; z)G(—z, —z; y) 

+16H(1; z)G(—z, —z, 0; y) + 30H(1; z)G(—z; y) + 16H(1; z)G(—z, 0,1 — z; y) 

+16H(1; z)G(-z, 0, -z; y) - 42H(1; z)G(-z, 0; y) + 8H(1; z)G(0,1 - z, -z; y) 

-18H(1; z)G(0,1 - z; y) - 16H(1; z)G(0,1 - z, 0; y) + 16H(1; z)G(0, -z, 1 - z; y) 

+16H(1; z)G(0, -z, -z; y) - 42H(1; z)G(0, -z; y) - 15H(1; z)G(0; y) - 16H(1; z)G(0,0,1 - z; y) 

-16H(1; z)G(0,0, -z; y) + 8H(1; z)G(0,1, 0; y) - 8H(1; z)G(l, 1 - z, -z; y) - 8H(1; z)G(l, 0, -z; y) 
+16H(1; z)G(l, 0, 0; y) + 13H(1, 0; z) + 24H(1,0; z)G(l - z; y) + 16H(1,0; z)G(-z, -z; y) 

-42H(1, 0; z)G(-z; y) + 16H(1,0; z)G(0,1 - z; y) - 16H(1, 0; z)G(0, -z; y) + 8H(1,0, 0,1; z) 

-6H(1, 0,1; z) - 24H(1,0,1; z)G(l - z; y) + 16H(1,0,1; z)G(-z; y) - 8H(1,0,1; z)G(0; y) 

+8H(1, 0,1; z)G(l; y) + 18H(1,1; z) + 32H(1,1; z)G(-z, -z; y) - 36H(1,1; z)G(-z; y) 

-16H(1,1; z)G(-z, 0; y) - 16H(1,1; z)G(0, -z; y) + 18H(1,1; z)G(0; y) + 16H(1,1; z)G(0, 0; y) 
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+12H(1,1, 0; z) - 8H(1,1, 0; z)G(l -z;y) + 16H(1,1,0,1; z) + 16H(1,1,1, 0; z) 

—32G(1 — z,l — z, —z, 1 — z;y) + 18G(1 — z,l — z;y) + 18G(1 — z,l — z,0; y) 

+16G(1 - z, 1 - z, 0, 0; y) - 24G(1 - z, -z, 1 - z; y) + 52G(1 - z; y) + 18G(1 - z, 0,1 - z; y) 

+16G(1 - z, 0,1 - z, 0; y) + 15G(1 - z, 0; y) + 16G(1 - z, 0,0,1 - z; y) - 8G(1 - z, 0,1, 0; y) 

+8G(1 - z, 1,1 - z, 0; y) + 8G(1 - z, 1,0,1 - z; y) - 42G(1 - z, 1, 0; y) - 16G(1 - z, 1,0, 0; y) 

— 16G(1 — z, 1,1,0; y) — 36G(—z, 1 — z, 1 — z; y) — 16G(—z, 1 — z, 1 — z, 0; y) 

+32G(—z, 1 — z, —z, 1 — z; y) — 30G(—z, 1 — z; y) — 16G(—z, 1 — z, 0,1 — z; y) + 42G(—z, 1 — z, 0; y) 
+16G(—z, 1 — z, 1,0; y) + 32G(—z, —z, 1 — z, 1 — z; y) — 84G(—z, —z, 1 — z; y) 

— 16G(—z, —z, 1 — z, 0; y) + 32G(—z, —z, —z, 1 — z; y) — 16G(—z, —z, 0,1 — z; y) 

— 16G(—z, 0,1 — z, 1 — z; y) + 42G(—z, 0,1 — z; y) — 16G(—z, 0, —z, 1 — z; y) + 16G(—z, 0,1, 0; y) 
+18G(0,1 - z, 1 - z; y) + 16G(0,1 - z, 1 - z, 0; y) - 8G(0,1 - z, -z, 1 - z; y) + 15G(0,1 - z; y) 
+16G(0,1 - z, 0,1 - z; y) - 8G(0,1 - z, 1,0; y) - 16G(0, -z, 1 - z, 1 - z; y) + 42G(0, -z, 1 - z; y) 
-16G(0, -z, -z, 1 - z; y) + 16G(0,0,1 - z, 1 - z; y) + 16G(0, 0, -z, 1 - z; y) - 8G(0,1,1 - z, 0; y) 
-8G(0,1, 0,1 - z; y) - 42G(0,1,0; y) + 16G(0,1,1,0; y) + 8G(1,1 - z, -z, 1 - z; y) 

-16G(1,1 - z, 0,0; y) + 8G(1,1 - z, 1,0; y) - 16G(1, 0,1 - z, 0; y) + 8G(1,0, -z, 1 - z; y) - 28G(1, 0; y) 
-16G(1,0,0,1 - z; y) + 8G(1,0,1,0; y) + 36G(1,1,0; y) + 16G(1,1,0,0; y) - 16G(1,1,1,0; y)] 

TT^ 

+ — [9 - 7H(0; z) + 6H(0; z)G(l - z; y) - 2H(0; z)G(l; y) + 4H(0,1; z) - 24H(1; z) 

6 

+12H(1; z)G(l - z; y) + 8H(1; z)G(-z; y) - 4H(1; z)G(0; y) - 6H(1; z)G(l; y) - 2H(1,0; z) - 6H(1,1; z) 
-20G(1 - z, 1 - z; y) + 34G(1 - z; y) + 6G(1 - z, 0; y) + 8G(1 - z, 1; y) - 8G(-z, 1 - z; y) 

+8G(0,1 - z; y) - 7G(0; y) - 4G(0,1; y) + 6G(1,1 - z; y) - 10G(1; y) - 4G(1,0; y)] 

+ i [40C3 + 40C3H(1; z) - 80C3G(1 - z; y) + 40C3G(1; y) - 2 - 29H(0; z) + 20H(0; z)G(l - z, 1 - z; y) 

-32H(0; z)G(l - z, 1 - z, 0; y) - 4H(0; z)G(l - z; y) + 32H(0; z)G(l - z, 0; y) 

+16H(0; z)G(l - z, 0,0; y) + 16H(0; z)G(l - z, 1,0; y) - 24H(0; z)G(-z, 1 - z, 1 - z; y) 

+52H(0; z)G(-z, 1 - z; y) - 32H(0; z)G(-z, 1 - z, 0; y) - 24H(0; z)G(-z, -z, 1 - z; y) 

+24H(0; z)G(-z, 0,1 - z; y) + 8H(0; z)G(0,1 - z, 1 - z; y) - 44H(0; z)G(0,1 - z; y) 

+16H(0; z)G(0,1 - z, 0; y) + 24H(0; z)G(0, -z, 1 - z; y) + 40H(0; z)G(0; y) - 24H(0; z)G(0, 0,1 - z; y) 
-20H(0; z)G(0, 0; y) + 16H(0; z)G(l, 1 - z, 0; y) - 8H(0; z)G(l, 0,1 - z; y) + 4H(0; z)G(l, 0; y) 

-16H(0; z)G(l, 0,0; y) - 16H(0; z)G(l, 1,0; y) + 20H(0, 0; z) + 16H(0, 0; z)G(l - z, 1 - z; y) 

-4H(0,0; z)G(l - z; y) + 16H(0,0; z)G(l - z, 0; y) - 20H(0, 0; z)G(0; y) + 36H(0,0,1; z) 

-16H(0,0,1; z)G(l - z; y) + 16H(0,0,1; z)G(l; y) - 16H(0,0,1, 0; z) + 16H(0,0,1,1; z) 

+48H(0,1; z)G(l - z, 1 - z; y) - 52H(0,1; z)G(l - z; y) - 24H(0,1; z)G(-z, 1 - z; y) 

-24H(0,1; z)G(-z, -z; y) + 84H(0,1; z)G(-z; y) + 8H(0,1; z)G(0,1 - z; y) - 8H(0,1; z)G(0; y) 
-16H(0,1; z)G(l, 1 - z; y) - 8H(0,1; z)G(l, 0; y) - 28H(0,1,0; z) + 32H(0,1,0; z)G(l - z; y) 
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-32H(0,1, 0; z)G{-z-, y) + 16H(0,1, 0; z)G(0; y) - 24H(0,1, 0; z)G(l; y) + 16H(0,1,0,1; z) 

-20H(0,1,1; z) + 24H(0,1,1; z)G{-z; y) - 16H(0,1,1; z)G(0; y) + 32H(0,1,1, 0; z) - 58H(1; z) 
+48H(1; z)G(l — z,l — z, —z; y) — 32H(1; z)G(l — z, —z] y) — 4H(1; z)G(l — z; y) 

-20H(1; z)G(l - z, 0; y) - 16H(1; z)G(l - z, 0,0; y) - 8H(1; z)G(l - z, 1,0; y) 

—48H(1; z)G(—z, 1 — z, —z; y) + 40H(1; z)G(—z, 1 — z; y) + 24H(1; z)G(—z, 1 — z, 0; y) 

—48H(1; z)G(—z, —z, 1 — z; y) — 48H(1; z)G(—z, —z, —z; y) + 136H(1; z)G(—z, —z; y) 

+24H(1; z)G(-z, -z, 0; y) - 4H(1; z)G(-z; y) + 24H(1; z)G(-z, 0,1 - z; y) + 24H(1; z)G(-z, 0, -z; y) 
-52H(1; z)G(-z, 0; y) + 16H(1; z)G(0,1 - z, -z; y) - 20H(1; z)G(0,1 - z; y) 

-16H(1; z)G(0,1 - z, 0; y) + 24H(1; z)G(0, -z, 1 - z; y) + 24H(1; z)G(0, -z, -z; y) 

-52H(1; z)G(0, -z; y) + 4H(1; z)G(0; y) - 16H(1; z)G(0, 0,1 - z; y) - 24H(1; z)G(0,0, -z; y) 

+4H(1; z)G(0,0; y) + 8H(1; z)G(0,1, 0; y) - 16H(1; z)G(l, 1 - z, -z; y) - 16H(1; z)G(l, 0, -z; y) 
+24H(1; z)G(l, 0; y) + 16H(1; z)G(l, 0, 0; y) + 60H(1, 0; z) - 40H(1,0; z)G(l - z, 1 - z; y) 

+80H(1, 0; z)G(l - z; y) + 16H(1, 0; z)G(l - z, 0; y) - 32H(1, 0; z)G(-z, 1 - z; y) 

+24H(1, 0; z)G(-z, -z; y) - 52H(1,0; z)G(-z; y) + 32H(1, 0; z)G(-z, 0; y) + 32H(1, 0; z)G(0,1 - z; y) 
-24H(1, 0; z)G(0, -z; y) - 36H(1, 0; z)G(0; y) - 16H(1, 0; z)G(0, 0; y) + 24H(1, 0; z)G(l, 1 - z; y) 
-16H(1, 0; z)G(l, 0; y) - 20H(1, 0, 0; z) + 16H(1,0, 0; z)G(l - z; y) + 16H(1,0, 0,1; z) + 28H(1,0,1; z) 
-40H(1,0,1; z)G(l - z; y) + 24H(1, 0,1; z)G(-z; y) - 8H(1, 0,1; z)G(0; y) + 16H(1,0,1; z)G(l; y) 
-24H(1, 0,1,0; z) + 4H(1,1; z) + 48H(1,1; z)G(-z, -z; y) - 40H(1,1; z)G(-z; y) 

-24H(1,1; z)G(-z, 0; y) - 24H(1,1; z)G(0, -z; y) + 20H(1,1; z)G(0; y) + 16H(1,1; z)G(0, 0; y) 

-64H(1,1,0; z) + 16H(1,1,0; z)G(l - z; y) + 32H(1,1, 0; z)G(-z; y) - 16H(1,1, 0; z)G(0; y) 

-24H(1,1,0; z)G(l; y) + 24H(1,1,0,1; z) - 8H(1,1,1, 0; z) - 48G(1 - z, 1 - z, -z, 1 - z; y) 

+4G(1 - z, 1 - z; y) + 20G(1 - z, 1 - z, 0; y) + 16G(1 - z, 1 - z, 0, 0; y) + 40G(1 - z, 1 - z, 1,0; y) 
+32G(1 - z, -z, 1 - z; y) + 58G(1 - z; y) + 20G(1 - z, 0,1 - z; y) + 16G(1 - z, 0,1 - z, 0; y) 

-4G(1 - z, 0; y) + 16G(1 - z, 0, 0,1 - z; y) - 4G(1 - z, 0,0; y) - 48G(1 - z, 0,1,0; y) 

+8G(1 - z, 1,1 - z, 0; y) + 8G(1 - z, 1,0,1 - z; y) - 100G(1 - z, 1, 0; y) - 32G(1 - z, 1, 0,0; y) 

-32G(1 - z, 1,1, 0; y) - 40G(-z, 1 - z, 1 - z; y) - 24G(-z, 1 - z, 1 - z, 0; y) 

+48G(—z, 1 — z, —z, 1 — z; y) + 4G(—z, 1 — z; y) — 24G(—z, 1 — z, 0,1 — z; y) 

+52G(—z, 1 — z, 0; y) + 56G(—z, 1 — z, 1,0; y) + 48G(—z, —z, 1 — z, 1 — z; y) 

— 136G(—z, —z, 1 — z; y) — 24G(—z, —z, 1 — z, 0; y) + 48G(—z, —z, —z, 1 — z; y) 

—24G(—z, —z, 0,1 — z; y) — 24G(—z, 0,1 — z, 1 — z; y) + 52G(—z, 0,1 — z; y) — 24G(—z, 0, —z, 1 — z; y) 
+56G(-z, 0,1, 0; y) + 20G(0,1 - z, 1 - z; y) + 16G(0,1 - z, 1 - z, 0; y) - 16G(0,1 - z, -z, 1 - z; y) 
-4G(0,1 - z; y) + 16G(0,1 - z, 0,1 - z; y) - 4G(0,1 - z, 0; y) - 24G(0,1 - z, 1,0; y) 

-24G(0, -z, 1 - z, 1 - z; y) + 52G(0, -z, 1 - z; y) - 24G(0, -z, -z, 1 - z; y) - 29G(0; y) 

+16G(0,0,1 - z, 1 - z; y) - 4G(0, 0,1 - z; y) + 24G(0,0, -z, 1 - z; y) + 20G(0, 0; y) - 16G(0, 0,1, 0; y) 
-8G(0,1,1 - z, 0; y) - 8G(0,1,0,1 - z; y) - 20G(0,1,0; y) + 16G(0,1,1,0; y) 
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+16G(1, z; y) - 24G(1,1 - z, 0; y) - 16G(1,1 - z, 0,0; y) - 8G(1,1 - z, 1, 0; y) 

-24G(1, 0,1 - z; y) - 16G(1,0,1 - z, 0; y) + 16G(1,0, -z, 1 - z; y) - 56G(1,0; y) - 16G(1, 0,0,1 - z; y) 
+24G(1, 0,0; y) + 24G(1,0,1,0; y) + 40G(1,1,0; y) + 32G(1,1, 0,0; y)] , (4.20) 

D2o{y,z) = 

^ [ - 2H(0; z) - 3H(1, 0; z) - 3G(1, 0; y)] + " 2H(1,0; z) - 2G(1,0; y)] + ^ [ - 74 

+15H(0; z) + 36H(1,0; z) + 15G(0; y) + 36G(1,0; y)] + y ^2 “ 3H(0; z)G(0; y) + 50G(0; y) 

-18G(0, 0; y) - 12G(1,0; y)] + [Stt" + 38 - 3H(0; z) - 9H(0; z)G(0; y) + 87G(0; y) 

-54G(0,0;y) - 36G(1,0;y)] + , f .3 [ - 2Tr^R{l; z) + 27r2G(l - z; y) + 12H(0; z)G(l - z,0;y) 

(y + zr 

+12H(0,1, 0; z) + 12H(1,0; z) + 12H(1,0; z)G(l - z-,y) - 12H(1,0; z)G(0; y) - 12H(1,1, 0; z) 

-12G(1 - z, 1,0; y) - 12G(0,1, 0; y) + 12G(1, 0; y)] + [2^^ + ^H(l; z) - ^G(l - z; y) 

-6H(0; z) - 8H(0; z)G(l - z, 0; y) + 12H(0; z)G(0; y) - 8H(0,1,0; z) + 4H(1,0; z) 

-8H(1, 0; z)G(l - z; y) + 8H(1, 0; z)G(0; y) + 8H(1,1,0; z) + 8G(1 - z, 1, 0; y) + 6G(0; y) 

2 

+8G(0,1,0; y) - 20G(1, 0; y)] + ^ [ - ^ + 2H(0; z) - 2H(0; z)G(0; y) - 2H(1,0; z) - 2G(0; y) 

+2G(1,0; y)] + [2^"H(1; z) - 27r2G(l - z; y) - 12H(0; z)G(l - z, 0; y) - 12H(0,1, 0; z) 

-12H(1, 0; z)G(l - z; y) + 12H(1,0; z)G(0; y) + 12H(1,1, 0; z) + 12G(1 - z, 1, 0; y) + 12G(0,1,0; y)] 

4:7r^ 47r^ 

^ - 2^2 _ H(l; z) + —G(1 - z; y) + 8H(0; z)G(l - z, 0; y) - 12H(0; z)G(0; y) 

(y + z)'="- 3 3 

+8H(0,1,0; z) - 12H(1,0; z) + 8H(1,0; z)G(l - z; y) - 8H(1,0; z)G(0; y) - 8H(1,1,0; z) 

-8G(1 - z, 1,0; y) - 8G(0,1,0; y) + 12G(1,0; y)] + + 6H(0; z)G(0; y) + 6H(1, 0; z) 

-6G(1, 0; y)] + + 6H(0; z)G(0; y) - 70G(0; y) + 36G(0, 0; y) + 24G(1,0; y)] 

1 r 27r^ 27r^ 

+ (^y + ^)2 J/) - 4H(0; z)G(l - z, 0; y) - 4H(0,1, 0; z) - 4H(1,0; z)G(l - z; y) 

1 tt^ 

+4H(1, 0; z)G(0; y) + 4H(1,1,0; z) + 4G(1 - z, 1,0; y) + 4G(0,1,0; y)] + ^ ^ - ^H(l; z) 

y z o o 

TT^ 

+—G(1 - z; y) + 2 - H(0; z) + 2H(0; z)G(l - z, 0; y) - 4H(0; z)G(0; y) + 2H(0,1, 0; z) - 4H(1, 0; z) 
o 

+2H(1, 0; z)G(l - z; y) - 2H(1, 0; z)G(0; y) - 2H(1,1,0; z) - 2G(1 - z, 1, 0; y) - G(0; y) - 2G(0,1, 0; y) 

+4G(1, 0; y)] + ^ [ + 431 - 12H(0; z) + 24G(1 - z; y) - 12G(0; y) - 24G(1; y)] + ^ - 38 C 3 

+31H(0; z) + 12H(0; z)G(l - z, 0; y) + 10H(0; z)G(0; y) - 18H(0; z)G(0, 0; y) + 3H(0; z)G(l, 0; y) 
-41H(0,0; z) - 18H(0, 0; z)G(0; y) - 3H(0,1,0; z) + 29H(1,0; z) + 12H(1, 0; z)G(l - z; y) 
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-15H(1, 0; 2 )G( 0 ; y) - 18H(1, 0,0; z) - 12G(1 - z, 1, 0; y) + 31G(0; y) - 41G(0,0; y) + 3G(0,1, 0; y) 
-29G(1,0; y) + 18G(1,0, 0; y) + 12G(1,1,0; 2/)] + y + 2H(0; z)G(0; y) + 2H(1, 0; z) - 2G(1, 0; y) , (4.21) 
E2o{y,z) = 

A [ _ 27r2G(l -z;y) + 12H(0; z)G(l - z, 1 - z; y) - 47H(0; z)G(l -z;y) + 3H(0; z)G(l - z, 0; y) 

9y 

-15H(0; z)G(-z, 1 - z; y) + 3H(0; z)G(0,1 - z; y) + 18H(0,0; z)G(l - z; y) + 9H(0,1; z) 

+3H(0,1; z)G(-z; y) + 12H(1; z)G(l - z, -z; y) + 12H(1; z)G(-z, 1 - z; y) - 12H(1; z)G(-z, -z; y) 
-38H(1; z)G(-z; y) + 3H(1; z)G(-z, 0; y) + 3H(1; z)G(0, -z; y) - 9H(1; z)G(0; y) + 9H(1,0; z) 

-12H(1, 0; z)G(l - z; y) + 15H(1,0; z)G(-z; y) - 12H(1,1; z)G(-z; y) - 12G(1 - z, -z, 1 - z; y) 

+9G(1 — z, 0; y) — 12G(—z, 1 — z, 1 — z; y) + 38G(—z, 1 — z; y) — 3G(—z, 1 — z, 0; y) 

+12G(-z, -z, 1 - z; y) - 3G(-z, 0,1 - z; y) + 9G(0,1 - z; y) - 3G(0, -z, 1 - z; y)] 

+ , [ - 2H(0; z)G(l - z; y) + 2H(0,1; z) - 2H(1; z)G(0; y) + 2H(1, 0; z) + 2G(1 - z, 0; y) 

y(y + z) 

+2G(0,1 - z; y)l + — r 87 r 2 G(l - z; y) - 38 + 3H(0; z) - 48H(0; z)G(l - z, 1 - z; y) 

18y 

+188H(0; z)G(l -z;y)- 12H(0; z)G(l - z, 0; y) + 60H(0; z)G(-z, 1 - z;y) - 12H(0; z)G(0,1 - z;y) 
-72H(0,0; z)G(l - z; y) - 36H(0,1; z) - 12H(0,1; z)G(-z; y) - 48H(1; z)G(l - z, -z; y) 

-48H(1; z)G(-z, 1 - z; y) + 48H(1; z)G(-z, -z; y) + 152H(1; z)G(-z; y) - 12H(1; z)G(-z, 0; y) 

-12H(1; z)G(0, -z; y) + 36H(1; z)G(0; y) - 36H(1,0; z) + 48H(1,0; z)G(l - z; y) 

-60H(1, 0; z)G(-z; y) + 48H(1,1; z)G(-z; y) + 48G(1 - z, -z, 1 - z; y) - 36G(1 - z, 0; y) 

+48G(—z, 1 — z, 1 — z; y) — 152G(—z, 1 — z; y) + 12G(—z, 1 — z, 0; y) — 48G(—z, —z, 1 — z; y) 

+12G(-z, 0,1 - z; y) - 36G(0,1 - z; y) + 12G(0, -z, 1 - z; y) + 15G(0; y)] + [ - 2^' 

+3H(0; z)G(0; y) - 50G(0; y) + 18G(0,0; y) + 12G(1, 0; y)] + - 38 + 3H(0; z) 

+9H(0; z)G(0; y) - 87G(0; y) + 54G(0,0; y) + 36G(1, 0; y)] + 

+38G(0; y) - 18G(0,0; y) - 12G(1,0; y)] + " 9H(0; z)G(l - z; y) - 3H(0,1; z) - 26H(1; z) 

+12H(1; z)G(l -z;y)- 12H(1; z)G(-z; y) + 9H(1; z)G(0; y) + 9H(1, 0; z) - 12H(1,1; z) 

-12G(1 - z, 1 - z; y) + 26G(1 - z; y) - 9G(1 - z, 0; y) + 12G(-z, 1 - z; y) - 9G(0,1 - z; y)] 

1 

+ + 12G(1 - z; y) - 9G(0; y)] + ^ [ - 7 + H(l; z) 

-5G(1 - z; y) + 4G(1; y)] + ^ [ - + 6 C 3 + 72H(0; z) + 72H(0; z)G(l - z, 1 - z; y) 

-147H(0; z)G(l - z; y) + 36H(0; z)G(l - z, 0; y) - 108H(0; z)G(-z, 1 - z; y) + 36H(0; z)G(0,1 - z; y) 
-18H(0; z)G(l, 0; y) + 108H(0,0; z)G(l - z; y) - 36H(0,0,1; z) - 201H(0,1; z) 

+72H(0,1; z)G(l - z; y) - 36H(0,1; z)G(-z; y) + 72H(0,1; z)G(0; y) + 18H(0,1, 0; z) - 72H(0,1,1; z) 
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+68H(1; z) + 144H(1; z)G(l - 2 , -z; y) - 108H(1; z)G(l -z;y)- 72H(1; z)G(l - z,0; y) 

+144H(1; z)G(-z, 1 - z; y) - 144H(1; z)G(-z, -z; y) - 348H(1; z)G(-z; y) + 108H(1; z)G(-z, 0; y) 
-72H(1; z)G(0,1 - z; y) + 108H(1; z)G(0, -z; y) + 147H(1; z)G(0; y) - 108H(1; z)G(0,0; y) 

-81H(1, 0; z) - 72H(1, 0; z)G(l - z; y) + 108H(1, 0; z)G(-z; y) - 18H(1,0; z)G(0; y) - 72H(1,0,1; z) 
+108H(1,1; z) - 144H(1,1; z)G(-z; y) + 72H(1,1; z)G(0; y) + 108G(1 - z, 1 - z; y) 

+72G(1 - z, 1 - z, 0; y) - 144G(1 - z, -z, 1 - z; y) - 68G(1 - z; y) + 72G(1 - z, 0,1 - z; y) 

-147G(1 - z, 0; y) + 108G(1 - z, 0, 0; y) - 144G(-z, 1 - z, 1 - z; y) + 348G(-z, 1 - z; y) 

-108G(-z, 1 - z, 0; y) + 144G(-z, -z, 1 - z; y) - 108G(-z, 0,1 - z; y) + 72G(0,1 - z, 1 - z; y) 
-147G(0,1 - z; y) + 108G(0,1 - z, 0; y) - 108G(0, -z, 1 - z; y) + 72G(0; y) + 108G(0,0,1 - z; y) 

-18G(0,1,0; y) + 228G(1,0; y) - 108G(1, 0,0; y) - 72G(1,1, 0; y)] + ^ [ 2 ^^ ^ 27r"H(l; z) 

- z; y) + y) + 19H(0; z) + 12H(0; z)G(l - z, 1 - z; y) - 29H(0; z)G(l - z; y) 

+6H(0; z)G(l - z, 0; y) - 18H(0; z)G(-z, 1 - z; y) + 6H(0; z)G(0,1 - z; y) - 3H(0; z)G(0; y) 

-3H(0; z)G(l, 0; y) - 9H(0, 0; z) + 18H(0,0; z)G(l - z; y) - 6H(0, 0,1; z) - 35H(0,1; z) 

+12H(0,1; z)G(l - z; y) - 6H(0,1; z)G(-z; y) + 12H(0,1; z)G(0; y) + 3H(0,1,0; z) - 12H(0,1,1; z) 
+38H(1; z) + 24H(1; z)G(l - z, -z; y) - 12H(1; z)G(l - z; y) - 12H(1; z)G(l - z, 0; y) 

+24H(1; z)G(-z, 1 - z; y) - 24H(1; z)G(-z, -z; y) - 64H(1; z)G(-z; y) + 18H(1; z)G(-z, 0; y) 

-12H(1; z)G(0,1 - z; y) + 18H(1; z)G(0, -z; y) + 29H(1; z)G(0; y) - 18H(1; z)G(0, 0; y) - 3H(1,0; z) 
-12H(1, 0; z)G(l - z; y) + 18H(1,0; z)G(-z; y) - 3H(1,0; z)G(0; y) - 12H(1, 0,1; z) + 12H(1,1; z) 
-24H(1,1; z)G(-z; y) + 12H(1,1; z)G(0; y) + 12G(1 - z, 1 - z; y) + 12G(1 - z, 1 - z, 0; y) 

-24G(1 - z, -z, 1 - z; y) - 38G(1 - z; y) + 12G(1 - z, 0,1 - z; y) - 29G(1 - z, 0; y) 

+18G(1 — z, 0,0; y) — 24G(—z, 1 — z, 1 — z; y) + 64G(—z, 1 — z; y) — 18G(—z, 1 — z, 0; y) 

+24G(-z, -z, 1 - z; y) - 18G(-z, 0,1 - z; y) + 12G(0,1 - z, 1 - z; y) - 29G(0,1 - z; y) 

+18G(0,1 - z, 0; y) - 18G(0, -z, 1 - z; y) + 19G(0; y) + 18G(0, 0,1 - z; y) - 9G(0,0; y) 

-3G(0,1, 0; y) + 32G(1, 0; y) - 18G(1,0, 0; y) - 12G(1,1,0; y)] , (4.22) 

F2o{y,z) = 

^ [ - 177r2 - 20H(0; z) + 3H(0; z)G(0; y) + 15H(0, 0; z) - 20G(0; y) + 15G(0,0; y)] , (4.23) 

G2oiy,z) = 

2 

^ [ - 9H(0; z)G(l - z; y) - 9H(1; z)G(-z; y) + 9G(-z, 1 - z; y)] + ^ [3H(0; z)G(l - z; y) 

+3H(1; z)G(-z; y) - 3G(-z, 1 - z; y)] + ^ [6H(0; z)G(l - z; y) + 6H(1; z)G(-z; y) - 6G(-z, 1 - z; y)] 

2 

ZTT^ 

+— [12H(1; z)G(l - z; y) - 12G(1 - z, 1 - z; y) - 2G(1 - z; y) + 12G(1 - z, 0; y) - 2G(1; y)] 

+ ^ [ - 72C3G(1 - z; y) - 9H(0; z) + 12H(0; z)G(l - z, 1 - z; y) - 18H(0; z)G(l - z; y) 

3y 
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-4H(0; z)G(l - z, 0; y) - 8H(0; z)G{-z, l-z;y)- 4H(0; z)G(l,0; y) + 24H(0, 0,1; z)G(l - z; y) 
+24H(0,1; z)G(l - z, 1 - z; y) - 8H(0,1; z)G(l - z;y) - 24H(0,1; z)G(l - z, 0; y) 

-24H(0,1; z)G(-z, 1 - z; y) - 24H(0,1; z)G(-z, -z; y) + 4H(0,1,1; z) + 24H(0,1,1; z)G(-z; y) 
+24H(1; z)G(l - z, 1 - z, -z; y) + 4H(1; z)G(l - z, -z; y) - 24H(1; z)G(l - z, 0, -z; y) 

+4H(1; z)G(l - z, 0; y) - 24H(1; z)G(-z, 1 - z, -z; y) + 8H(1; z)G(-z, 1 - z; y) 

+24H(1; z)G(-z, 1 - z, 0; y) - 24H(1; z)G(-z, -z, 1 - z; y) - 24H(1; z)G(-z, -z, -z; y) 

-8H(1; z)G(-z, -z; y) + 24H(1; z)G(-z, -z, 0; y) - 18H(1; z)G(-z; y) + 24H(1; z)G(-z, 0,1 - z; y) 
-8H(1; z)G(-z, 0; y) + 4H(1; z)G(0,1 - z; y) - 4H(1; z)G(l, 0; y) - 12H(1,0; z)G(l - z; y) 

+8H(1,0; z)G(-z; y) + 4H(1,0,1; z) - 24H(1,0,1; z)G(l - z; y) + 24H(1,0,1; z)G(-z; y) 

+24H(1,1; z)G(-z, -z; y) - 8H(1,1; z)G(-z; y) - 24H(1,1; z)G(-z, 0; y) - 4H(1,1; z)G(0; y) 
+4H(1,1, 0; z) - 24G(1 - z, 1 - z, -z, 1 - z; y) - 4G(1 - z, 1 - z, 0; y) + 24G(1 - z, 1 - z, 1,0; y) 
-4G(1 - z, -z, 1 - z; y) - 4G(1 - z, 0,1 - z; y) + 24G(1 - z, 0, -z, 1 - z; y) - 24G(1 - z, 0,1, 0; y) 
—8G(—z, 1 — z, 1 — z; y) — 24G(—z, 1 — z, 1 — z, 0; y) + 24G(—z, 1 — z, —z, 1 — z; y) 

+18G(—z, 1 — z; y) — 24G(—z, 1 — z, 0,1 — z; y) + 8G(—z, 1 — z, 0; y) + 24G(—z, 1 — z, 1,0; y) 
+24G(—z, —z, 1 — z, 1 — z; y) + 8G(—z, —z, 1 — z; y) — 24G(—z, —z, 1 — z, 0; y) 


+24G(—z, —z, —z, 1 — z; y) — 24G(—z, —z, 0,1 — z; y) — 24G(—z, 0,1 — z, 1 — z; y) + 8G(—z, 0,1 — z; y) 
+24G(-z, 0,1, 0; y) - 4G(0,1 - z, 1 - z; y) + 4G(1,1 - z, 0; y) + 4G(1, 0,1 - z; y) + 4G(1,1, 0; y)] 

r2 1 27r^ 

+~ [3H(0; z)G(l - z; y) + 3H(1; z)G(-z; y) - 3G(-z, 1 - z; y)] + ^ [~G(1; y) 

+4H(0; z)G(-z, 1 - z; y) + 4H(0; z)G(l, 0; y) + 4H(1; z)G(-z, -z; y) + 4H(1; z)G(-z, 0; y) 

+4H(1; z)G(l, 0; y) - 4H(1, 0; z)G(-z; y) - 4G(-z, 1 - z, 0; y) - 4G(-z, -z, 1 - z; y) 

1 27r^ 

-4G(-z, 0,1 - z; y) - 4G(1,1 - z, 0; y) - 4G(1, 0,1 - z; y) - 4G(1,1,0; y)] + ——-- [—G(l; y) 

3y(l - z) 3 

+4H(0; z)G(-z, 1 - z; y) + 4H(0; z)G(l, 0; y) + 4H(1; z)G(-z, -z; y) + 4H(1; z)G(-z, 0; y) 

+4H(1; z)G(l, 0; y) - 4H(1, 0; z)G(-z; y) - 4G(-z, 1 - z, 0; y) - 4G(-z, -z, 1 - z; y) 

-4G(-z, 0,1 - z; y) - 4G(1,1 - z, 0; y) - 4G(1, 0,1 - z; y) - 4G(1,1,0; y)] 


1 


■ [8H(0; z)G(l - z, 1 - z; y) - 8H(0,1; z)G(l - z; y) + 8H(0,1,1; z) + 8H(1; z)G(l - z, 0; y) 


3y(y + z) 

+8H(1; z)G(0,1 - z; y) - 8H(1, 0; z)G(l - z; y) + 8H(1, 0,1; z) - 8H(1,1; z)G(0; y) + 8H(1,1,0; z) 
-8G(1 - z, 1 - z, 0; y) - 8G(1 - z, 0,1 - z; y) - 8G(0,1 - z, 1 - z; y)] + ^ [ - 6H(1; z)G(l - z; y) 

+6G(1 - z, 1 - z; y) + G(1 - z; y) - 6G(1 - z, 0; y) + G(l; y)] + ^ [36C3G(1 - z; y) + 18H(0; z) 

-12H(0; z)G(l - z, 1 - z; y) + 11H(0; z)G(l - z; y) - 4H(0; z)G(l - z, 0; y) - 2H(0; z)G(-z, 1 - z; y) 
+6H(0; z)G(0,1 - z; y) + 2H(0; z)G(l, 0; y) - 12H(0, 0,1; z)G(l - z; y) - 12H(0,1; z)G(l - z, 1 - z; y) 
+4H(0,1; z)G(l - z; y) + 12H(0,1; z)G(l - z, 0; y) + 12H(0,1; z)G(-z, 1 - z; y) 
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+12H(0,1; z)G{-z, -z] y) + 6H(0,1; z)G{-z] y) - 8H(0,1,1; z) - 12H(0,1,1; z)G{-z-, y) + 9H(1; z) 
-12H(1; z)G(l -z,l-z,-z-y) - 8H(1; 2)G(1 - z, -z; y) + 12H(1; z)G(l - z, 0, -z; y) 

-8H(1; z)G(l - z, 0; y) + 12H(1; z)G(-z, 1 - z, -z; y) - 4H(1; z)G(-z, 1 - z; y) 

— 12H(1; z)G(—z, 1 — z, 0; y) + 12H(1; z)G(—z, —z, 1 — z; y) + 12H(1; z)G(—z, —z, —z; y) 

+4H(1; z)G(-z, -z; y) - 12H(1; z)G(-z, -z, 0; y) + 11H(1; z)G(-z; y) - 12H(1; z)G(-z, 0,1 - z; y) 
-2H(1; z)G(-z, 0; y) - 8H(1; z)G(0,1 - z; y) + 6H(1; z)G(0, -z; y) + 2H(1; z)G(l, 0; y) 

+6H(1, 0; z)G(l - z; y) + 2H(1, 0; z)G(-z; y) - 8H(1, 0,1; z) + 12H(1,0,1; z)G(l - z; y) 

-12H(1, 0,1; z)G(-z; y) - 12H(1,1; z)G(-z, -z; y) + 4H(1,1; z)G(-z; y) + 12H(1,1; z)G(-z, 0; y) 
+8H(1,1; z)G(0; y) - 8H(1,1,0; z) + 12G(1 - z, 1 - z, -z, 1 - z; y) + 8G(1 - z, 1 - z, 0; y) 

— 12G(1 — z, 1 — z, 1,0; y) + 8G(1 — z, —z, 1 — z; y) — 9G(1 — z; y) + 8G(1 — z, 0,1 — z; y) 

— 12G(1 — z, 0, —z, 1 — z; y) + 12G(1 — z, 0,1, 0; y) + 4G(—z, 1 — z, 1 — z; y) 

+12G(—z, 1 — z, 1 — z, 0; y) — 12G(—z, 1 — z, —z, 1 — z; y) — 11G(—z, 1 — z; y) 

+12G(-z, 1 - z, 0,1 - z; y) + 2G(-z, 1 - z, 0; y) - 12G(-z, 1 - z, 1,0; y) 

— 12G(—z, —z, 1 — z, 1 — z; y) — 4G(—z, —z, 1 — z; y) + 12G(—z, —z, 1 — z, 0; y) 


— 12G(—z, —z, —z, 1 — z; y) + 12G(—z, —z, 0,1 — z; y) + 12G(—z, 0,1 — z, 1 — z; y) 

+2G(-z, 0,1 - z; y) - 12G(-z, 0,1, 0; y) + 8G(0,1 - z, 1 - z; y) - 6G(0, -z, 1 - z; y) 

z 2'k’^ 27r^ 

-2G(1,1 - z, 0; y) - 2G(1, 0,1 - z; y) - 2G(1,1,0; y)] + —-- [ - —H(l; z) - —G(l; y) 

3(1 — y — z) 3 3 

-8H(0; z)G(-z, 1 - z; y) + 4H(0; z)G(0,1 - z; y) - 4H(0; z)G(l, 0; y) + 8H(0, 0,1; z) 

+8H(0,1; z)G(-z; y) - 4H(0,1; z)G(0; y) + 4H(0,1,0; z) - 8H(1; z)G(-z, 0; y) - 4H(1; z)G(l, 0; y) 
+8H(1, 0; z)G(-z; y) - 4H(1,0; z)G(0; y) - 4H(1, 0,1; z) - 8H(1,1, 0; z) + 8G(-z, 1 - z, 0; y) 

+8G(-z, 0,1 - z; y) - 4G(0,1, 0; y) + 4G(1,1 - z, 0; y) + 4G(1, 0,1 - z; y) + 4G(1,1, 0; y)] 

2 

+^(r^ [ - 3 - H(l; z) + G(1 - z; y)] + [SCs - 2H(0; z)G(0,1 - z; y) 

+2H(0, 0,1; z) + 6H(0,1; z) - 2H(0,1; z)G(l - z; y) - 2H(0,1,0; z) - 2H(1; z)G(l - z, -z; y) 

+6H(1; z)G(-z; y) - 2H(1; z)G(0, -z; y) - 6H(1,0; z) + 2H(1, 0; z)G(l - z; y) + 2H(1, 0,1; z) 

-2H(1,1,0; z) + 2G(1 - z, -z, 1 - z; y) - 6G(-z, 1 - z; y) + 2G(0, -z, 1 - z; y)] 

+ -^ [3H(0; z) - 2H(0; z)G(l - z; y) + 3H(1; z) - 2H(1; z)G(-z; y) - 3G(1 - z; y) + 2G(-z, 1 - z; y)] 
1 - y 


■ ‘ - 8H(0; z)G(l - z, 1 - z; y) - 24H(0; z)G(l - z; y) + 24H(0,1; z) + 8H(0,1; z)G(l - z; y) 


3(y + z)3 

-8H(0,1,1; z) - 8H(1; z)G(l - z, 0; y) - 8H(1; z)G(0,1 - z; y) - 24H(1; z)G(0; y) + 24H(1,0; z) 
+8H(1, 0; z)G(l - z; y) - 8H(1, 0,1; z) + 8H(1,1; z)G(0; y) - 8H(1,1,0; z) + 8G(1 - z, 1 - z, 0; y) 
+8G(1 - z, 0,1 - z; y) + 24G(1 - z, 0; y) + 8G(0,1 - z, 1 - z; y) + 24G(0,1 - z; y)] 

■ [ - 24H(0; z) + 4H(0; z)G(l - z, 1 - z; y) + 16H(0; z)G(l - z; y) - 16H(0,1; z) 


3(y + z)2 

-4H(0,1; z)G(l - z; y) + 4H(0,1,1; z) + 4H(1; z)G(l - z, 0; y) + 4H(1; z)G(0,1 - z; y) 
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+16H(1; z)G(0; y) - 16H(1,0; z) - 4H(1, 0; 2 )G (1 -z;y) + 4H(1, 0,1; z) - 4H(1,1; z)G(0; y) 

+4H(1,1,0; z) - 4G(1 - z,l- z,0;y) - 4G(1 - z,0,1- z;y) - 16G(1 - 2 ,0; y) - 4G(0, l-z,l-z;y) 
-16G(0, l-z-,y) + 24G(0; y)] + [8H(0; z) - 8H(0; 2 )G (1 -z,l-z;y) + 8H(0,1; z)G(l - z; y) 

-8H(0,1,1; z) - 8H(1; z)G(l - z, 0; y) - 8H(1; z)G(0, l-z;y) + 8H(1, 0; z)G(l - z;y) - 8H(1, 0,1; z) 
+8H(1,1; z)G(0; y) - 8H(1,1,0; z) + 8G(1 - z, 1 - z, 0; y) + 8G(1 - z, 0,1 - z; y) + 8G(0,1 - z, 1 - z; y) 

-8G(0; y)] + ^ ^^3 [ 7 r 2 H( 0 ; z) + 7 r^H(l; z) - 7 r 2 G(l - z; y) - I 8 C 3 + 6H(0; z)G(0,1 - z; y) 

-6H(0, 0,1; z) + 6H(0,1; z)G(l - z; y) + 6H(0,1, 0; z) + 6H(1; z)G(l - z, -z; y) + 6H(1; z)G(0, -z; y) 

-6H(1, 0; z)G(l - z; y) - 6H(1, 0,1; z) + 6H(1,1, 0; z) - 6G(1 - z, -z, 1 - z; y) - 6G(0, -z, 1 - z; y)] 

2 2 
+ yY [6H(0; z)G(l - z; y) + 6H(1; z)G(-z; y) - 6 G(-z, 1 - z; y)] + [3H(0; z)G(l - z; y) 

1 27r^ 

+3H(1; z)G(-z; y) - 3G(-z, 1 - z; y)] + —-—-- [-^H(l; z) - 4H(0; z)G(-z, 1 - z; y) 

3(1 - y - z)(l - y)3 

+4H(0; z)G(0,1 - z; y) + 8H(0, 0,1; z) + 8H(0,1; z)G(-z; y) - 4H(0,1; z)G(0; y) + 4H(0,1,0; z) 

+4H(1; z)G(-z, -z; y) - 4H(1; z)G(-z, 0; y) + 4H(1,0; z)G(-z; y) - 4H(1,0; z)G(0; y) - 4H(1,0,1; z) 
-8H(1,1,0; z) + 4G(-z, 1 - z, 0; y) - 4G(-z, -z, 1 - z; y) + 4G(-z, 0,1 - z; y) - 4G(0,1, 0; y)] 

1 27r^ TT^ TT^ 

+ [— + yH(1; z) - yG( 1; y) - 3H(0; z)G(0,1 - z; y) + 4H(0; z)G(0; y) 

-3H(0; z)G(l, 0; y) - 6H(0, 0,1; z) - 6H(0,1; z)G(-z; y) + 3H(0,1; z)G(0; y) - 3H(0,1,0; z) 

-6H(1; z)G(-z, -z; y) - 3H(1; z)G(l, 0; y) + 4H(1,0; z) + 3H(1, 0; z)G(0; y) + 3H(1,0,1; z) 

+6H(1,1, 0; z) + 6 G(-z, -z, 1 - z; y) + 3G(0,1,0; y) + 3G(1,1 - z, 0; y) + 3G(1,0,1 - z; y) 

1 TT^ 27r^ 27r^ 

—4G(1, 0; y) + 3G(1,1, 0; y)] + y [-^ - ~ ~ 

+4H(0; z)G(0,1 - z; y) - H(0; z)G(0; y) - 4H(0,0,1; z) + H(0,1; z) + 4H(0,1; z)G(l - z; y) 

+4H(0,1,0; z) + 4H(1; z)G(l - z, -z; y) + H(l; z)G(-z; y) + 4H(1; z)G(0, -z; y) - H(l; z)G(0; y) 

-H(l, 0; z) - 4H(1,0; z)G(l - z; y) - 4H(1,0,1; z) + 4H(1,1,0; z) - 4G(1 - z, -z, 1 - z; y) 

+G(1 - z, 0; y) - G(-z, 1 - z; y) + G(0,1 - z; y) - 4G(0, -z, 1 - z; y) + 3G(0; y) + 2G(1,0; y)] 

+ 3 (y + ^)2 [ “ z) + 7 r 2 G(l - z; y) - 6H(0; z)G(l - z, 1 - z; y) - 6H(0; z)G(l - z; y) 

+6H(0; z)G(l - z, 0; y) - 6H(0; z)G(-z, 1 - z; y) + 6H(0; z)G(0,1 - z; y) - 6H(0,1; z) 

-6H(0,1; z)G(l - z; y) - 6H(0,1; z)G(-z; y) + 6H(0,1,0; z) + 6H(0,1,1; z) + 24H(1; z) 

-12H(1; z)G(l - z, -z; y) - 8H(1; z)G(l - z; y) + 6H(1; z)G(l - z, 0; y) - 12H(1; z)G(-z, 1 - z; y) 
-12H(1; z)G(-z, -z; y) - 12H(1; z)G(-z; y) + 6H(1; z)G(-z, 0; y) + 6H(1; z)G(0,1 - z; y) 

+6H(1; z)G(0, -z; y) + 6H(1; z)G(0; y) + 6H(1, 0; z) + 6H(1,0; z)G(l - z; y) + 6H(1,0; z)G(-z; y) 
-6H(1, 0; z)G(0; y) + 6H(1, 0,1; z) + 8H(1,1; z) + 12H(1,1; z)G(-z; y) - 6H(1,1; z)G(0; y) 

-6H(1,1,0; z) + 8G(1 - z, 1 - z; y) - 6G(1 - z, 1 - z, 0; y) + 12G(1 - z, -z, 1 - z; y) - 24G(1 - z; y) 
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—6G(1 — z, 0,1 — z; y) — 6G(1 — z, 0; y) + 12G(—z, 1 — z, 1 — z; y) + 12G(—z, 1 — z; y) 

-6G(-z, 1 - z, 0; y) + 12G(-z, -z, 1 - z; y) - 6G(-z, 0,1 - z; y) - 6G(0,1 - z, 1 - z; y) 

-6G(0,1 - z; y) - 6G(0, -z, 1 - z; y)] + [tt^ + yH(1; z) - yG(1 - z; y) - 12H(0; z) 

-12H(0; z)G(l - z; y) - 3H(0; z)G(l - z, 0; y) + 6H(0; z)G(0; y) - 18H(0,1; z) - 3H(0,1,0; z) 

-32H(1; z) - 8H(1; z)G(l - z; y) - 30H(1; z)G(-z; y) + 12H(1; z)G(0; y) + 12H(1, 0; z) 

-3H(1, 0; z)G(l - z; y) + 3H(1, 0; z)G(0; y) + 8H(1,1; z) + 3H(1,1, 0; z) + 8G(1 - z, 1 - z; y) 

+32G(1 - z; y) - 12G(1 - z, 0; y) + 3G(1 - z, 1, 0; y) + 30G(-z, 1 - z; y) - 12G(0,1 - z; y) 

T 27r^ 27r^ 

— 12G(0; y) + 3G(0,1, 0; y)] + "^ [-~ 2/) ”1—^)G(1 — z, 1 — z; y) 

-4H(0; z)G(l - z, 0; y) + 4H(0; z)G(l, 0; y) - 4H(0,1,1; z) - 9H(1; z) + 4H(1; z)G(l - z, -z; y) 

-4H(1; z)G(l - z, 0; y) + 8H(1; z)G(-z, 1 - z; y) - 4H(1; z)G(0,1 - z; y) + 4H(1; z)G(l, 0; y) 

-4H(1, 0; z)G(l - z; y) - 4H(1, 0,1; z) - 8H(1,1; z)G(-z; y) + 4H(1,1; z)G(0; y) - 4H(1,1,0; z) 

+4G(1 — z, 1 — z, 0; y) — 4G(1 — z, —z, 1 — z; y) + 9G(1 — z; y) + 4G(1 — z, 0,1 — z; y) 

-8G(-z, 1 - z, 1 - z; y) + 4G(0,1 - z, 1 - z; y) - 4G(1,1 - z, 0; y) - 4G(1,0,1 - z; y) - 4G(1,1,0; y)] 

TT^ 

+— [25 - 12H(0; z) + 6H(0; z)G(l - z; y) - 6H(0; z)G(l; y) - 23H(1; z) + 12H(1; z)G(l - z; y) 

-6H(1; z)G(l; y) - 6H(1, 0; z) - 6H(1,1; z) - 12G(1 - z, 1 - z; y) + 22G(1 - z; y) + 6G(1 - z, 0; y) 
-12G(0; y) + 6G(1,1 - z; y) + G(l; y)] + ^ [72C3 + 18C3H(1; z) - 36C3G(1 - z; y) + 18C3G(1; y) 

-5H(0; z) + 11H(0; z)G(l - z, 1 - z; y) + 5H(0; z)G(l - z; y) + 6H(0; z)G(l - z, 0,1 - z; y) 

+4H(0; 2)G(1 — z, 0; y) — 6H(0; z)G(—z, 1 — z, 1 — z; y) + 9H(0; z)G(—z, 1 — z; y) 

-6H(0; z)G(-z, -z, 1 - z; y) + 6H(0; z)G(-z, 0,1 - z; y) + 6H(0; z)G(0,1 - z, 1 - z; y) 

-22H(0; z)G(0,1 - z; y) + 6H(0; z)G(0, -z, 1 - z; y) + H(0; z)G(0; y) - 6H(0; z)G(0,0,1 - z; y) 
-6H(0; z)G(l, 0,1 - z; y) - 2H(0; z)G(l, 0; y) - 2H(0,0,1; z) + 6H(0, 0,1; z)G(l; y) - 11H(0,1; z) 
+12H(0,1; z)G(l - z, 1 - z; y) - 9H(0,1; z)G(l - z; y) - 6H(0,1; z)G(l - z, 0; y) 

-6H(0,1; z)G(-z, 1 - z; y) - 6H(0,1; z)G(-z, -z; y) + 7H(0,1; z)G(-z; y) + 13H(0,1; z)G(0; y) 
-6H(0,1; z)G(l, 1 - z; y) - 10H(0,1,0; z) + 6H(0,1, 0; z)G(l - z; y) - 6H(0,1,0; z)G(l; y) 

-11H(0,1,1; z) + 6H(0,1,1; z)G(-z; y) - 10H(1; z) + 12H(1; z)G(l - z, 1 - z, -z; y) 

+2H(1; z)G(l - z, -z; y) + 16H(1; z)G(l - z; y) - z)G(l - z, 0; y) 

— 12H(1; z)G(—z, 1 — z, —z; y) + 22H(1; z)G(—z, 1 — z; y) + 6H(1; z)G(—z, 1 — z, 0; y) 

-12H(1; z)G(-z, -z, 1 - z; y) - 12H(1; z)G(-z, -z, -z; y) + 16H(1; z)G(-z, -z; y) 

+6H(1; z)G(-z, -z, 0; y) - 6H(1; z)G(-z; y) + 6H(1; z)G(-z, 0,1 - z; y) + 6H(1; z)G(-z, 0, -z; y) 
-9H(1; z)G(-z, 0; y) + 6H(1; z)G(0,1 - z, -z; y) - 11H(1; z)G(0,1 - z; y) + 6H(1; z)G(0, -z, 1 - z; y) 
+6H(1; z)G(0, -z, -z; y) - 9H(1; z)G(0, -z; y) - 5H(1; z)G(0; y) - 6H(1; z)G(0, 0, -z; y) 

-6H(1; z)G(l, 1 - z, -z; y) - 6H(1; z)G(l, 0, -z; y) - 2H(1; z)G(l, 0; y) + H(l, 0; z) 
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-6H(1, 0; z)G(l -z,l-z]y) + 13H(1, 0; z)G(l -z]y) + 6H(1, 0; z)G{-z, -Z] y) - 9H(1, 0; z)G{-z] y) 
-6H(1, 0; z)G(0, -z; y) - 2H(1, 0; z)G(0; y) + 6H(1, 0; z)G(l, 1 - z; y) + 6H(1,0, 0,1; z) + 0,1; z) 

-12H(1, 0,1; z)G(l -z-y) + 6H(1, 0,1; z)G(-z; y) + 6H(1, 0,1; z)G(l; y) - 6H(1,0,1,0; z) 

-16H(1,1; z) + 12H(1,1; z)G(-z, -z; y) - 22H(1,1; z)G(-z; y) - 6H(1,1; z)G(-z, 0; y) 

-6H(1,1; z)G(0, -z; y) + 11H(1,1; z)G(0; y) - 12H(1,1,0; z) + 6H(1,1, 0; z)G(l - z; y) 

-6H(1,1,0; z)G(l; y) + 6H(1,1,0,1; z) - 6H(1,1,1,0; z) - 12G(1 - z, 1 - z, -z, 1 - z; y) 

— 16G(1 — z, 1 — z; y) + 11G(1 — z, 1 — z, 0; y) + 6G(1 — z, 1 — z, 1, 0; y) — 2G(1 — z, —z, 1 — z; y) 
+10G(1 - z; y) + 11G(1 - z, 0,1 - z; y) + 5G(1 - z, 0; y) - 6G(1 - z, 0,1, 0; y) - 24G(1 - z, 1, 0; y) 

—22G(—z, 1 — z, 1 — z; y) — 6G(—z, 1 — z, 1 — z, 0; y) + 12G(—z, 1 — z, —z, 1 — z; y) 

+6G(—z, 1 — z; y) — 6G(—z, 1 — z, 0,1 — z; y) + 9G(—z, 1 — z, 0; y) + 6G(—z, 1 — z, 1,0; y) 


+12G(—z, —z, 1 — z, 1 — z; y) — 16G(—z, 


-z, 1 - z; y) - 6G(-z, -z, 1 - z, 0; y) 


+12G(—z, —z, —z, 1 — z; y) — 6G(—z, —z, 0,1 — z; y) — 6G(—z, 0,1 — z, 1 — z; y) 

+9G(—z, 0,1 — z; y) — 6G(—z, 0, —z, 1 — z; y) + 6G(—z, 0,1,0; y) + 11G(0,1 — z, 1 — z; y) 

—6G(0,1 — z, —z, 1 — z; y) + 5G(0,1 — z; y) — 6G(0, —z, 1 — z, 1 — z; y) + 9G(0, —z, 1 — z; y) 

-6G(0, -z, -z, 1 - z; y) - 5G(0; y) + 6G(0, 0, -z, 1 - z; y) + G(0,1,0; y) + 6G(1,1 - z, -z, 1 - z; y) 
+2G(1,1 - z, 0; y) + 2G(1,0,1 - z; y) + 6G(1,0, -z, 1 - z; y) - 6G(1,0; y) - G(l, 1,0; y)] . (4.24) 

4.2.2 One-loop contribution to 

The finite remainder of the self-interference of the one-loop amplitude is decomposed as 


with 


y, z) 


V 


N'^ {Aii{y,z) + Aii{z,y)) + {Bii{y, z) + Bii{z,y)) 


+ ^ iCii{y,z) + Cii{z,y)) + NNp {Dii{y, z) + Dii{z,y)) 
+ ^ {Eiiiy, z) + Eii{z,y)) + Np {Fii{y, z) + Fii{z,y)) 


(4.25) 


Aiiiy,z) = 


[7r2G(0; y) + 10H(0; z)G(0; y) + 12H(0; z)G(0,0; y) + 6H(1, 0; z)G(0; y) + 21G(0; y) 


6(1-y)2 


■ V -h 37r2G(0; y) + 21 + 10H(0; z) 


+23G(0, 0; y) - 6G(0,1,0; y) - 12G(1, 0, 0; y)] + —_ 

-h36H(0; z)G(0; y) -f 36H(0; z)G(0, 0; y) -f 6H(1,0; z) -h 18H(1, 0; z)G(0; y) -h 73G(0; y) -f 33G(0, 0; y) 


■;-27r2G(0;y)-20H(0;z)G(0;y) 


-18G(0,1,0; y) - 6G(1,0; y) - 36G(1,0,0; y)] + ^ 

-24H(0; z)G(0,0; y) - 12H(1,0; z)G(0; y) - 42G(0; y) - 28G(0,0; y) + 12G(0,1,0; y) + 24G(1,0, 0; y)] 
Ttt^ 

+— [tt^ -P 169 -h 20H(0; z) -h 12H(0; z)G(0; y) -h 12H(1, 0; z) -h 20G(0; y) - 12G(1, 0; y)] 
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+ - [ + 72 + 60H(0; z) + 61H(0; z)G(0; y) + 30H(0; z)G(0, 0; y) - 9H(0; 2)G(0,1, 0; y) 

9 

-15H(0; z)G(l, 0; y) - 18H(0; 0)G(1,0, 0; y) + 25H(0,0; z) + 30H(0,0; z)G(0; y) + 18H(0,0; z)G(0,0; y) 
+15H(0,1, 0; z) + 9H(0,1, 0; z)G(0; y) + 36H(1, 0; z) + 15H(1,0; z)G(0; y) - 9H(1, 0; 2)G(1, 0; y) 

+30H(1,0, 0; z) + 18H(1,0, 0; z)G(0; y) + 9H(1,0,1, 0; z) + 18H(1,1, 0,0; z) + 60G(0; y) 

+25G(0,0; y) - 15G(0,1, 0; y) - 36G(1, 0; y) - 30G(1,0, 0; y) + 9G(1, 0,1, 0; y) + 18G(1,1,0, 0; y)] (4.26) 
Bii{y,z) = 

+ 7 ^ [ - 27r^H(0; 2)G(1 - z;y) - 27r^H(l; z)G{-z; y) + 27r^G(-z, 1 - z; y) + 3 - 42H(0; z)G(l - z; y) 

6y 

-20H(0; z)G(l - z, 0; y) + 12H(0; z)G(l - z, 1, 0; y) + 2H(0; z)G(-z, 1 - z; y) 

+12H(0; z)G(-z, 1 - z, 0; y) + 12H(0; z)G(-z, 0,1 - z; y) - 20H(0; z)G(0,1 - z; y) 

+12H(0; z)G(0, -z, 1 - z; y) + 12H(0; z)G(l, 1 - z, 0; y) + 12H(0; z)G(l, 0,1 - z; y) 

-4H(0, 0; z)G(l - z; y) - 24H(0, 0; z)G(l - z, 0; y) - 24H(0, 0; z)G(0,1 - z; y) - 2H(0,1; z)G(-z; y) 
-12H(0,1; z)G(-z, 0; y) - 12H(0,1; z)G(0, -z; y) - 12H(0,1, 0; z)G(l - z; y) - 42H(1; z)G(-z; y) 
-20H(1; z)G(-z, 0; y) + 12H(1; z)G(-z, 1, 0; y) - 20H(1; z)G(0, -z; y) + 12H(1; z)G(l, -z, 0; y) 

+12H(1; z)G(l, 0, -z; y) + 12H(1, 0; z)G(-z, 1 - z; y) - 2H(1, 0; z)G(-z; y) - 12H(1,0; z)G(-z, 0; y) 
-12H(1,0; z)G(0, -z; y) - 24H(1, 0,0; z)G(l - z; y) - 12H(1, 0,1; z)G(-z; y) - 24H(1,1,0; z)G(-z; y) 
+42G(-z, 1 - z; y) + 20G(-z, 1 - z, 0; y) - 12G(-z, 1 - z, 1, 0; y) + 20G(-z, 0,1 - z; y) 

-12G(-z, 1,1-z, 0; y) - 12G(-z, 1,0,1 - z; y) + 20G(0, -z, 1 - z; y) - 12G(1, -z, 1 - z, 0; y) 

-12G(1, -z, 0,1 - z; y) - 12G(1, 0, -z, 1 - z; y)] + ^ [ - z)G(l - z; y) 

6y 

+47r2H(l; z)G(-z; y) - Att‘^G{-z, 1 - z; y) - 10H(0; z) + 90H(0; z)G(l - z; y) + 40H(0; z)G(l - z, 0; y) 
-24H(0; z)G(l - z, 1, 0; y) - 16H(0; z)G(-z, 1 - z; y) - 24H(0; z)G(-z, 1 - z, 0; y) 

-24H(0; z)G(-z, 0,1 - z; y) + 40H(0; z)G(0,1 - z; y) - 24H(0; z)G(0, -z, 1 - z; y) - 6H(0; z)G(0; y) 
-24H(0; z)G(l, 1 - z, 0; y) - 24H(0; z)G(l, 0,1 - z; y) + 32H(0, 0; z)G(l - z; y) 

+48H(0, 0; z)G(l - z, 0; y) + 48H(0,0; z)G(0,1 - z; y) + 6H(0,1; z) + 16H(0,1; z)G(-z; y) 

+24H(0,1; z)G(-z, 0; y) + 24H(0,1; z)G(0, -z; y) + 24H(0,1,0; z)G(l - z; y) - 9H(1; z) 

+96H(1; z)G(-z; y) + 40H(1; z)G(-z, 0; y) - 24H(1; z)G(-z, 1, 0; y) + 40H(1; z)G(0, -z; y) 

-6H(1; z)G(0; y) - 24H(1; z)G(l, -z, 0; y) - 24H(1; z)G(l, 0, -z; y) - 6H(1,0; z) 

-24H(1, 0; z)G(-z, 1 - z; y) + 16H(1,0; z)G(-z; y) + 24H(1, 0; z)G(-z, 0; y) + 24H(1, 0; z)G(0, -z; y) 
+48H(1, 0,0; z)G(l - z; y) + 24H(1,0,1; z)G(-z; y) + 48H(1,1,0; z)G(-z; y) + 9G(1 - z; y) 

+6G(1 - z, 0; y) - 96G(-z, 1 - z; y) - 40G(-z, 1 - z, 0; y) + 24G(-z, 1 - z, 1,0; y) 

-40G(-z, 0,1 - z; y) + 24G(-z, 1,1 - z, 0; y) + 24G(-z, 1, 0,1 - z; y) + 6G(0,1 - z; y) 

-40G(0, -z, 1 - z; y) - 10G(0; y) + 24G(1, -z, 1 - z, 0; y) + 24G(1, -z, 0,1 - z; y) 

+24G(1,0, -z, 1 - z; y)] + ^(^4^ [ " " ^^(O; ^)G(0,1 - z; y) 

-10H(0; z)G(0; y) - 12H(0; z)G(0, 0; y) - 6H(0,1; z)G(0; y) - 12H(1; z)G(-z, 0; y) 
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-12H(1; z)G{0, -z; y) + 9H(1; 0 )G(O; y) + 12H(1; 2)G(0, 0; y) - 6H(1,0; z)G(0; y) - 9G(1 - 2 ,0; y) 
-12G(1 - z, 0,0; y) + 12G(-2,1 - z, 0; y) + 12G(-z, 0,1 - z; y) - 9G(0,1 - z; y) - 12G(0,1 - z, 0; y) 
+12G(0, -z, 1 - z; y) - 42G(0; y) - 12G(0,0,1 - z; y) - 26G(0,0; y) + 12G(0,1,0; y) + 24G(1,0, 0; y)] 

+ - 3^2G(0; y) - 42 - 10H(0; z) - 6H(0; z)G(l - z; y) - 18H(0; z)G(l - z, 0; y) 

-18H(0; z)G(0,1 - z; y) - 36H(0; z)G(0; y) - 36H(0; z)G(0, 0; y) - 6H(0,1; z) - 18H(0,1; z)G(0; y) 
+9H(1; z) - 12H(1; z)G(-z; y) - 36H(1; z)G(-z, 0; y) - 36H(1; z)G(0, -z; y) + 33H(1; z)G(0; y) 
+36H(1; z)G(0,0; y) - 6H(1,0; z) - 18H(1, 0; z)G(0; y) - 9G(1 - z; y) - 33G(1 - z, 0; y) 

-36G(1 - z, 0,0; y) + 12G(-z, 1 - z; y) + 36G(-z, 1 - z, 0; y) + 36G(-z, 0,1 - z; y) - 33G(0,1 - z; y) 
-36G(0,1 - z, 0; y) + 36G(0, -z, 1 - z; y) - 136G(0; y) - 36G(0,0,1 - z; y) - 6G(0,0; y) 


+36G(0,1, 0; y) + 12G(1,0; y) + 72G(1, 0,0; y)] + 

1 


■;-H(l;z)G(0;y) + G(l-z,0;y) 


+G(0,1 — z; y)] + 


(l-y)2(y + z) 


(l-y)2(y+ z)2 

[2H(1; z)G(0; y) - 2G(1 - z, 0; y) - 2G(0,1 - z; y) + G(0; y)] 


[ ~ H(l; z)G(0; y) + G(1 — z, 0; y) + G(0,1 — z; y) — G(0; y)] + —- - ———^ [ — H(l; z) 


(1-y) 

-2H(1; z)G(0; y) + G(1 - z; y) + 2G(1 - z, 0; y) + 2G(0,1 - z; y)] + 


{^-y){y + zf 


(l-y)(y + z) 


[l + 2H(l;z) 

■[7r2G(0;y) 


+2H(1; z)G(0; y) - 2G(1 - z; y) - 2G(1 - z, 0; y) - 2G(0,1 - z; y) + 2G(0; y)] + —_ 

-3 + 12H(0; z)G(l - z, 0; y) + 12H(0; z)G(0,1 - z; y) + 10H(0; z)G(0; y) + 12H(0; z)G(0,0; y) 

+6H(0,1; z)G(0; y) - 3H(1; z) + 18H(1; z)G(-z, 0; y) + 18H(1; z)G(0, -z; y) - 18H(1; z)G(0; y) 
-12H(1; z)G(0, 0; y) + 6H(1, 0; z)G(0; y) + 3G(1 - z; y) + 18G(1 - z, 0; y) + 12G(1 - z, 0, 0; y) 
-18G(-z, 1 - z, 0; y) - 18G(-z, 0,1 - z; y) + 18G(0,1 - z; y) + 12G(0,1 - z, 0; y) 

-18G(0, -z, 1 - z; y) + 66G(0; y) + 12G(0,0,1 - z; y) + 14G(0,0; y) - 12G(0,1,0; y) - 24G(1,0, 0; y)] 


+ ^ - G y) + 11H(0; z)G(l - z; y) + 12H(0; z)G(l - z, 0; y) 

+12H(0; z)G(0,1 - z; y) - 11H(0,1; z) - 12H(0,1; z)G(0; y) - 42H(1; z) - 11H(1; z)G(0; y) 
+12H(1; z)G(l, 0; y) - 0; z) + 12H(1, 0; z)G(l - z; y) - 12H(1, 0; z)G(0; y) - 12H(1, 0,1; z) 

-24H(1,1,0; z) + 42G(1 - z; y) + 11G(1 - z, 0; y) - 12G(1 - z, 1, 0; y) + 11G(0,1 - z; y) 

-12G(1,1 - z, 0; y) - 12G(1,0,1 - z; y)] + [2^' + 42 + 11H(0; z) + 12H(0; z)G(0; y) 

+12H(1, 0; z) + 11G(0; y) - 12G(1,0; y)] + — [ - 8 - 2H(0; z)G(l - z; y) - 2H(0,1; z) + 3H(1; z) 
-4H(1; z)G(-z; y) + 2H(1; z)G(0; y) - 3G(1 - z; y) - 2G(1 - z, 0; y) + 4G(-z, 1 - z; y) 


-2G(0,1 - z; y) + 2G(1, 0; y)l + - [ - 96 - 40H(0; z) - 39H(0; z)G(l - z; y) - 29H(0; z)G(l - z, 0; y) 

■' 6 

-12H(0; z)G(l - z, 0, 0; y) + 6H(0; z)G(l - z, 1,0; y) + 20H(0; z)G(-z, 1 - z; y) 

+12H(0; z)G(-z, 1 - z, 0; y) + 12H(0; z)G(-z, 0,1 - z; y) - 29H(0; z)G(0,1 - z; y) 
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-12H(0; z)G(0,1 - 0 , 0; y) + 12H(0; z)G(0, -z, 1 - z; y) - 24H(0; z)G(0; y) - 12H(0; z)G(0, 0,1 - z; y) 
+6H(0; z)G(0,1,0; y) + 6H(0; z)G(l, 1 - z, 0; y) + 6H(0; z)G(l, 0,1 - z; y) + 10H(0; z)G(l, 0; y) 
+12H(0; z)G(l, 0,0; y) - 20H(0,0; z)G(l -z;y)- 12H(0,0; z)G(l - z, 0; y) - 12H(0, 0; z)G(0,1 - z; y) 
-20H(0,0,1; z) - 12H(0,0,1; z)G(0; y) - 9H(0,1; z) - 20H(0,1; z)G(-z; y) - 12H(0,1; z)G(-z, 0; y) 
-12H(0,1; z)G(0, -z; y) + 9H(0,1; z)G(0; y) + 12H(0,1; z)G(0, 0; y) + 6H(0,1; z)G(l, 0; y) 

-10H(0,1, 0; z) - 6H(0,1, 0; z)G(l - z;y) - 6H(0,1,0; z)G(0; y) - 6H(0,1,0,1; z) - 12H(0,1,1,0; z) 
+36H(1; z) - 48H(1; z)G(-z; y) - 20H(1; z)G(-z, 0; y) + 12H(1; z)G(-z, 1, 0; y) - 20H(1; z)G(0, -z; y) 
+39H(1; z)G(0; y) + 20H(1; z)G(0, 0; y) - 6H(1; z)G(0,1,0; y) + 12H(1; z)G(l, -z, 0; y) 
+12H(l;z)G(l,0,-z;y)-9H(l;z)G(l,0;y)-12H(l;z)G(l,0,0;y)-9H(l,0;z) 

-9H(1,0; z)G(l - z; y) - 6H(1, 0; z)G(l - z, 0; y) + 12H(1, 0; z)G(-z, 1 - z; y) - 20H(1, 0; z)G(-z; y) 
-12H(1,0; z)G(-z, 0; y) - 6H(1, 0; z)G(0,1 - z; y) - 12H(1, 0; z)G(0, -z; y) + 19H(1,0; z)G(0; y) 
+12H(1,0; z)G(0,0; y) + 6H(1,0; z)G(l, 0; y) - 12H(1, 0,0; z)G(l - z; y) - 12H(1,0,0,1; z) 

+9H(1,0,1; z) - 12H(1,0,1; z)G(-z; y) + 6H(1,0,1; z)G(0; y) - 6H(1,0,1, 0; z) + 18H(1,1,0; z) 
-24H(1,1, 0; z)G(-z; y) + 12H(1,1,0; z)G(0; y) - 36G(1 - z; y) - 39G(1 - z, 0; y) - 20G(1 - z, 0, 0; y) 
+6G(1 - z, 0,1,0; y) + 9G(1 - z, 1,0; y) + 12G(1 - z, 1, 0,0; y) + 48G(-z, 1 - z; y) 

+20G(-z, 1 - z, 0; y) - 12G(-z, 1 - z, 1, 0; y) + 20G(-z, 0,1 - z; y) - 12G(-z, 1,1-z, 0; y) 

-12G(-z, 1,0,1 - z; y) - 39G(0,1 - z; y) - 20G(0,1 - z, 0; y) + 6G(0,1 - z, 1,0; y) 

+20G(0, -z, 1 - z; y) - 40G(0; y) - 20G(0,0,1 - z; y) + 6G(0, 1,1-z, 0; y) + 6G(0,1, 0,1 - z; y) 
+10G(0,1,0; y) + 9G(1,1 - z, 0; y) + 12G(1,1 - z, 0,0; y) - 12G(1, -z, 1 - z, 0; y) 

-12G(1, -z, 0,1 - z; y) + 9G(1, 0,1 - z; y) + 12G(1, 0,1 - z, 0; y) - 12G(1,0, -z, 1 - z; y) 

+48G(1,0; y) + 12G(1,0, 0,1 - z; y) + 20G(1, 0, 0; y) - 12G(1,0,1,0; y) - 24G(1,1, 0, 0; y)] 

[H(0; z) - 2H(0; z)G(l - z; y) - 2H(0,1; z) + 2H(1; z) - 4H(1; z)G(-z; y) + 2H(1; z)G(0; y) 

-2G(1 - z; y) - 2G(1 - z, 0; y) + 4G(-z, 1 - z; y) - 2G(0,1 - z; y) + G(0; y) + 2G(1, 0; y)] 

+ i [21H(0; z) - 53H(0; z)G(l - z; y) - 52H(0; z)G(l - z, 0; y) - 24H(0; z)G(l - z, 0, 0; y) 

+12H(0; z)G(l - z, 1,0; y) + 22H(0; z)G(-z, 1 - z; y) + 24H(0; z)G(-z, 1 - z, 0; y) 

+24H(0; z)G(-z, 0,1 - z; y) - 52H(0; z)G(0,1 - z; y) - 24H(0; z)G(0,1 - z, 0; y) 

+24H(0; z)G(0, -z, 1 - z; y) + 20H(0; z)G(0; y) - 24H(0; z)G(0,0,1 - z; y) + 12H(0; z)G(0,0; y) 
+12H(0; z)G(0,1,0; y) + 12H(0; z)G(l, 1 - z, 0; y) + 12H(0; z)G(l, 0,1 - z; y) + 14H(0; z)G(l, 0; y) 
+24H(0; z)G(l, 0,0; y) + 2H(0,0; z) - 4H(0, 0; z)G(l - z; y) - 24H(0, 0; z)G(l - z, 0; y) 

-24H(0,0; z)G(0,1 - z; y) + 12H(0, 0; z)G(0; y) - 40H(0,0,1; z) - 24H(0, 0,1; z)G(0; y) - 31H(0,1; z) 
-22H(0,1; z)G(-z; y) - 24H(0,1; z)G(-z, 0; y) - 24H(0,1; z)G(0, -z; y) + 30H(0,1; z)G(0; y) 

+24H(0,1; z)G(0,0; y) + 12H(0,1; z)G(l, 0; y) - 14H(0,1, 0; z) - 12H(0,1,0; z)G(l - z; y) 

-12H(0,1, 0; z)G(0; y) - 12H(0,1, 0,1; z) - 24H(0,1,1, 0; z) + 42H(1; z) - 84H(1; z)G(-z; y) 
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-22H(1; z)G{-z, 0; y) + 24H(1; z)G{-z, 1,0; y) - 22H(1; 2)G(0, - 0 ; y) + 53H(1; z)G(0; y) 

+4H(1; z)G(0,0; y) - 12H(1; z)G(0,1, 0; y) + 24H(1; z)G(l, -z, 0; y) + 24H(1; z)G(l, 0, - 0 ; y) 

-12H(1; z)G(l, 0; y) - 24H(1; z)G(l, 0, 0; y) + 11H(1, 0; z) - 12H(1, 0; z)G(l - z; y) 

-12H(1,0; z)G(l - z, 0; y) + 24H(1,0; z)G(-z, l-z;y)- 22H(1,0; z)G(-z; y) 

-24H(1,0; z)G(-z, 0; y) - 12H(1,0; z)G(0,1 - z; y) - 24H(1,0; z)G(0, -z; y) + 38H(1,0; z)G(0; y) 
+24H(1,0; z)G(0,0; y) + 12H(1,0; z)G(l, 0; y) + 12H(1,0, 0; z) - 24H(1, 0,0; z)G(l - z; y) 

-24H(1,0, 0,1; z) + 12H(1,0,1; z) - 24H(1, 0,1; z)G(-z; y) + 12H(1, 0,1; z)G(0; y) - 12H(1,0,1,0; z) 
+24H(1,1, 0; z) - 48H(1,1, 0; z)G(-z; y) + 24H(1,1, 0; z)G(0; y) - 42G(1 - z; j/) - 53G(1 - z, 0; y) 
-4G(1 - z, 0,0; y) + 12G(1 - z, 0,1,0; y) + 12G(1 - z, 1,0; y) + 24G(1 - z, 1, 0,0; y) 

+84G(-z, l-z;y) + 22G(-z, 1 - z, 0; y) - 24G(-z, 1 - z, 1, 0; y) + 22G(-z, 0,1 - z; y) 

-24G(-z, 1,1 - z, 0; y) - 24G(-z, 1, 0,1 - z; y) - 53G(0,1 - z; y) - 4G(0,1 - z, 0; y) 

+12G(0,1 - z, 1,0; y) + 22G(0, -z, 1 - z; y) + 21G(0; y) - 4G(0, 0,1 - z; y) + 2G(0,0; y) 

+12G(0,1,1 - z, 0; y) + 12G(0,1, 0,1 - z; y) - 4G(0,1, 0; y) + 12G(1,1 - z, 0; y) + 24G(1,1 - z, 0, 0; y) 
-24G(1, -z, 1 - z, 0; y) - 24G(1, -z, 0,1 - z; y) + 12G(1,0,1 - z; y) + 24G(1, 0,1 - z, 0; y) 

-24G(1,0, -z, 1 - z; y) + 42G(1,0; y) + 24G(1,0, 0,1 - z; y) - 8G(1,0, 0; y) - 24G(1, 0,1, 0; y) 
-48G(l,l,0,0;y)] , (4.27) 

Ciiiy,z) = 

3 ^^ [ - 2H(0; z)G(l - z, -z, 1 - z;y) - 4H(0; z)G(-z, 1 - z, 1 - z;y) 

yi 

+4H(0,0; z)G(l - z, 1 - z; y) + 2H(0,1; z)G(l - z, -z; y) + 2H(0,1; z)G(-z, 1 - z; y) 

-2H(1; z)G(-z, 1 - z, -z; y) - 4H(1; z)G(-z, -z, 1 - z; y) + 2H(1, 0; z)G(l - z, -z; y) 

+2H(1,0; z)G(-z, 1 - z; y) + 4H(1,1; z)G(-z, -z; y) + 2G(-z, 1 - z, -z, 1 - z; y) 

+4G(-z, -z, 1 - z, 1 - z; y)] + ^ [ - 8H(0; z)G(l - z, 1 - z; y) + 12H(0; z)G(l - z, -z, 1 - z; y) 

yZ 

+24H(0; z)G(-z, 1 - z, 1 - z; y) - 2H(0; z)G(-z, 1 - z; y) - 24H(0, 0; z)G(l - z, 1 - z; y) 

+4H(0, 0; z)G(l - z; y) - 12H(0,1; z)G(l - z, -z; y) + 4H(0,1; z)G(l - z; y) 

-12H(0,1; z)G(-z, 1 - z; y) + 2H(0,1; z)G(-z; y) - 4H(1; z)G(l - z, -z; y) 

+12H(1; z)G(—z, 1 — z, —z; y) — 8H(1; z)G(—z, 1 — z; y) + 24H(1; z)G(—z, —z, 1 — z; y) 

-12H(1, 0; z)G(l - z, -z; y) + 4H(1,0; z)G(l - z; y) - 12H(1,0; z)G(-z, 1 - z; y) 

+2H(1,0; z)G(-z; y) - 24H(1,1; z)G(-z, -z; y) + 8H(1,1; z)G(-z; y) + 4G(1 - z, -z, 1 - z; y) 

+8G(—z, 1 — z, 1 — z; y) — 12G(—z, 1 — z, —z, 1 — z; y) — 24G(—z, —z, 1 — z, 1 — z; y)] 

z2 

H —- [4H(0; z)G(l — z, 1 — z; y) — 8H(0; z)G(l — z, —z, 1 — z; y) — 16H(0; z)G(—z, 1 — z, 1 — z; y) 

+2H(0; z)G(-z, 1 - z; y) + 16H(0,0; z)G(l - z, 1 - z; y) - 4H(0,0; z)G(l - z; y) 

+8H(0,1; z)G(l - z, -z; y) - 2H(0,1; z)G(l - z; y) + 8H(0,1; z)G(-z, 1 - z; y) 

-2H(0,1; z)G(-z; y) + 2H(1; z)G(l - z, -z; y) - 8H(1; z)G(-z, 1 - z, -z; y) + 4H(1; z)G(-z, 1 - z; y) 
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-16H(1; z)G{-z, -z, l-z-,y) + 8H(1, 0; z)G(l - 2 , -z; y) - 2H(1, 0; z)G(l - z; y) 

+8H(1, 0; z)G(-z, l-z;y)- 2H(1,0; z)G(-z; y) + 16H(1,1; z)G(-z, -z; y) - 4H(1,1; z)G(-z; y) 
-2G(1 - z, -z, 1- z;y)- 4G(-z, 1 - z, 1 - z; y) + 8 G(-z, 1 - z, -z, 1 - z; y) 

+16G(—z, —z, 1 — z, 1 — z; y)] + ^ [ + 4H(0; z:)G(l — z, 1 — z; y) — 4H(0; z:)G(l — z, —z, 1 — z; y) 

-8H(0; z)G(-z, 1 - z, 1 - z; y) + 8H(0,0; z)G(l - z, 1 - z; y) + 4H(0,1; z)G(l - z, -z; y) 

-2H(0,1; z)G(l - z; y) + 4H(0,1; z)G(-z, 1 - z; y) + 2H(1; z)G(l - z, -z; y) 

-4H(1; z)G(-z, 1 - z, -z; y) + 4H(1; z)G(-z, 1 - z; y) - 8H(1; z)G(-z, -z, 1 - z; y) 

+4H(1,0; z)G(l - z, -z; y) - 2H(1,0; z)G(l - z; y) + 4H(1, 0; z)G(-z, 1 - z; y) 

+8H(1,1; z)G(-z, -z; y) - 4H(1,1; z)G(-z; y) - 2G(1 - z, -z, 1 - z; y) - 4G(-z, 1 - z, 1 - z; y) 

+4G(-z, 1 - z, -z, 1 - z; y) + 8 G(-z, -z, 1 - z, 1 - z; y)] + - [7 + 14H(0; z)G(l - z, 1 - z; y) 

V 4 

+4H(0; ^)G(1 — z, 1 — z, 0; y) — 14H(0; z:)G(l — z, —z, 1 — z; y) + 5H(0; z)G(l — z; y) 

+4H(0; z)G(l - z, 0,1 - z; y) - 2H(0; z)G(l - z, 1,0; y) - 28H(0; z)G(-z, 1 - z, 1 - z; y) 

+7H(0; z)G(-z, 1 - z; y) + 4H(0; z)G(0,1 - z, 1 - z; y) - 2H(0; z)G(l, 1 - z, 0; y) 

-2H(0; z)G(l, 0,1 - z; y) + 2H(0, 0; z) + 24H(0,0; z)G(l - z, 1 - z; y) - 14H(0,0; z)G(l - z; y) 

+4H(0, 0,1; z)G(l - z; y) + 2H(0,1; z) + 14H(0,1; z)G(l - z, -z; y) - 7H(0,1; z)G(l - z; y) 

-2H(0,1; z)G(l - z, 0; y) + 12H(0,1; z)G(-z, 1 - z; y) - 7H(0,1; z)G(-z; y) - 2H(0,1; z)G(0,1 - z; y) 
+2H(0,1,0; z)G(l - z; y) + 4H(0,1,1; z)G(-z; y) + 7H(1; z)G(l - z, -z; y) + 2H(1; z)G(l - z, -z, 0; y) 
-H(l;z)G(l - z;y) + 2H(l;z)G(l - z,0,-z;y) - 16H(1; z)G(-z, 1 - z,-z;y) 

+14H(1; z)G(—z, 1 — z; y) + 4H(1; z)G(—z, 1 — z, 0; y) — 32H(1; z)G(—z, —z, 1 — z; y) 

+7H(1; z)G(-z; y) + 4H(1; z)G(-z, 0,1 - z; y) - 2H(1; z)G(-z, 1,0; y) + 2H(1; z)G(0,1 - z, -z; y) 
+4H(1; z)G(0, -z, 1 - z; y) - 2H(1; z)G(l, -z, 0; y) - 2H(1; z)G(l, 0, -z; y) + 2H(1, 0; z) 

+14H(1, 0; z)G(l - z, -z; y) - 7H(1,0; z)G(l - z; y) - 2H(1,0; z)G(l - z, 0; y) 

+14H(1, 0; z)G(-z, 1 - z; y) - 7H(1,0; z)G(-z; y) - 2H(1,0; z)G(0,1 - z; y) 

+2H(1, 0,1; z)G(-z; y) + H(l, 1; z) + 32H(1,1; z)G(-z, -z; y) - 14H(1,1; z)G(-z; y) 

-4H(1,1; z)G(-z, 0; y) - 4H(1,1; z)G(0, -z; y) + G(1 - z, 1 - z; y) - 7G(1 - z, -z, 1 - z; y) 

-2G(1 - z, -z, 1 - z, 0; y) - 2G(1 - z, -z, 0,1 - z; y) - 2G(1 - z, 0, -z, 1 - z; y) 

— 14G(—z, 1 — z, 1 — z; y) — 4G(—z, 1 — z, 1 — z, 0; y) + 16G(—z, 1 — z, —z, 1 — z; y) 

—7G(—z, 1 — z; y) — 4G(—z, 1 — z, 0,1 — z; y) + 2G(—z, 1 — z, 1, 0; y) + 32G(—z, —z, 1 — z, 1 — z; y) 
-4G(-z, 0,1 - z, 1 - z; y) + 2G(-z, 1,1 - z, 0; y) + 2G(-z, 1,0,1 - z; y) - 2G(0,1 - z, -z, 1 - z; y) 
-4G(0, -z, 1 - z, 1 - z; y) + 2G(1, -z, 1 - z, 0; y) + 2G(1, -z, 0,1 - z; y) + 2G(1,0, -z, 1 - z; y)] 

+ - [ + 4 - 20H(0; z)G(l - z, 1 - z; y) - 8H(0; z)G(l - z, 1 - z, 0; y) + 16H(0; z)G(l - z, -z, 1 - z; y) 

y 

-13H(0; z)G(l -z-,y)- 8H(0; z)G(l - z, 0,1 - z; y) + 4H(0; z)G(l - z, 1,0; y) 

+32H(0; z)G(-z, 1 - z, 1 - z; y) - 4H(0; z)G(-z, 1 - z; y) - 8H(0; z)G(0, 1 - z, 1 - z; y) 
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+4H(0; z)G(l, 1 - z, 0; y) + 4H(0; z)G(l, 0,1 - z; y) - 24H(0,0; z)G(l - z, 1 - z; y) 

+8H(0, 0; z)G(l - z;y) - 8H(0, 0,1; z)G(l -z;y)- H(0,1; z) - 16H(0,1; z)G(l - z, -z; y) 

+10H(0,1; z)G(l - z; y) + 4H(0,1; z)G(l - z, 0; y) - 12H(0,1; z)G(-z, 1 - z; y) + 4H(0,1; z)G(-z; y) 

+4H(0,1; z)G(0,1 - z; y) - 4H(0,1,0; z)G(l - z; y) - 8H(0,1,1; z)G(-z; y) - ^H(l; z) 

-10H(1; z)G(l - z, -z; y) - 4H(1; z)G(l - z, -z, 0; y) + 2H(1; z)G(l - z; y) - 4H(1; z)G(l - z, 0, -z; y) 
+20H(1; z)G(-z, 1 - z, -z; y) - 20H(1; z)G(-z, 1 - z; y) - 8H(1; z)G(-z, 1 - z, 0; y) 

+40H(1; z)G(-z, -z, 1 - z; y) - 14H(1; z)G(-z; y) - 8H(1; z)G(-z, 0,1 - z; y) + 4H(1; z)G(-z, 1,0; y) 
-4H(1; z)G(0,1 - z, -z; y) - 8H(1; z)G(0, -z, 1 - z; y) - H(l; z)G(0; y) + 4H(1; z)G(l, -z, 0; y) 

+4H(1; z)G(l, 0, -z; y) - 2H(1,0; z) - 16H(1, 0; z)G(l - z, -z; y) + 10H(1, 0; z)G(l - z; y) 

+4H(1, 0; z)G(l - z, 0; y) - 16H(1, 0; z)G(-z, 1 - z; y) + 4H(1, 0; z)G(-z; y) + 4H(1, 0; z)G(0,1 - z; y) 
-4H(1, 0,1; z)G(-z; y) - 2H(1,1; z) - 40H(1,1; z)G(-z, -z; y) + 20H(1,1; z)G(-z; y) 

+8H(1,1; z)G(-z, 0; y) + 8H(1,1; z)G(0, -z; y) - 2G(1 - z, 1 - z; y) + 10G(1 - z, -z, 1 - z; y) 

+4G(1 - z, -z, 1 - z, 0; y) + 4G(1 - z, -z, 0,1 - z; y) + ^G(l - z; y) + 4G(1 - z, 0, -z, 1 - z; y) 

+G(1 — z, 0; y) + 20G(—z, 1 — z, 1 — z; y) + 8G(—z, 1 — z, 1 — z, 0; y) — 20G(—z, 1 — z, —z, 1 — z; y) 
+14G(—z, 1 — z; y) + 8G(—z, 1 — z, 0,1 — z; y) — 4G(—z, 1 — z, 1,0; y) — 40G(—z, —z, 1 — z, 1 — z; y) 
+8G(—z, 0,1 — z, 1 — z; y) — 4G(—z, 1,1 — z, 0; y) — 4G(—z, 1,0,1 — z; y) + 4G(0,1 — z, —z, 1 — z; y) 
+G(0,1 - z; y) + 8G(0, -z, 1 - z, 1 - z; y) - 4G(1, -z, 1 - z, 0; y) - 4G(1, -z, 0,1 - z; y) 

-4G(1, 0, -z, 1 - z; y) - G(l, 0; y)] + ^(^4^ [2H(0; z)G(l - z, 0; y) + 2H(0; z)G(0,1 - z; y) 

+2H(0,1; z)G(0; y) + 4H(1; z)G(-z, 0; y) + 4H(1; z)G(0, -z; y) - 3H(1; z)G(0; y) - 4H(1; z)G(0,0; y) 
+3G(1 - z, 0; y) + 4G(1 - z, 0, 0; y) - 4G(-z, 1 - z, 0; y) - 4G(-z, 0,1 - z; y) + 3G(0,1 - z; y) 

+4G(0,1 - z, 0; y) - 4G(0, -z, 1 - z; y) + 7G(0; y) + 4G(0,0,1 - z; y) + 1G(0, 0; y) - 2G(0,1, 0; y) 
-4G(1,0, 0; y)] + 2/) + 6H(0; z)G(l - z, 0; y) + 6H(0; z)G(0,1 - z; y) 

+2H(0,1; z) + 6H(0,1; z)G(0; y) - 3H(1; z) + 4H(1; z)G(-z; y) + 12H(1; z)G(-z, 0; y) 

+12H(1; z)G(0, -z; y) - z)G(0; y) - 12H(1; z)G(0,0; y) + 3G(1 - z; y) + 11G(1 - z, 0; y) 

+12G(1 - z, 0,0; y) - 4G(-z, 1 - z; y) - 12G(-z, 1 - z, 0; y) - 12G(-z, 0,1 - z; y) + 11G(0,1 - z; y) 
+12G(0,1 - z, 0; y) - 12G(0, -z, 1 - z; y) + 21G(0; y) + 12G(0,0,1 - z; y) - 9G(0, 0; y) - 6G(0,1,0; y) 

-2G(1,0; y) - 12G(1,0,0; y)] + ^ [H(1; z)G( 0; y) - G(1 - z, 0; y) - G(0,1 - z; y)] 

+ + t - 2H(1; z)G(0; y) + 2G(1 - z, 0; y) + 2G(0,1 - z; y) - G(0; y)] 

^)G(0; y) — G(1 — z, 0; y) — G(0,1 — z; y) + G(0; y)] 

^(l-y)(y + z )2 [ + + (l-y)(y + z) ^ ^ ^ 
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+ -^ [1 - 2H(0; 2 )G (1 - 0; y) - 2H(0; z)G(0,1 - z; y) + H(l; z) - 2H(1; z)G(-z, 0; y) 

-2H(1; z)G(0, -z; y) + 3H(1; z)G(0; y) - G(1 - z; y) - 3G(1 - z, 0; y) + 2G(-z, 1 - z, 0; y) 

+2G(-z, 0,1 - z; y) - 3G(0,1 - z; y) + 2G(0, -z, 1 - z; y) - 8G(0; y) - 2G(0,0; y)] 

+ 2 (ylz )2 [ - z)G(l - z, 1 - z; y) + H(0; z)G(l - z; y) - H(0,1; z) + 4H(0,1,1; z) + 14H(1; z) 

-4H(1; z)G(l - z, -z; y) + 8H(1; z)G(l - z; y) + 4H(1; z)G(l - z, 0; y) - 8H(1; z)G(-z, 1 - z; y) 
+4H(1; z)G(0,1 - z; y) - H(l; z)G(0; y) - 2H(1; z)G(l, 0; y) - H(l, 0; z) + 2H(1, 0; z)G(l - z; y) 

+2H(1, 0,1; z) - 8H(1,1; z) + 8H(1,1; z)G(-z; y) - 4H(1,1; z)G(0; y) - 8G(1 - z, 1 - z; y) 

-4G(1 - z, 1 - z, 0; y) + 4G(1 - z, -z, 1 - z; y) - 14G(1 - z; y) - 4G(1 - z, 0,1 - z; y) + G(1 - z, 0; y) 
+2G(1 - z, 1,0; y) + 8 G(-z, 1 - z, 1 - z; y) - 4G(0,1 - z, 1 - z; y) + G(0,1 - z; y) + 2G(1,1 - z, 0; y) 

+2G(1,0,1 - z; y)] + . [ - 14 + H(0; z) - 2H(0; z)G(l - z; y) - 2H(0,1; z) + 6H(1; z) 

2 (y + z)'- 

-4H(1; z)G(-z; y) + 2H(1; z)G(0; y) - 6G(1 - z; y) - 2G(1 - z, 0; y) + 4G(-z, 1 - z; y) 

-2G(0,1 - z; y) + G(0; y) + 2G(1, 0; y)] + T [8 + 3H(0; z)G(l - z, 1 - z; y) 

+2H(0; z)G(l - z, 1 - z, 0; y) - 2H(0; z)G(l - z, -z, 1 - z; y) + 4H(0; z)G(l - z; y) 

+2H(0; z)G(l - z, 0,1 - z; y) - H(0; z)G(l - z, 1, 0; y) - 4H(0; z)G(-z, 1 - z, 1 - z; y) 

+2H(0; z)G(0,1 - z, 1 - z; y) - H(0; z)G(l, 1 - z, 0; y) - H(0; z)G(l, 0,1 - z; y) 

+2H(0, 0; z)G(l - z, 1 - z; y) + 2H(0, 0,1; z)G(l - z; y) + 2H(0,0,1,1; z) + 4H(0,1; z) 

+2H(0,1; z)G(l - z, -z; y) - H(0,1; z)G(l, 0; y) + H(0,1, 0; z)G(l - z; y) + H(0,1,0,1; z) 

-3H(0,1,1; z) + 4H(0,1,1; z)G(-z; y) - 2H(0,1,1; z)G(0; y) - 6H(1; z) + 3H(1; z)G(l - z, -z; y) 
+2H(1; z)G(l - z, -z, 0; y) + 2H(1; z)G(l - z, 0, -z; y) - 3H(1; z)G(l - z, 0; y) 

-2H(1; z)G(l - z, 0,0; y) - 4H(1; z)G(-z, 1 - z, -z; y) + 6H(1; z)G(-z, 1 - z; y) 

+4H(1; z)G(-z, 1 - z, 0; y) - 8H(1; z)G(-z, -z, 1 - z; y) + 8H(1; z)G(-z; y) 

+4H(1; z)G(-z, 0,1 - z; y) - 2H(1; z)G(-z, 1,0; y) + 2H(1; z)G(0,1 - z, -z; y) - 3H(1; z)G(0,1 - z; y) 
-2H(1; z)G(0,1 - z, 0; y) + 4H(1; z)G(0, -z, 1 - z; y) - 4H(1; z)G(0; y) - 2H(1; z)G(0,0,1 - z; y) 
+H(1; z)G(0,1, 0; y) - 2H(1; z)G(l, -z, 0; y) - 2H(1; z)G(l, 0, -z; y) + 2H(1; z)G(l, 0, 0; y) 

+2H(1, 0; z)G(l - z, -z; y) - H(l, 0; z)G(l - z, 0; y) + 2H(1,0; z)G(-z, 1 - z; y) 

-H(l, 0; z)G(0,1 - z; y) + 2H(1, 0,1; z)G(-z; y) - H(l, 0,1; z)G(0; y) + 8H(1,1; z)G(-z, -z; y) 

-6H(1,1; z)G(-z; y) - 4H(1,1; z)G(-z, 0; y) - 4H(1,1; z)G(0, -z; y) + 3H(1,1; z)G(0; y) 

+2H(1,1; z)G(0, 0; y) + 3G(1 - z, 1 - z, 0; y) + 2G(1 - z, 1 - z, 0,0; y) - 3G(1 - z, -z, 1 - z; y) 

-2G(1 - z, -z, 1 - z, 0; y) - 2G(1 - z, -z, 0,1 - z; y) + 6G(1 - z; y) + 3G(1 - z, 0,1 - z; y) 

+2G(1 - z, 0,1 - z, 0; y) - 2G(1 - z, 0, -z, 1 - z; y) + 4G(1 - z, 0; y) + 2G(1 - z, 0, 0,1 - z; y) 

-G(l - z, 0,1,0; y) - 2G(1 - z, 1,0, 0; y) - 6 G(-z, 1 - z, 1 - z; y) - 4G(-z, 1 - z, 1 - z, 0; y) 

+4G(-z, 1 - z, -z, 1 - z; y) - 8G(-z, 1 - z; y) - 4G(-z, 1 - z, 0, 1 - z; y) + 2G(-z, 1 - z, 1, 0; y) 
+8G(—z, —z, 1 — z, 1 — z; y) — 4G(—z, 0, 1 — z, 1 — z; y) + 2G(—z, 1, 1 — z, 0; y) + 2G(—z, 1,0,1 — z; y) 


42 



+3G(0, 1- z,l- z;y) + 2G(0 ,1 - z,l - z,0;y) - 2G(0, l-z,-z,l- Z] y) + 4G(0, 1-z]y) 

+2G(0,1 - z, 0,1 - z; j/) - G(0,1 - z, 1,0; y) - 4G(0, -z, 1 - z, 1 - z; y) + 2G(0,0,1 - z, 1 - z; y) 
-G(0,1,1 - z, 0; y) - G(0,1,0,1- z; y) - 2G(1,1 - z, 0,0; y) + 2G(1, -z, 1 - z, 0; y) 

+2G(1, -z, 0,1 - z; y) - 2G(1,0,1 - z, 0; y) + 2G(1,0, -z, 1 - z; y) - 4G(1,0; y) - 2G(1, 0,0,1 - z; y) 

+G(1,0, l,0;y) + 2G(1,1,0,0; y) + ^H(l,l;z) - ^H(l;z)G(l - z;?/) + ^G(l - z, 1 - z;?/) 

+ ^H(l;z)G(l,0;y) - ^H(l, 0; z)G(l - z; y) - ^H(1,0, l;z)- ^G(l - z, 1, 0; y) - ^G(l, 1 - z, 0; y) 

-^G(l, 0,1 - z; y)] - ^ - ^H(0; z) + 10H(0; z)G(l - z, 1 - z; y) + 8H(0; z)G(l - z, 1 - z, 0; y) 

-8H(0; z)G(l - z, -z, 1 - z; y) + 4H(0; z)G(l - z; y) + 8H(0; z)G(l - z, 0,1 - z; y) 

-2H(0; z)G(l - z, 0; y) - 4H(0; z)G(l - z, 1,0; y) - 16H(0; z)G(-z, 1 - z, 1 - z; y) 

+5H(0; z)G(-z, 1 - z; y) + 8H(0; z)G(0,1 - z, 1 - z; y) - 2H(0; z)G(0,1 - z; y) 

-4H(0; z)G(l, 1 - z, 0; y) - 4H(0; z)G(l, 0,1 - z; y) + H(0; z)G(l, 0; y) + 3H(0, 0; z) 

+8H(0, 0; z)G(l - z, 1 - z; y) - 8H(0, 0; z)G(l - z; y) - 2H(0, 0,1; z) + 8H(0,0,1; z)G(l - z; y) 

+8H(0, 0,1,1; z) + 10H(0,1; z) + 8H(0,1; z)G(l - z, -z; y) - 5H(0,1; z)G(-z; y) - 3H(0,1; z)G(0; y) 
-4H(0,1; z)G(l, 0; y) - H(0,1,0; z) + 4H(0,1,0; z)G(l - z; y) + 4H(0,1, 0,1; z) - 10H(0,1,1; z) 
+16H(0,1,1; z)G(-z; y) - 8H(0,1,1; z)G(0; y) - 7H(1; z) + 10H(1; z)G(l - z, -z; y) 

+8H(1; z)G(l - z, -z, 0; y) - 6H(1; z)G(l - z; y) + 8H(1; z)G(l - z, 0, -z; y) - 10H(1; z)G(l - z, 0; y) 

-8H(1; z)G(l - z, 0,0; y) - 16H(1; z)G(-z, 1 - z, -z; y) + 20H(1; z)G(-z, 1 - z; y) 

+16H(1; z)G(-z, 1 - z, 0; y) - 32H(1; z)G(-z, -z, 1 - z; y) + 14H(1; z)G(-z; y) 

+16H(1; z)G(-z, 0,1 - z; y) - 5H(1; z)G(-z, 0; y) - 8H(1; z)G(-z, 1,0; y) + 8H(1; z)G(0,1 - z, -z; y) 
-10H(1; z)G(0,1 - z; y) - 8H(1; z)G(0,1 - z, 0; y) + 16H(1; z)G(0, -z, 1 - z; y) - 5H(1; z)G(0, -z; y) 
-4H(1; z)G(0; y) - 8H(1; z)G(0, 0,1 - z; y) + 8H(1; z)G(0,0; y) + 4H(1; z)G(0,1,0; y) 

-8H(1; z)G(l, -z, 0; y) - 8H(1; z)G(l, 0, -z; y) + 5H(1; z)G(l, 0; y) + 8H(1; z)G(l, 0,0; y) + 3H(1, 0; z) 
+8H(1, 0; z)G(l - z, -z; y) - 5H(1,0; z)G(l - z; y) - 4H(1, 0; z)G(l - z, 0; y) 

+8H(1, 0; z)G(-z, 1 - z; y) - 5H(1,0; z)G(-z; y) - 4H(1,0; z)G(0,1 - z; y) + H(l, 0; z)G(0; y) 

-5H(1, 0,1; z) + 8H(1, 0,1; z)G(-z; y) - 4H(1,0,1; z)G(0; y) + 6H(1,1; z) + 32H(1,1; z)G(-z, -z; y) 
-20H(1,1; z)G(-z; y) - 16H(1,1; z)G(-z, 0; y) - 16H(1,1; z)G(0, -z; y) + 10H(1,1; z)G(0; y) 

+8H(1,1; z)G(0, 0; y) + 6G(1 - z, 1 - z; y) + 10G(1 - z, 1 - z, 0; y) + 8G(1 - z, 1 - z, 0, 0; y) 

— 10G(1 — z, —z, 1 — z; y) — 8G(1 — z, —z, 1 — z, 0; y) — 8G(1 — z, —z, 0,1 — z; y) + 7G(1 — z; y) 
+10G(1 - z, 0,1 - z; y) + 8G(1 - z, 0,1 - z, 0; y) - 8G(1 - z, 0, -z, 1 - z; y) + 4G(1 - z, 0; y) 

+8G(1 - z, 0,0,1 - z; y) - 8G(1 - z, 0,0; y) - 4G(1 - z, 0,1,0; y) - 5G(1 - z, 1, 0; y) 

-8G(1 - z, 1,0, 0; y) - 20G(-z, 1 - z, 1 - z; y) - 16G(-z, 1 - z, 1 - z, 0; y) 

+16G(—z, 1 — z, —z, 1 — z; y) — 14G(—z, 1 — z; y) — 16G(—z, 1 — z, 0,1 — z; y) 

+5G(—z, 1 — z, 0; y) + 8G(—z, 1 — z, 1,0; y) + 32G(—z, —z, 1 — z, 1 — z; y) — 16G(—z, 0,1 — z, 1 — z; y) 
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+5G(—z, 0,1 — z; y) + 8 G(—z, 1,1 — z, 0; y) + 8 G(—z, 1,0,1 — z; y) + 10G(0,1 — z, 1 — z; y) 

+8G(0,1 - z, 1 - z, 0; y) - 8G(0,1 - z, -z, 1 - z; y) + 4G(0,1 - z; y) + 8G(0,1 - z, 0,1 - z; y) 

-8G(0,1 - z, 0; y) - 4G(0,1 - z, 1,0; y) - 16G(0, -z, 1 - z, 1 - z; y) + 5G(0, -z, 1 - z; y) 

-^G(0; y) + 8G(0,0,1 - z, 1 - z; y) - 8G(0, 0,1 - z; y) + 3G(0, 0; y) - 4G(0, 1,1-z, 0; y) 

-4G(0,1,0,1 - z; y) + 4G(0,1, 0; y) - 5G(1,1 - z, 0; y) - 8G(1,1 - z, 0,0; y) + 8G(1, -z, 1 - z, 0; y) 
+8G(1, -z, 0,1 - z; y) - 5G(1, 0,1 - z; y) - 8G(1, 0,1-z, 0; y) + 8G(1,0, -z, 1 - z; y) - 7G(1, 0; y) 
-8G(1, 0,0,1 - z; y) + 8G(1,0, 0; y) + 4G(1, 0,1, 0; y) + 8G(1,1,0, 0; y) , (4.28) 

Dn(y,z) = 

1 [H(0; z) + G(0; y)] + [ - H(0; z)G(0; y) - 2G(0,0; y)] + [ - H(0; z) 


-3H(0; z)G(0; y) - G(0; y) - 6G(0, 0; y)] 


3(1-2/) 


[2H(0;z)G(0;y) + 4G(0,0;y)] 


Tjr2 T 

+— [ - 22 - H(0; z) - G(0; 2/)] + ^^ [ - 12H(0; z) - 10H(0; z)G(0; y) - 6H(0; z)G(0, 0; y) 

+3H(0; z)G(l, 0; y) - 10H(0, 0; z) - 6H(0, 0; z)G(0; y) - 3H(0,1, 0; z) - 3H(1, 0; z)G(0; y) 
-6H(1, 0,0; z) - 12G(0; y) - 10G(0, 0; y) + 3G(0,1,0; y) + 6G(1, 0,0; y)] , 

Eii{y,z) = 

[H(0; z)G(l - z, 0; y) - H(0; z)G(-z, 1 - z; y) + H(0; z)G(0,1 - z; y) + 2H(0, 0; z)G(l - z; y) 


(4.29) 


3y 

+H(0,1; z)G(-z; y) + H(l; z)G(-z, 0; y) + H(l; z)G(0, -z; y) + H(l, 0; z)G(-z; y) - G(-z, 1 - z, 0; y) 


-G(-z, 0,1 - z; y) - G(0, -z, 1 - z; y)] + ^ [H(0; z) - 4H(0; z)G(l - z, 0; y) 

+4H(0; z)G(-z, 1 - z; y) - 4H(0; z)G(0,1 - z; y) - 8H(0,0; z)G(l - z; y) - 4H(0,1; z)G(-z; y) 
-4H(1; z)G(-z, 0; y) - 4H(1; z)G(0, -z; y) - 4H(1,0; z)G(-z; y) + 4G(-z, 1 - z, 0; y) 

+4G(—z, 0,1 — z; y) + 4G(0, —z, 1 — z; y) + G(0; y)] + v) + 2G(0,0; y)] 

+ 3(1 [ - H(0; 2)G(0; y) 


-2G(0, 0; y)] + I " ^(0; ^)G(1 - z; y) + H(0,1; z) + H(l; z)G(0; y) + H(l, 0; z) 

-G(l - z, 0; y) - G(0,1 - z; y)] + ^(^^ [ " H(0; z) - G(0; y)] + ^ [8H(0; z) 

+3H(0; z)G(l - z; y) + 4H(0; z)G(l - z, 0; y) - 4H(0; z)G(-z, 1 - z; y) + 4H(0; z)G(0,1 - z; y) 
-2H(0; z)G(l, 0; y) + 4H(0, 0; z)G(l - z; y) + 4H(0, 0,1; z) - 3H(0,1; z) + 4H(0,1; z)G(-z; y) 
+2H(0,1,0; z) + 4H(1; z)G(-z, 0; y) + 4H(1; z)G(0, -z; y) - 3H(1; z)G(0; y) - 4H(1; z)G(0, 0; y) 
-3H(1,0; z) + 4H(1, 0; z)G(-z; y) - 2H(1, 0; z)G(0; y) + 3G(1 - z, 0; y) + 4G(1 - z, 0, 0; y) 
-4G(-z, 1 - z, 0; y) - 4G(-z, 0,1 - z; y) + 3G(0,1 - z; y) + 4G(0,1 - z, 0; y) - 4G(0, -z, 1 - z; y) 

+8G(0; y) + 4G(0,0,1 - z; y) - 2G(0,1, 0; y) - 4G(1,0, 0; y)] + ^ [H(0; z)G( 1 - z; y) 


44 



+2H(0; z)G(l - 0; y) - 2H(0; z)G{-z, l-z;y) + 2H(0; z)G(0, l-z;y)- H(0; z)G(0; y) 

-H(0; z)G(l, 0; y) - H(0, 0; z) + 2H(0, 0; z)G(l -z;y) + 2H(0, 0,1; z) - H(0,1; z) 

+2H(0,1; z)Gi-z-, y) + H(0,1,0; z) + 2H(1; z)G{-z, 0; y) + 2H(1; z)G(0, -z; y) - H(l; z)G(0; y) 

-2H(1; z)G(0, 0; y) - H(l, 0; z) + 2H(1,0; z)G(-z; y) - H(l, 0; z)G(0; y) + G(1 - z, 0; y) 

+2G(1 - z, 0,0; y) - 2G(-z, 1 - z, 0; y) - 2G(-z, 0,1 - z; y) + G(0,1 - z; y) + 2G(0,1 - z, 0; y) 

-2G(0, -z, 1 - z; y) + 2G(0,0,1 - z; y) - G(0, 0; y) - G(0,1,0; y) - 2G(1, 0,0; y)] , (4.30) 

Fii(y,z) = 

^ [27r2 + H(0; z)G(0; y) + H(0,0; z) + G(0, 0; y)] . (4.31) 


5 Conclusions and Outlook 


In this paper, we have presented analytic formulae for the two-loop virtual corrections to the process 7 * ^ qqg, 
which arise from the interference of the two-loop with the tree amplitude and from the self-interference of 
the one-loop amplitude. Together with the contribution from the self-interference of the one-loop amplitudes 
for 7 * ^ ggg these form the full 0{al) corrections to the three-parton subprocess contribution to 

e+e" ^ 3 jets at NNLO. 

Similar results can in principle be obtained for (2 -|- l)-jet production in deep inelastic ep scattering or 
{V + l)-jet production in hadron-hadron collisions. However, the complexity of the cut structure of the 
non-planar graphs together with the rather different domains of convergence of the one- and two-dimensional 
harmonic polylogarithms makes this a non-trivial task, and we defer this to a later paper. 

It must also be kept in mind that these virtual corrections form only part of a full NNLO calculation, 
which also has to include the one-loop corrections to 7 * —> 4 partons mg where one of the partons becomes 
collinear or soft, as well as tree-level processes 7 * ^ 5 partons Effl] with two soft or collinear partons. Only 
after summing all these contributions (and including terms from the renormalization of parton distributions 
for processes with partons in the initial state), do the divergent terms cancel among one another. The 
factorization properties of both the one-loop, one-unresolved-parton contribution p 0 [-|^ and the tree-level, 
two-unresolved-parton contributions || 6 ^-|^ have been studied, but a systematic procedure for isolating the 
infrared singularities has not been established. Although this is still an open and highly non-trivial issue, 
significant progress is anticipated in the near future. 

The remaining finite terms must then be combined into a numerical program implementing the experi¬ 
mental definition of jet observables and event-shape variables. A first calculation involving the above features 
was presented for the case of photon-plus-one-jet final states in electron-positron annihilation in |^,71|, thus 
demonstrating the feasibility of this type of calculations. A prerequisite for such a numerical program is a 
stable and efficient next-to-leading order four-jet program, where the infrared singularities for the one-loop 
7 * —> 4 partons are combined with the tree-level 7 * —> 5 parton with one parton unresolved. Four such 
programs currently exist [f72|-|^, each of which could be used as a starting point for a full 0{al) NNLO 
three-jet program. 
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A One-loop master integrals 


In this appendix, we list the expansions for the one-loop master integrals appearing in and 

These squared amplitudes can be expressed in terms of only two master integrals evaluated at 

d = 4-2e, 


Bub(/) = 



d^k 


{2TrY kP'ik — pY ’ 


(A.l) 


Box(s23,Si3,Si23) = 


-►1 

Pus 


P2 

Pi 


Ps 

f d^ 




(27r)'^ fc2(fc - p2Y{k - P2- PzYik - Pi- P2- PsY 


(A.2) 


Note that in Eqs. (4.13) and (4.14) we have written the one-loop functions /i and /2 in terms of the one-loop 
box integral in d = 6 — 2e. This is straightforwardly related to the box integral in d = 4 — 2e dimensions by 


Box®(s 23, Si3 , S123) = -Box(s 23, S13, S123) 


2 (d — 3 )si2 
2 

si 2 (d- 4) 


(Bub(si3) -I- Bub(s23) - Bub(si23)) 


(A.3) 


Closed expressions for these integrals for symbolic d in terms of T-functions and the 2 F 1 hypergeometric 
function have been known for a long time (see e.g. The bubble integral reads 


^ , .(47rrr(l + e)r2(l-e)^ 

Bub(si 2 ) = -(-S 12 ) 


1 


167r2 r(l-2e) ' " e(l - 2e) 

In the present context, an expansion of the box integral to the second order in e = (4 — d)/2 is required. 

( Sl3 S23 


(A.4) 


o . ^ .(4^)^r(l + e)r2(l-e) ^ 

Box(s23, si3, S 123 ) = -(-S 123 ) 


1 2 ^4 

-2-e ^ ’ V®123’ S 123 J 


with 

k.iAv^z) 

l4.ip{y,z) 

k.i,o{y,z) 

k.i-ii.y,z) 


167 r 2 r(l - 2 e) 


= 2 , 

= -2H(0;z)-2G(0;y) 


E 


i=-2 


+ 0{e^), (A.5) 


(A. 6 ) 

(A.7) 


TT 

y 


2H(0; 2 )G( 0 ; y) + 2H(0, 0; z) + 2H(1,0; z) + 2G(0,0; y) - 2G(1, 0; y) + 

2H(0; z)G(l - z, 0; y) - 2H(0; z)G(0,0; y) - 2H(0,0; z)G(0; y) - 2H(0, 0,0; z) 
-2H(0,1, 0; z) + 2H(1, 0; z)G(l - z\y) - 2H(1, 0; z)G(0; y) - 2H(1,0, 0; z) 

-2H(1,1, 0; z) - 2G(1 - z, 1, 0; y) - 2G(0, 0,0; y) + 2G(0,1,0; y) + 2G(1,0, 0; y)) 

+ y [-H(0; z) - H(l; z) -f G(1 - z;y) - G(0; y)] , 

/ 4 .i,- 2 (y, z) = 2H(0; z)G(l - z, 1 - z, 0; y) - 2H(0; z)G(l - z, 0, 0; y) - 2H(0; z)G(0,1 - z, 0; y) 
+2H(0; z)G(0,0, 0; y) - 2H(0, 0; z)G(l - z, 0; y) + 2H(0, 0; z)G(0, 0; y) 

-h2H(0,0,0; z)G(0; y) -f 2H(0,0,0, 0; z) -f 2H(0,0,1, 0; z) - 2H(0,1,0; z)G(l - z; y) 


(A. 8 ) 


(A.9) 
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+2H(0,1,0; z)G(0; y) + 2H(0,1,0, 0; z) + 2H(0,1,1, 0; z) + 2H(1,0; z)G(l -z,l-z;y) 
-2H(1,0; z)G(l - z, 0; y) - 2H(1,0; z)G(0,1 - z; y) + 2H(1, 0; z)G(0, 0; y) 

-2H(1,0, 0; z)G(l -z;y) + 2H(1,0, 0; z)G(0; y) + 2H(1, 0,0, 0; z) + 2H(1, 0,1, 0; z) 

-2H(1,1, 0; z)G(l -z;y) + 2H(1,1, 0; z)G(0; y) + 2H(1,1,0, 0; z) + 2H(1,1,1, 0; z) 

-2G(1 - z, 1 - z, 1, 0; y) + 2G(1 - z, 0,1,0; y) + 2G(1 - z, 1,0, 0; y) + 2G(0,1 - z, 1, 0; y) 

+2G(0, 0,0,0; y) - 2G(0, 0,1, 0; y) - 2G(0,1,0, 0; y) - 2G(1, 0,0, 0; 2/) + 

+ y [ - H(0; z)G(l -z;y) + H(0; z)G(0; y) + H(0,0; z) + H(0,1; z) - H(l; z)G(l - z; y) 
+H(1; z)G(0; y) + H(l, 0; z) + H(l, 1; z) + G(1 - z, 1 - z; y) - G(1 - z, 0; y) 

-G(0,1-z;z/) + G(0,0;2/)] . (A.IO) 


B Harmonic polylogarithms 

The generalized polylogarithms S„_p(x) of Nielsen turn out to be insufficient for the computation of multi¬ 
scale integrals beyond one loop. To overcome this limitation, one has to extend generalized polylogarithms 
to harmonic poly logarithms 

Harmonic polylogarithms are obtained by the repeated integration of rational factors. If these rational 
factors contain, besides the integration variable, only constants, the resulting functions are one-dimensional 
harmonic polylogarithms (or simply harmonic poly logarithms, HPLs) If the rational factors depend on 

a further variable, one obtains two-dimensional harmonic polylogarithms (2dHPLs) Bil- In the following, 
we define both classes of functions, and summarize their properties. 


B.l One-dimensional harmonic poly logarithms 

The HPLs, introduced in j^, are one-variable functions H(d;x) depending, besides the argument x, on a 
set of indices, grouped for convenience into the vector a, whose components can take one of the three values 
(1,0,—1) and whose number is the weight w of the HPL. More explicitly, the three HPLs with w = 1 are 
defined as 


r dx' 

H(l;x) = / ^^=-ln(l-x), 

H(0; x) = Inx , 

r dx’ 

H(-l;x) = / =ln(l + x); 

Jo J + X 

their derivatives can be written as 


dx 


H(a; x) = f(a; x) 


a = 1 , 0 , —1 , 


where the 3 rational fractions f(a; x) are given by 


f(l;x) = 

1 

1 — X 

f( 0 ;x) = 

1 

f(-l;x) = 

X 

1 

1 -1- X 


(B.l) 


(B.2) 


(B.3) 
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For weight w larger than 1, write a = {a,b), where a is the leftmost component of a and b stands for the 
vector of the remaining {w — 1) components. The harmonic polylogarithms of weight w are then defined as 
follows: if all the w components of a take the value 0 , a is said to take the value Ou, and 


H( 0 uj; x) = —;■ In™ X 
w\ 


while, if a ^ duj, 


In any case the derivatives can be written in the compact form 

dx 


H(a; a;) = / dx' f{a;x') . 

JQ 

n the compact form 
H(a; x) = f(a; a:)H( 6 ; x) , 


(B.4) 

(B.5) 

(B. 6 ) 


where, again, a is the leftmost component of a and b stands for the remaining {w — 1 ) components. 

It is immediate to see, from the very definition Eq. (B_b), that there are 3™ HPLs of weight w, and that 
they are linearly independent. The HPLs are generalizations of Nielsen’s poly logarithms The function 

Sn,p(x), in Nielsen’s notation, is equal to the HPL whose first n indices are all equal to 0 and the remaining 
p indices all equal to I: 

S„,p(x) =H(d„,lp;x) ; (B.7) 

in particular the Euler polylogarithms Li„(x) = S„_ip(x) correspond to 

Li„(x) = H(d„_i, 1;x) . (B. 8 ) 

As shown in the product of two HPLs of a same argument x and weights p, q can be expressed as a 
combination of HPLs of that argument and weight r = p + q, according to the product identity 

H(p;x)H(g;x) = ^ H(f;x), (B.9) 

r—p[iiq 

where p, q stand for the p and q components of the indices of the two HPLs, while p W q represents all mergers 
of p and q into the vector f with r components, in which the relative orders of the elements of p and q are 
preserved. 

The explicit formulae relevant up to weight 4 are 

H(a ; x) H( 6 ; x) = H(a, 6 ; x) + H( 6 , a; x) , 


and 


H(a; x) H( 6 , c; x) = H(a, 6 , c; x) + H( 6 , a, c; x) + H( 6 , c, a; x) , 

H(a; x) H( 6 , c, d; x) = H(a, 5, c, d; x) + H(&, a, c, d; x) + H( 6 , c, a, d; x) + H( 6 , c, d, a; x) 

H(a, 5; x) H(c, d; x) = H(a, 6 , c, d; x) + H(a, c, 5, d; x) + H(a, c, d, 6 ; x) 

+ H(c, a, b, d; x) + H(c, a, d, 5; x) + H(c, d, a, b\ x) , 


(B.IO) 

(B.ll) 

(B.12) 


where a,b,c,d are indices taking any of the values (1,0,—1). The formulae can be easily verified, one at a 
time, by observing that they are true at some specific point (such as x = 0, where all the HPLs vanish except 
in the otherwise trivial case in which all the indices are equal to 0), then taking the x-derivatives of the two 
sides according to Eq. (B.6) and checking that they are equal (using when needed the previously established 
lower-weight formulae). 

Another class of identities is obtained by integrating (B.4) by parts. These integration-by-parts (IBP) 
identities read: 

H(mi,..., TOq; x) = H(mi; x)H(to 2 , ..., TOq; x) — H(m 2 , mi; x)H(m 3 ,..., x) 

-I- ... -I- (-l)'5+^H(m,,..., mi; x) . (B.I3) 

These identities are not fully linearly independent from the product identities. 

A numerical implementation of the HPLs up to weight w = 4 is available [0 . 
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B.2 Two-dimensional harmonic poly logarithms 

The 2dHPLs family is obtained by the repeated integration, in the variable y, of rational factors chosen, in 
any order, from the set 1/y, l/{y—l), l/(j/ + z—1), l/(y+ z), where 2 ; is another independent variable (hence 
the ‘two-dimensional’ in the name). In full generality, let us define the rational factor as 

g(a;y) = ^— , (B.14) 

y-a 


where a is the index, which can depend on z, a = a(z); the rational factors which we consider for the 2dHPLs 
are then 


g(0;2/) 
g(i;y) 
g(l - z-,y) 
g{-z]y) 


1 

y 


1 



1 

y + z-l ’ 
1 

y + z 


(B.15) 


With the above definitions the index takes one of the values 0,1, (1 — z) and {—z). 

Correspondingly, the 2dHPLs at weight ic = 1 (i.e. depending, besides the variable y, on a single further 
argument, or index) are defined to be 


G(0;y) 

G(l;y) 

G(1 - z;y) 
G{-z;y) 


In y, 

In (1 - y) , 




(B.16) 


The 2dHPLs of weight w larger than 1 depend on a set of w indices, which can be grouped into a ic-dimensional 
vector of indices a. By writing the vector as a = (a, b), where a is the leftmost component of a and b stands 
for the vector of the remaining {w — 1) components, the 2dHPLs are then defined as follows: if all the w 
components of a take the value 0, a is written as 0^; and 


G(du,;y) = ^In^y , 

Wl 


(B.17) 


while, if a ^ du,, 

G(a;y)=/ dy'g(a;y') G(6;y') . (B.18) 

Jo 

In any case the derivatives can be written in the compact form 

^G(a; y) = g(a; y)G(6; y) , (B.19) 

dy 

where, again, a is the leftmost component of a and b stands for the remaining {w — 1) components. 

It should be noted that the notation for the 2dHPLs employed here is the notation of , which is different 
from the original definition proposed in [^ . Detailed conversion rules between different notations, as well as 
relations to similar functions in the mathematical literature (hyperlogarithms and multiple polylogarithms) 
can be found in the appendix of Q. 
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Algebra and reduction equations of the 2dHPLs are the same as for the ordinary HPLs. The product of 
two 2dHPLs of a same argument y and weights p, q can be expressed as a combination of 2dHPLs of that 
argument and weight r = p + q, according to the product identity 

G{p;x)G{q]x) = ^ G{f-,x) , (B.20) 


where p,q stand for the p and q components of the indices of the two 2dHPLs, while p W (f represents all 
possible mergers of p and (f into the vector r with r components, in which the relative orders of the elements 
of p an d q a re pres erved. The explicit product identities up to weight w = 4 are identical to those for the 
HPLs ( B.1C| )-( B.12 ), with all H replaced by G. 

The integration-by-parts identities read: 


G{mi,... ,mq]x) = G{mi;x)G{m2,---,mq;x) — G{m2,rni;x)G(rn3,...,mq]x) 

-I- ... -I- ...,mi;x) . 

A numerical implementation of the 2dHPLs up to weight w = 4 is available M- 


(B.21) 
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